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Abstarct 


As considerable portion of the energy required to drive the global atmosphere 
comes from tropical region in the form of precipitational heating, the tropical 
dynamics play a major role in setting the pattern of world-climate. This heating 
and associated dynamics in the tropics exhibit a, wide range of both temporal and 
spatial variabilities. It is imperative to understand these variabilities as it leads to 
a better predictive skill. Over the past few decades much progress ha.s been made 
with models of different complexities to delineate the main physical processes that 
give rise to them. Despite these efforts, the general skill in simulating the detailed 
structure of the tropical variabilities is still limited. One possible reason for this 
may be the difficulty in modeling of convectively induced internal dynamics which 
dominates the tropical atmosphere. 

Several studies have demonstrated the significant roles played by the moist feed- 
backs and internal dynamics in the genesis and structure of the tropical vari- 
abilities. In particular, it was shown that the introduction of convective tim< 
lag (CTL) brings about selective destabilization of the tropical normal modes a.i 
observed intraseasonal frequencies. However, these studies involving CTL hav< 
been done in an analytical setting. The assumption in these studies that trop 
ical variabilities are largely controlled by convectively forced internal dynamics 
needs to be examined in a more realistic non-linear setting, for a wider range c 
observed structures. This thesis examines this question from a number of angle 
and for a number of time scales ranging from monthly to subseasona! to intei 
annual. The primary objective here is to formulate and develop a model of th 



tropical anomaly circulation incorporating convectivelv forced internal dynamics. 
For this purpose, the concept and formulation of CTL has been modified to model 
the non-linear convective forcing and related circulation anomalies in the lower 
troposphere of tropical atmosphere. 

A brief outline of the overall content and the organization of the thesis is given 
below 

CHAPTER I: Introduction 

The tropical atmosphere exhibits a. number of variabilities at widely different 
spatial and temporal scales. Thanks to the large-scale and persistent field ex- 
periments of the last few decades, there now exist a. number of data sets on 
the structure of these variabilities. Using these data sets like C'OADS. it is now 
possible to examine the spatio-temporal structure of the tropical anomaly circu- 
lation for a number of years. These observations impose powerful constraints for 
modeling and provide the only tool for an objective evaluation of various model 
simulations. This introductory chapter describes the salient features of tropical 
variabilities at various scales and discusses their implications for modeling studies. 
CHAPTER II: The Driving Mechanisms of Tropical Variabilities 
Considerable efforts have been devoted in the past few decades to understand 
and simulate the tropical variabilities. In particular, a number of studies have 
explored different feedback mechanisms between convection and large-scale dy- 
namics that can generate the observed structure of tropical variabilities. These 
modeling efforts can be broadly divided into the three categories: process models, 
intermediate models and general circulation models (GCM). While the process 
models help to identify specific physical processes responsible for certain observed 
variability, they cannot address the question of detailed spatio-temporal struc- 
ture of the flow. Another limitation of the simple process models is that, being 
often linear, they cannot accommodate the non-linear and threshold processes. 
The GCMs, on the other hand, do not suffer from such limitations. However, 
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the general skill of the GCM in simulating the detailed structure of the tropical 
variabilities is still limited. The intermediate model stands somewhere in be- 
tween the process models and the GCMs in their complexity and completeness. 
In particular, models that address the question of only anomaly circulation for a 
prescribed (space-time dependent) mean state can be termed intermediate mod- 
els. Thus the model developed and used in this study is an intermediate model 
of tropical anomaly circulation. This chapter discusses the general mechanisms 
of convective forcings in the tropics and their effect on the tropical circulation. 
CHAPTER III: The Model Formulation 

'rids chapter describes the physical basis and the formulation of the model- 
equations. In particular, this chapter describes the concept and the formulation 
of the CTL for parameterizing the effect of organized convection on tropical cir- 
culation. A major portion of this chapter is devoted to the derivation of the 
model-equations. A discussion on the implemented numerical scheme and initial 
and boundary conditions is also included. 

CHAPTER IV: The Monthly Anomaly Fields: Observed vs Simulated 
Structures 

Next, the model’s ability to simulate the observed structure of tropical variabili- 
ties is examined. For this we examine a quantity 


A a^(x,y,m) 




(0.1) 


a Zo( X l ?/) m ) 

where <j£ 0 (x,y,m) is the standard deviation of the observed data for the month 
to (to= 1,12 for January to December). The observed data set consisted of the 
anomaly fields for 12 months for 8 years (1979-1987). The corresponding simu- 
lated quantity, a£ s (x,y,m) represents the standard deviation from eight one-year 
model simulations with initial conditions from respective years. Here, £ denotes 
either zonal or meridional component of the wind. If was found that both CTL 
and the structure of moist feedback (EVVF), arising from the structure of under- 



lying SST, have profound effect on the structure of the variabilities. In particular, 
it was found that a convective timelag of about 0.5 hour generates fields which 
can account for a large part of the observed structure of variabilities. 

CHAPTER V: The Annual Cycle of the Composit Anomaly Circula- 
tion 

We next examine the model simulations for the structure related to annual cy- 
cle. The sixteen-simulation composite model fields were analyzed to examine the 
annual cycle for different locations and dynamical conditions. In particular, the 
results were analyzed for two values of CTL which could generate the structure 
of the monthly anomaly fields in agreement with the observed structures. It is 
shown that the model annual cycle exhibits a large part of the spatial (zonal) 
and temporal characteristics of the observed fields. These findings have impor- 
tant implications for long- and medium-range forecasting in the tropics. 
CHAPTER VI: The Spectrum of Oscillations: The Total Spectrum 
A prominent feature of the tropical circulation is the existence of a number of 
oscillations at different temporal and spatial scales. In the intraseasonal scale the 
most prominent one is the 30-50 day eastward propagating oscillation, while in 
the interannual scale the most well documented oscillation is the ENSO related 
variabilities. Hence, one of the requirements of a good tropical-tropospheric model 
is its ability to generate these variabilities. We have, therefore, examined the 
existence and structure of oscillations at both intraseasonal and interannual time 
scales in our model. In particular, a detailed power spectrum analysis of the model 
simulations for various dynamical conditions as well as for different values of CTL. 
The model results considered for analysis of ISO were generated by compositing 
a number of one-year model simulations for different initial conditions, while 
the existence of longer period oscillations were analyzed using a number of 38- 
year model simulations. We first investigated the total spectra for both annual 
and interannual simulations using red noise spectra for identifying significant 
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peaks. It was found that both the mean conditions (in terms of distribution 
of SST and mean wind) and CTL affect the genesis of' variabilities at different 
periods. A significant finding was that convective forcing ran excite oscillations 
with interannual period even without ocean-atmosphere coupling. A detailed 
analysis of each of the significant oscillations is then presented in subsequent 
chapters. 

CHAPTER VII: The Spectrum of Oscillations: The 30-50 Day Oscil- 
lation 

Perhaps the most prominent and the most well studied intrascasonal oscillation 
in the tropics is the 30-50 day oscillation. Following its identification in the 
early seventies, a large number of studies have examined various aspects of this 
oscillation in detail. Hence, the properties of 30-50 day oscillation provide an 
excellent testing ground for our model. This chapter, therefore, examines the 
existence and various properties of this oscillation in detail using a composite 
one-year model simulation. It is shown that the model can simulate most of 
the established characteristics of the 30-50 day oscillation. Moreover, a few new 
features had been revealed in our analysis, which may be verified in a detailed 
observational study. 

CHAPTER VIII: The Spectrum of Oscillations: The Quasi-Biweekly 
Waves 

Another prominent, but less studied oscillation of the tropical circulation is the 
10-20 day or quasi-biweekly (QBW) oscillation. Observationally, this oscillation 
is most prominent over the summer monsoon region; it has not been generally 
reported over other regions of the globe. The QBW has significant influence on 
the monsoon circulation and precipitation. So the existence and structure of 
QBW in our composite model simulation has been examined. Once again it was 
found that the power spectra of simulated QBW are in close agreement with that 
of observed. Certain new characteristics to be verified observational)' are also 



mentioned. 

CHAPTER IX: The Spectrum of Oscillations:i Interannual Varaibili- 
ties 

After establishing the ability of the model to simulate the structure of the ob- 
served intraseasonal oscillations, we address the issue of longer period oscillations. 
One of the prominent low-frequency oscillations in the tropics is the quasi-biennial 
oscillation (QBO). Due to the recent observational analyses of the structure of 
the QBO, it is now possible to address these questions with our model. In partic- 
ular, we could show that the model simulates a QBO with spatial and temporal 
characteristics in close agreement with observations. 

An important climate variability in the tropics is the El Nino-Southern (ENSO) 
scale variabilities. Although the most widely held view on the genesis of the 
ENSO variabilities is ocean- atmosphere interaction, it is possible that the at- 
mospheric dynamics alone can generate interannual oscillations, although their 
structures are modified by ocean-atmosphere coupling. It was indeed found that 
the CTL induced convective forcing can generate ENSO scale and longer timescale 
variabilities. This chapter also discusses the structure of the interannual variabil- 
ities with reference to the observed structure for different mean and dynamical 
conditions. In addition, we address and investigate the important question of 
inter-relationships among oscillations at annual, biennial and interannual time 
scales. 

CHAPTER X: Convective Couplng and Tropical Variabilities 

It has long been held that the most prominent interannual variability in the trop- 
ics, viz. ENSO, is an ocean-atmosphere coupled variability. However, our findings 
with the effect of SST-induced distribution of moist feedback on the characteris- 
tics of ISO and monthly fields imply that ocean dynamics has considerable impact 
even on shorter time scale variabilities like ISO. This is also supported by recent 
observational studies. To investigate the effect of CTL in a coupled environment, 
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we have formulated ocean-atmosphere coupling involving CTL in terms of convec- 
tive coupling. The main principle of convective coupling is that ocean-atmosphere 
coupling, and in particular SST, can affect the atmospheric heating only by mod- 
ulating atmospheric convection (or column precipitation)-!, lie atmosphere in the 
large scale is insensitive to the underlying SST as a direct heat source. Since' 
many of the observed tropical variabilities exhibit a predominantly Kelvin wave 
character, the implications of convective coupling are investigated for the dynam- 
ics of Kelvin waves in an analytical setting. It was found that, depending on the 
mean conditions, convective coupling with a convective timelag, can selectively 
excite waves over a wide spectrum, from intraseasonal to annual to interannual 
with observed characteristics. This sets the stage for incorporating CTL in a. 
more realistic coupled model. 

CHAPTER XI: Conclusioing Remarks 

The last chapter presents the conclusions and scope for future works. The main 
conclusion that emerges is that a significant part of the tropical variabilities can 
be understood as forced by convection induced dynamics. The success of this 
formulation in the present work indicates that it may be useful in parameterizing 
convection in more complex models. 

The present work shows that SST can significantly affect the tropical variabilities 
even at monthly and other intraseasonal scales. This has been seen in Chapter 
V to IX, where SST effects the circulation implicitely through its effect on the 
strength of the EWF, and also in Chapter XI, where SST is included in a cou- 
pled and dynamical framework. These findings have important consequences for 
medium range weather forecasting in the tropics. 

It is worth noting that the first-baroclinic mode model developed by us can sup- 
port a rich spectrum of physics. This chapter also discuss the certain possible 
improvements in the formulation of moist fee.dba.cks which can improve' the qual- 
ity of simulation of the anomaly fields. 
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Chapter 1 

Introduction 


As humanity faces problems like anthropogenic climate change and depletion of 
the Earth’s protective ozone layer, the need for a careful and integrated long term 
planning of various industrial and developmental activities in order to preserve 
a delicately balanced bio-geo-ecosphere is becoming more acute. Such long-term 
policy forecast requires as accurate and reliable as possible a description of the 
dynamics of the climate system- the role of mathematical models of climate in 
such decision and policy making thus ranges from helpful to enabling to critical. 
Accurate forecast of the weather and rainfall patterns, on the other hand, can 
make significant contribution to a country’s economy and the quality of life of its 
people through better and advance planning of agricultural strategy. This is par- 
ticularly true for countries like India with an agro-based economy. The advent of 
high performance computing now allows numerical integration of complex three- 
dimensional models of the atmospheric general circulation. However, in spite of 
the tremendous progress made in numerical techniques and computing abilities, 
accurate simulation and forecasting of the complex atmospheric processes at user 
specified range and scales still remains a distant goal. This is particularly true in 
the tropics, where a complex interplay of a. number of physical processes makes 
prediction a challenging task. The tropics, however, play a major role also in the 
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global context. As considerable portion of the energy required to drive the global 
atmosphere comes from the tropical region in the form of precipitational heating, 
the tropical dynamics play a major role in setting the pattern of world-climate. 
This heating and associated dynamics in the tropics exhibit a wide range of both 
temporal and spatial variabilities. It is imperative to understand these vari- 
abilities for developing better predictive skill. Over the past few decades much 
progress has been made with models of different complexities to delineate the 
main physical processes that give rise to them. Despite these efforts, the general 
skill in simulating the detailed structure of the tropical variabilities is still limited. 
One possible reason for this difficulty is the dominance of convectively induced 
internal dynamics in the tropical atmosphere. Several studies have demonstrated 
the significant role played by the moist feedbacks and internal dynamics in the 
genesis and structure of the tropical variabilities. These studies have brought out 
the important role moist processes can play in exciting the observed spectrum 
of tropical variabilities (Neelin et al. , 1987; Lau and Shen, 1988; Goswarni and 
Goswami 1992). With appropriate parameterization of the convective heating 
process it was found that several observed intraseasonal oscillations in the tropics 
can be understood as intrinsic, modes of the tropical atmosphere. In particular, 
Goswami and Rao (1993) and Goswami and Mathew (1994) proposed and ex- 
plored the concept of a convective time lag (CTL) in the genesis and structure 
of tropical intraseasonal oscillations. A significant finding was that introduction 
of CTL brings about selective destabilization of the tropical normal modes at 
observed intraseasonal frequencies. 

However, these studies involving CTL have been done in an analytical setting. 
The implicit assumption in these studies, viz. tropical variabilities are largely 
controlled by convectively forced internal dynamics, needs to be investigated in 
a more realistic non-linear setting, for a wider range of observed structures. In 
particular, whether or not the formulation of CTL can be incorporated into a 
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more detailed model of the tropical atmosphere, with relevant, non-linearities 
needs to be examined. Such a model then can provide the basis for an improved 
prediction model for the tropics. This question has been examined from a number 
of angles and for a number of time scales ranging from monthly to subseasonal 
to interannual. The primary objective here is to formulate and develop a tropical 
anomaly circulation model incorporating convectively forced internal dynamics. 
In the present version of the model, we shall focus on a model of intermediate 
complexity. Thus a relatively simpler dynamical situation with more emphasis 
on the formulation of convective processes is considered. In particular, the model 
uses the concept and formulation of CTL to model (non-linear) convective forcing 
in the lower tropospheric anomaly circulation in the tropics. 

A significant contribution to the advancement of understanding and the mod- 
eling of the tropical atmosphere came from the large number of observational 
programmes in the last few decades. Along with these observational programmes, 
the modeling of atmospheric general circulation made enough advance to be able 
to supplement the data generated by the observational programmes through com- 
prehensive and objective analysis. Thanks to these large-scale and persistent field 
experiments combined with model analysis, there now exists a number of data 
sets on the structure of the atmospheric circulation. Using data sets like the 
Comprehensive Ocean Atmospheric Data Sets (GOADS), it is now possible to 
examine the spatio-temporal structure of the tropical anomaly circulation for a 
number of years. Analyses of these data sets enabled the scientists to unravel 
structures hitherto unknown in the tropical atmosphere. A typical example is 
the identification of 3-6 day westward propagating wave with a zonal wavelength 
of a, bout 6000 km. by Liebmann and Hendon (1000) from an analysis of the 7-year 
reanalysed data from European Center for Medium Range Weather Forecasting 
(ECMWF). Another example, involving long-term oscillations, is the discovery of 
quasi-biennial oscillation (QBO) and its interconnections with El Nino-Southern 
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Oscillation (HNSO). As we shall see, the detailed spatio-temporal structures now 
exist for most of the variabilities in the tropics. These observations and analyses 
put powerful constraints on modeling and provide the only tool for an objective 
evaluation of various model simulations. 

Fairly in pace with the observational studies, enormous efforts have been de- 
voted in the past few decades to understand and simulate the tropical variabil- 
ities. In particular, a number of studies have explored different feedback mech- 
anisms between convection and large-scale dynamics that can generate the ob- 
served structure of tropical variabilities. These modeling efforts can be broadly 
divided into the three categories of process models, intermediate models and gen- 
eral circulation models. While the process models help to identify specific physi- 
cal processes responsible for certain observed variability, they cannot address the 
question of detailed spatio-temporal structure of the" flow. Another limitation of 
the simple process models, which are often linear, is that they cannot accom- 
modate the non-linear processes. The General Circulation Models (GCM), on 
the other hand do not suffer from such limitations. However, the general skill of 
the GCM in simulating the detailed structure of the tropical variabilities is still 
limited. 'The intermediate models are somewhere in between the process mod- 
els and the GCM in their complexity and completeness. In particular, models 
that address the question of only anomaly circulation for a prescribed (space- 
time dependent) mean state can be termed intermediate models. Thus the model 
developed and used in this study is an intermediate model of tropical anomaly 
circulation. The organization of this thesis is given below. 

After providing a brief discussion of the various forcings important for the 
variabilities of the tropical atmosphere, we turn towards developing a model for 
the anomaly circulation in the tropics incorporating convection induced internal 
dynamics governed by CTL. The third chapter describes the physical basis of our 
model and the formulation of the model equations. In particular, this chapter 
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describes the concept and the formulation of the CTL for parameterizing the effect 
of organized convection on tropical circulation. A discussion on the numerical 
scheme implemented and initial and boundary conditions is also included. 

The crucial step of model calibration and validation was carried out by com- 
paring model simulation of the monthly anomaly fields with the corresponding 
observed structure. In particular, we consider composite-simulations for a num- 
ber of values of CTL to determine the range of values that would provide the 
best simulation of the observed anomaly fields. In addition, this chapter also 
investigates the effect of different dynamical conditions, like the spatio-temporal 
structure of evaporation- wind feedback (EWF), on the quality of model simu- 
lations. As will be shown, both CTL and the structure of EWF, arising from 
the structure of underlying SST, have profound effect on the structure of the 
variabilities. 

Next, the model simulations are examined for the annual cycle. For this purpose 
sixteen one-year model simulations were generated with initial conditions from 
observed anomaly fields for 1979 to 1987. The simulated fields were then analyzed 
for the nature of variation through the year. A significant finding was that SST- 
induced spatio-temporal distribution of moist feedback has profound impact, on 
the structure of the fields. In particular, spatio-temporal variation of the strength 
of EWF, determined by distribution of SST field, is necessary to simulate the 
correct distribution of the convective heating. These findings have important 
implications for long- and medium-range forecasting in the tropics. 

A prominent feature of the tropical circulation is the existence of a number 
of oscillations at different temporal and spatial scales. In the intraseasonal scale 
the most prominent one is the 30-50 day eastward propagating wave, while in 
the interannual scale the most well documented oscillation is the ENSO related 
variabilities. Hence, one of the requirements of a good model of the tropical 
tropospheric circulation is its ability to generate these variabilities. While it was 



shown in earlier studies that moist feedbacks and CTL can give rise to selective 
excitation of the tropical waves at observed intraseasonal frequencies (Goswami 
and Ran, 1993; Goswami and Mathew 199-1), these studies have been on an 
analytical setting. In particular, these studies could not investigate the effect of 
non-linear convective heating. The existence and structure of oscillations at both 
intraseasonal and interannual timescales in the model simulation are examined. 
We first investigate the total spectra for both annual and interannual simulations 
using red noise spectra for identifying significant peaks in the power spectra. 
It was found that both mean conditions (in terms of distribution of SST and 
mean wind) as well as CTL affect the genesis of variabilities with intraseasonal 
period. A significant finding was that convective forcing can excite oscillations 
with interannual period even without ocean-atmosphere coupling. A detailed 
analysis of each of the significant oscillations is then presented in subsequent 
chapters. 

The available information on the structure of various tropical oscillations can be 
exploited to critically evaluate model simulations against the observed structure 
of these oscillations. Perhaps one of the most prominent and well studied intrasea- 
sona.1 signal in the tropics is the 30-50 day oscillation. Following its identification 
in the early seventies, a large number of studies have examined various aspects 
of this oscillation in detail. The properties of 30-50 day oscillation thus provide 
an excellent testing ground for our model. The existence and various properties 
of this oscillation are examined using a composite one-year model simulation. It 
is shown that the model can simulate most of the established characteristics of 
the 30-50 day oscillation. In addition, model simulations with special dynamical 
conditions like only-winter and only-summer mean winds have helped to identify 
the relative roles of the annual cycle of the mean wind and SST. 

Another prominent, but less studied oscillation of the tropical circulation is the 
10-20 day or quasi-biweekly (QBW) oscillation. Observationally, this oscillation 
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is most prominent over the summer monsoon region. The QBW lias significant 
influence on the monsoon circulation and precipitation. Therefore the existence 
and structure of QBW in our composite model simulation are examined. 

Apart from the intraseasonal oscillations, the tropical climate is also charac- 
terized by a number of long term oscillations. Among these QBO with a period 
between 1 9-36 months and ENSO with a, period between 3-6 years are most im- 
portant. The next chapter, therefore, examines the existence and characteristics 
of these low frequency oscillations in our model simulation. The power spectrum 
analysis of simulated time series at selected locations shows considerable power 
at low frequencies close to that observed, it is also shown that the existence of 
interannual variabilities requires the interannual variabilities in the SST fields. 

Many of observational analyses brought out the interconnection between os- 
cillations at different time scales (Murakami et al. ,' 1986; Philander, 1984) In 
order to successfully simulate the climate, a model not only has to simulate the 
individual oscillations, but also their interconnections. This aspect of the model 
simulations was also investigated. It is shown that the simulated and the observed 
phase relationships among annual cycle, QBO and ENSO are in good agreement. 

It has long been held that the most prominent interannual variability in the 
tropics, viz. ENSO, is an ocean-atmosphere coupled variability. However, O.ur 
findings with the effect of SST-induced distribution of moist feedback on the char- 
acteristics of ISO and monthly fields imply that ocean dynamics has considerable 
impact even on shorter time scale variabilities like ISO. This is also supported 
by observational studies. To investigate the effect of CTL in a coupled environ- 
ment, we have formulated ocean-atmosphere coupling involving CTL in terms of 
convective coupling. The main principle of convective coupling is that ocean- 
atmosphere coupling, and in particular SST, can afFect the atmospheric heating 
only by modulating atmospheric convection (or column precipitation); the atmo- 
sphere on the large scale is insensitive to underlying SST as a direct heat source. 
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Since many of the observed tropical variabilities exhibit, a predominantly Kelvin 
wave character, the implications of convective coupling are investigated for the 
dynamics of Kelvin waves in an analytical setting. It was found that, depending 
on the mean conditions, convective coupling with a CTL, can selectively excite 
waves over a wide spectrum, from intrasea.sona.1 to annual to interannual with 
observed characteristics. This sets the stage for incorporating CTL in a more 
realistic coupled model. 

The last chapter presents the conclusions and scope for future work. The main 
conclusion is that a significant part of the tropical variabilities can be understood 
as forced by convetion induced dynamics. We have characterized this convec- 
tive dynamics through CTL, which represents a timescale for interaction between 
organized convection and large-scale dynamics. The success of this formulation 
in the present work indicates that it may be useful in parameterizing convec- 
tive processes in more complex GCMs. A significant finding from the present 
work is that SST can significantly affect the tropical variabilities even at monthly 
and other intraseasonal scales. These findings have important consequences for 
medium range weather forecasting in the tropics. It is worth noting that the first- 
baroclinic mode model developed by us can support a rich spectrum of physics. 
While there are several directions in which further improvement is possible and 
necessary, the ability of the model to simulate a wide range of observed phe- 
nomena, with time scale ranging from a few days to several years, indicates its 
versatility. 
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Chapter 2 


The Driving Mechanisms of Tropical 

Variabilities 


Large-scale convection and circulation anomalies in the tropical atmosphere have, 
been observed to operate with certain characteristic time scales and spatial pat- 
terns. In particular, a large part of the tropical variability is determined by the so 
called intraseasonal oscillations (ISO). Prominent examples of these ISO are the 
30-50 day eastward propagating oscillation or the Madden and Julian oscillation 
(Madden and Julian, 1971,1972), the 10-20 day westward propagating mode (Kr- 
ishnamurthi cl al., 1973; Murakami and Frydrych 1974; Chang and Chen, 1993) 
and 3-6 day westward propagating wave (Liebmann and Hendon, 1992). A pri- 
mary question that needs to be answered concerns the mechanism(s) that govern 
the excitation and dynamics of these variabilities at the observed spatial and tem- 
poral scales. The significant role played by these variabilities in the regional as 
well as global weather systems explains the persistent effort at simulating them 
through models of varying degrees of complexify. In spite of these efforts, the 
detailed regional and global manifestations of ISO are still evasive to modeling 
studies. However, considerable advances in understanding the mechanism of gen- 
esis of several of the variabilities have taken place in the recent years through 
investigation of roles the of moist processes and the interaction of convection 
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and the large-scale dynamics associated with it. In the longer term scales, much 
attention has been focussed on the variabilities related to El Nino-Southern Os- 
cillation (ENSO) and, on the role of ocean-atmosphere coupling. Here, an outline 
of the various mechanisms that may be operative in the genesis of variabilities in 
the tropical atmosphere is provided. 

Tropical Wave Dynamics and Intraseasonal Oscillations 

hollowing Mastsuno (1966), a large number of the theoretical investigations 
have examined the ISO as manifestation of the internal dynamics of the trop- 
ical atmosphere. The idea that the intraseasonal oscillations could be intrinsic 
unstable modes of that part of the tropical atmosphere situated over the warm 
ocean pool, and not a globally propagating mode, is also strongly suggested by 
observational studies. 

The most distinctive aspect of the 30-50 day oscillation (namely its slow east- 
ward movement, and equatorial emplacement of its maxima) prompted earlier 
investigators to search for an analytical explanation of the phenomenon in terms 
of Kelvin wave dynamics. Parker (1973) interpreted the oscillation as the Kelvin 
wave mode studied by Matsuno (1966), although the vertical scale of low- fre- 
quency Kelvin wave is very small in the absence of heating. Chang (1977) sug- 
gested that the vertical-scale can be explained by the effect of strong viscosity. 

A fundamental problem in invoking the Kelvin wave dyhamics to explain MJO, 
the apparent disparity in the phase speed, was addressed by Swinbank et al. 
(1986). They have come up with some evidences that this may be understood in 
terms of moist Kelvin wave dynamics, if eTect of latent heat release on stability is 
taken into account properly . A simple model-calculation leads to the conclusion 
that precipitation associated with the Kelvin wave reduces static stability and so 
decreases the phase speed. They further confirmed these conclusions by GGM 
experiments in which latent heat-release wa.s increased, first by increasing SST 
and second by doubling the value of latent heat constant. In both the cases the 
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period of wave was increased. 

Another serious drawback with many studies that introduced tropical wave dy- 
namics to explain ISO is the lack of scale selection at observed (intraseasonal) 
frequencies. Only wave with highest growth rate will be observable, since by 
interpretation of linear theory this is the one that draws maximum energy. How- 
ever, in most studies involving Kelvin wave such a selective excitation (or scale 
selection) was not achieved. The highest growth rates still occurred at the high- 
est frequency. This resulted in a search for mechanisms that would provide a 
high-frequency cut-off in agreement with the observed spectra. 

CISK and wave-CISK 

The concept of conditional instability of the second kind (CISK) was first dis- 
cussed by Charney and Eliassen and Murakami, and is a simple representation of 
interaction between small scale convection and large'scale dynamics (and hence 
second kind) in presence of moist processes (and hence conditional). Since its 
formulation almost 30 years ago, CISK has been employed extensively as an ex- 
planation for various tropical phenomena that involve convection. The central 
point is that CISK cannot work in an atmosphere that is not conditionally un- 
stable. 

Recently the commonly held assumption that the mean tropical atmosphere is 
conditionally unstable has come into question. Betts (1982) noted the correspon- 
dence between the observed atmospheric profile and the a virtual moist-adiabatic 
lapse rate. The virtual moist adiabatic takes into account the effect on buoyancy 
of the condensed vapor, whereas a moist adiabatic assumes condensate immedi- 
ately rains out. Xu and Emanuel (1989) analyzed numerous tropical profiles and 
found that the atmosphere did indeed approximate to a virtual moist adiabatic. 
Convective towers probably carry their condensate to at least some level so that 
appropriate adiabatic is a reversible virtual adiabatic rather than an irreversible 
moist adiabatic. 
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if atmosphere is neutrally stable, intensification of disturbance will rely criti- 
cally on the surface fluxes and convection will perform the role of communicating 
the perturbations in the boundary layer moist entropy to the rest of the atmo- 
spheric column. The large scale dynamics would then respond to the rearrange- 
ment of the atmospheric mass that surface fluxes and convection supplied. In 
contrast to the CISK scenario, in this case there is no alternation in the free at- 
mospheric temporal, ure (and hence the pressure distribution) unless the boundary 
layer moist entropy changes. In recognition of this range of interactions, and to 
avoid further confusion from use of the term CISK, Neelin and Yu (1993) have 
suggested the term convective interaction with dynamics (CID) to cover all the 
possibilities. 

The CISK has failed to achieve what it was created for. Because the theory 
predicted highest growth rate at the smallest scale. In other words, it predicted 
that the only structure that should be seen in the tropical atmosphere are the 
cumulus clouds, and not, say the tropical cyclones. A modification of the old 
CISK idea by including the wave dynamics was proposed by Lindzen (1973), and 
came to be known as wave- CISK. The wave-CISK has been often questioned, 
since it invokes physically unreasonable negative heating in a linear setting. This 
happens in the linear theory due to the proportionality of heating to the wave 
divergence irrespective of their sign. 

In several papers, K.M Lau and his co-workers have used a convergent de- 
pendent heating scheme to argue that the intraseasonal oscillation is wave-CISK 
mode that propagate east due to moisture convergence (Lau and Peng 1987; Sui 
and ban 1989; Lau at al., 1989). Lau and Lau (198(5) have analyzed the 12-yea.r 
output of a GCM experiment with the objective of describing the spatial and tem- 
poral behavior of this tropical variability. The GCM used for this integration is a 
9-level. 15- wave number spectral model with realistic distribution of continents 
and orography. The model physics consists of radiative transfer process, ground 
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hydrology, moist, convective adjustment and ice-albedo feedback. In the course 
of integration, SST and solar heating were constrained to evolve through succes- 
sive annual cycles which are identical to each other. Even though some of the 
characteristics of the simulated oscillation ( ie . spatial structure, phase relation- 
ship between the upper and lower level variables, propagation characteristics and 
seasonal dependence) are in agreement with the observations, the period of wave 
number I component of model disturbances is shorter than the corresponding 
observed value. 

Hendon (1988), using a two- level model and a. ChSK-type cumulus heating 
scheme, identified its structure as that of a. nonlinear coupled Rossby- Kelvin 
mode. The mechanism of propagation is essentially one where the large-scale 
equaforially trapped CiSK disturbance forced by the convection neutralizes the 
stability that generates it preferentially to the west of its center, thus giving rise 
to slow eastward propagation. 

Evaporation-Wind Feedback 

Another mechanism that has been proposed to significantly affect the convec- 
tive forcing in the tropics is evaporation-wind feedback (EWF). The evaporation 
anomaly are generally small, but in the delicate balance tropical atmosphere mod- 
est gradient in the evaporation can have significant effects on the evolution of the 
convective activity. The first theories taking into account the effect of surface 
evaporation were proposed by Emanuel (1987) and Neelin et al. (1987) in the 
context of 30-50 day oscillation. In particular, they had considered the dynamics 
of Kelvin wave in an equatorial /3-plane for a moist neutral atmosphere. Consid- 
ering a. case in which the mean background winds are (implicitly) easterlies, they 
showed that inclusion of EWF destabilizes the Kelvin wave. However, a major 
problem was that the growth rate was found to be monotonically increasing with 
frequency, implying selection at the smallest scale. If the wind speed dependence 
of evaporation is changed by replacing its instantaneous value by the c.limato- 
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logical value, the amplitude of the power in the eastward propagating waves is 
substantially reduced. 

Prescribed vs Dynamical Heating 

Lau and Peng (1987) conducted two numerical experiments. The focus was the 
validation of two widely used heating parameterizations against the key features 
of the phenomenon. In one of the experiments, the heating was assumed to be 
localized and pulsating with a. period of 10 days. No feedback with the internal 
dynamics was permitted. Apparently, the horizontal scale of the main circulation 
was determined both by the scale of the heating and the dissipation. The major 
disagreement between the response of the model and reality is that both eastward 
and westward expansion of the wind perturbation are seen in the former, whereas 
only eastward propagation is seen in the latter. 

In their numerical experiment to determine the relative roles of localized heat- 
ing and wave-CISK, Lau and Peng (1987) removed the unphysical negative heat- 
ing by prescribing a positive-only heating. Although in their simulation the east- 
ward propagation was captured, the phase speed was about a factor two faster 
than what is observed. The amplitude was seen to be strongly modulated by the 
variation in SST. The ( MSK-iuduoed responses wen 1 more strong over the warm 
ocean. But the geographical preference of the convection was not seen. Strong 
similarity of the horizontal structure of wind anomalies in the deep tropics to that 
of the Kelvin wave had prompted them to propose so called ryiobile wave-CISK 
mechanism by which Kelvin waves were selectively destabilized. Kelvin wave 
consists of solely divergent winds. Therefore, associated low-level convergence is 
stronger than that of Rossby waves. Hence, heating source moves eastward along 
with the Kelvin wave convergence which in turn again excites Kelvin wave further 
east. 

The effect of vertical structure of heating on the phase speed was investigated 
by Sui and ban (1989). The convection was parameterized by using a simple 
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one-dimensional cloud model which took into account the available moisture sup- 
ply in the lower troposphere and the mean thermodynamic state for the entire 
troposphere. Consequently, the spatial distribution of convective heating could 
be determined internally as a function of SST. The major disturbances were 
seen to propagate eastward with phase speeds of 18 ms -1 and 9 ms" 1 when the 
maximum heating was near 500 mb and between 700-500 mb respectively. The 
corresponding time periods were 20 and 52 days. In all the experiments a uni- 
form SST distribution was assumed. However, for the case when heating was 
related only to boundary layer processes, the value of phase speed (13 ms" 1 ) was 
in between the phase speeds of two wave- CISK modes discussed earlier. But the 
wind response in this case was dominated by wave number 2. It was also stated 
that the boundary layer processes enhance the convection to the east of the deep 
convection thus reinforcing the eastward propagation. 

The oscillations simulated in another study (Lau et al. 1987) with the help of a 
spectral-GCM resemble observations to some extent. This model has same reso- 
lution as the model analyzed by Lau and Lau (1986), but the boundary conditions 
were more idealized. The prescribed SST-distribution was zonally symmetric and 
function of only latitude, convection was modeled by moist convective adjustment 
parameterization scheme. The annual mean solar input was prescribed at the top 
of the boundary and clouds were prescribed as a function of latitude only. There 
was a sharp contrast between the simulated dynamical and precipitational fields; 
larger portion of the total variance in the former resided in the wave number l 
and 2, whereas in the latter it resided in the higher inodes. Moreover, the small 
scab' precipitation did not propagate steadily at the same phase speed, like those 
presented in the models where Kuo type scheme was used to calculate healing. 
Instead, in the course of their eastward journey, they redeveloped in the regions 
favored by large-scale convection. The simulated circulation patterns at both 200 
and 850 mb had strong zonal components. There was a much stronger meridional 
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component at 850 mb, which presumably produced the frictional stress near the 
surface. 

Teleinduction and Genesis of Super Cloud Clusters 

According to Winston the super cloud clusters and their eastward movement 
originate from what they called as the supercloud clusters teleinduction mech- 
anism of initiating new cloud clusters by the existing cloud clusters at some 
distance (800-1200 km) to its east. Locations close 1 to the cloud clusters, because 
of the compensating downward motion associated with the existing cloud clus- 
ters, are not favored for the emergence of the new cloud cluster. The reason the 
eastern side is more favored than the western side for new cloud cluster to emerge 
is due to the fact that the basic flow in the boundary layer, being easterlies, is 
strengthened by the circulation induced by the existing cloud cluster. Once new 
cloud is formed it compete with existing cloud clusters for moisture. Since mois- 
ture mainly comes from the east in tropical easterly flow, new cloud has the first 
chance to harness it. He argues that the initiation of the cloud is not caused 
by the Wavc-CISK, but by the circulation features induced by the existing cloud 
clusters 800-1200 km to the west of it. The wave- CIS K is an instability respon- 
sible 1 for the growth of the cloud clusters, and cloud cluster teleinduction is the 
mechanism that brings the dynamical and thermodynamical fields at a particular 
location towards the critical point beyond which the wave instability starts. 
Extra-tropical Influence 

The assumption that each new intraseasonal episode of convective activity is 
initiated by the remnants of the previous oscillation approaching from the west 
has been the underlying assumption in many theoretical studies. There are also 
studies which report the potentiality of midlatitude baroclinic eddies in trigger- 
ing the convective activity. Blade and Hartmann (1993) examine these possible 
mechanisms responsible for the selection of the preferred period of MJO. For this 
purpose, a. global two-layer nonlinear model with a positive-only CISK type eu- 
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niuhis healing scheme is used. By making the SST in the tropics as a function 
of longitude, they investigate the relative importance of these two mechanisms in 
determining the time scale of the oscillation. When the cold SST is sufficiently 
stable, the CISK-wave propagates efficiently through the stable region in the form 
of damped moist Kelvin wave, and reemerges in the unstable region where its am- 
plitude grows again. When SST in the stable sector is set close to the instability 
threshold, the moist Kelvin waves (MKW) slow down and decay before reenter- 
ing the unstable region. However, the GISK perturbation periodically regenerates 
over the warm water in response to the local buildup of instability. This shows 
that MKW in this case is a forced response, and has no impact on the genesis of 
the next convection. This implies a new mechanism for setting the timescale' of 
the oscillation, alternative to that of simple zonal propagation around the globe. 
A recharge- cl targe theory is proposed whereby, the time scale of the disturbances 
is set by the growth and duration times of the convective episode together with 
the recharging time of the static stability in the convective region. Also shown 
is how midlatitude baroclinic eddies provide the quasi-stochastic force necessary 
to excite each new intraseasonal episode by organizing a region of subtropical 
con vecliou, which then grow and expand towards equator due to the effect of 
latent heating. 

Ocean-Land- Atmosphere Interaction 

The question, as to whether atmospheric variabilities at intraseasonal and 
shorter time scales are entirely forced by internal atmospheric instability or 
whether coupling between ocean and atmosphere play a role is still not fully 
answered. Most of the linear theories assume the availability of steady source of 
energy for the maintenance of the basic state. The source of energy for the sus- 
tained cumulus heating (on the subseason al time scale) is usually not addressed. 
Convection for such long period requires a steady supply of moisture from the 
ocean. It is found that in this timescale the latent heat flux from the ocean to 
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atmosphere is of considerable amount. The study by Krislmainurthi et al. (1988) 
indicate that oscillation of both the wind and the SST contribute to a strong 
coupling of atmosphere and the ocean at this time scale. The amplitude of SST 
variation at t his time scab' are of the order of 0.5" to 1.0"C and most, pronounced 
SST oscillation in t he SST occurs over the equatorial western pacific ocean and 
the Bay ol Bengal. If is conceivable that, one degree enhancement of SST over 
t hese climatological warm oceans can ailed the convection and their propagation. 

Oscillations of intraseasonal time scale have been observed in the tropical 
oceans. For example Mysak and Mcrtz (1984) reported a 40-60 day oscillation of 
long shore current and sea surface temperature in the western Indian Ocean. It 
is worth noting that 40-50 day signals have been found in many oceanographic- 
parameters, e.g., surface current over the Indian Ocean (McPhaden, 1982); sea 
level in the coastal regions of South America (Enfield and Lukas, 1984) and sea 
surface temperature in the Somali current (Mysak and Mertz, 1984). The dy- 
namical model study by Webster and Chou (1980) hints at the possibility that 
this intraseasonal variability could be an outcome of interaction of the ocean or 
the ground hydrological cycle with atmosphere. 

At annual, interannual and longer time scales the importance of ocean-atmosphere 
interaction has been recognized for a long time. A conceptual model for the 
KNSO based on sustained ocean-atmosphere interaction was proposed by Bjerk- 


ness (1969). A significant finding was that a coupled equatorial ocean-atmosphere 
system can sustain unstable oscillations, although the linear equatorial ocean and 
the atmosphere may by themselves be stable. In the conventional explanation, 
ENSO can be understood as a manifestation of equatorial Kelvin wave dynamics 
and associated delayed oscillator mechanism. In this scenario (Cane and Zebaik, 
1985; Schopf and Suarez, 1988; Battisli and Hirst, 1989; Philander, 1992), oceanic- 
wave dynamics and the reflection of the wave from boundary determine the struc- 
ture and growth rate of coupled instability. Hence, it was argued that ENSO type 
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oscillations should he absent over the tropical Atlantic and Indian oceans. Never- 
theless, several observational studies suggest that ENSO related variabilities are 
present also outside the tropical Pacific. A global relation of the surface pressure 
anomalies to the ENSO cycles was demonstrated by Krishna, murthi r.t al. (198(1). 
Barnett (1985) presented the evidence of eastward propagation of surface wind 
and pressure anomalies from equatorial Indian Ocean to Central Pacific. Yasunari 
(1987) suggested that ENSO should be considered as a global phenomena, in a 
land-ocean-at mosphere coupled environment. 

The relationship between these features observed in the ocean and their coun- 
terparts in the atmosphere remains to be ascertained. Thus, coupled atmospheric 
model may necessary for a complete understanding of this complex global phe- 
nomena. 'Phe question of role of ocean- atmosphere coupling in the genesis of 
sub-seasonal scale variabilities will be taken up in more detail in the present 
study. 

From an analysis of the various mechanisms presented above, and their relative 
success in explaining the genesis of tropical variabilities, the most universal one 
appears to be convection induced internal dynamics through interaction between 
convection and large scale dynamics. However, the scope and the representation 
of interaction between convection and large scale dynamics is rather broad, en- 
compassing mechanisms like CISK and EWF. Indeed, as described below, the 
current, most successful representation of interaction between convection and 
huge scale dynamics is a, combination of (MSK, EWF and convergence feedback 
in presence of an explicit moisture equation. 

Interaction between Convection and large-scale Dynamics 

Observations indicate that deep convection often occurs over the regions of large 
scale moisture convergence. It is, therefore, reasonable to assume that evapora- 
tion only regulates the amount of moisture available for the precipitation. Seager 
and Zebia.k's study (1991) is related to the interaction between convection and 
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dynamics in a. new linear primitive equation model wherein a new method of 
integration had been employed. An idealized sea surface temperature distribu- 
tion. representing the climatological warm pool area with a maximum at 10Q°E 
and a. cold region with maximum at HOW 0 , is used when simulations were per- 
formed for two different types of heating schemes. The simulated behaviors were 
changed radically when heating was evaluated by Betts-Miller convection scheme. 
Contrary to what is common in wave ClSK-sludies, but in agreement with the 
observations, eastward propagating disturbances do not remain coherent around 
the globe. They are seen only over the warm-pool regions and die out over the 
cold region. The phase speed of convection in this case is less than one third of 
that in the wave-CISK case. If the evaporation wind feedback were removed, then 
a moisture convergent mode is seen with a single disturbance that would encircle 
the globe in a, bout (i() days. If the ('fleet of moisture convergence is suppressed, 
then almost all deep convection and propagation is suppressed and the weak 
modes p rest'! it have only half the amplitude of those present in the experiment 
with full moisture' ('((nation. The model characteristics wore not altered much 
when the dominance of moisture' con vergemce' over the', evaporation wind feedback 
was reduced by assuming a higher value of the surface exchange eoeflicient. This 
is not true when one of the feedbacks is absent. No eastward propagation is ap- 
parent in the presence of only moisture convergence. The conclusion is that for 
different values of the surface evaporation time scale, either the EWF or low-level 
convergence of moisture can create eastward propagating deep convective mode. 
The periodicity is longer when full moisture equation is used. 

A successful application of the formulation of interaction between moist feed- 
back induced convection and large scale dynamics to explain the genesis and 
structure of ISO was made by Goswami and Goswami (1992). In particular, it 
was shown that the combined effect of EWF, convergence feedback and dissipa- 
tion can sele'ctively destabilize the mixed Rossby-gravity wave at 3-6 day time 
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scale. The spectral characteristics of this selectively destabilized wave matched 
very well with the observed 3-7 day westward propagating wave (Liebmann and 
liendcm 1991 ). However, the formulation of Goswami and Goswami, (1992) could 
not account for the scale selection at two other major ISO observed in the tropics 
viz. the 30-30 day oscillation and the 10-20 day oscillation. 

Moisture Dynamics and Convective Time Lag 

As mentioned earlier, other major drawback of most of the studies involving 
CID to explain ISO is their inability to produce selective excitation of the tropical 
wave with observed spectral characteristics. Besides, often only a particular ISO is 
addressed. If one lias to build up a consistent model of anomaly circulation in the 
tropics, a mechanism that can generate under appropriate dynamical conditions 
all the major oscillations in the tropics is needed. A reason for this failure might 
be due to the diagnostic status assigned to moisture variable. In particular, most 
of these studies consider a balance equation for the moisture variable so that the 
procipitational lira, ting is determined by a balance between surface evaporation, 
column convergence and column precipitation. Such a formulation is not likely to 
be valid in a general situation, where the storage of moisture will determine the 
scab's at which the heating will be most effective. In other words, while the moist 
feedbacks can bring about the essential destabilization, the dynamics of moisture 
variable determines the selective destabilization. The role of an explicit moisture 
dynamics in the excitation of tropical waves was formulated and examined by 
Kmamiel ( I OSS) a.nd Goswami and Kao ( 1991). Goswami and Kan, by introducing 
the concept of a convective time lag in an explicit moisture equation, showed that 
convection induced internal dynamics can selectively destabilize Kelvin wave with 
observed spectral characteristics. Subsequently, Goswami and Mathew (1994) 
investigated the dynamics of equatorial waves in presence of an explicit moisture 
equation with a convetive time lag. Goswami and Mathew (1994) showed that the 
3-1 day westward propagating wave ( Liebmann and Mathew, 1991) and the 10-20 
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day westward propagating wave observed over the summer monsoon region could 
bo understood as selectively destabilized tropical normal modes. However, the 
total anomaly circulation in the tropical lower troposphere contains a much wide 
range of variabilities than the these ISO. The question arises therefore, whether 
the analytical model of G os warn i and Rao (1994) and Goswami and Mathew 
( 1991) can be generalized to provide an acceptable representation of the low-level 
variabilities in the tropics. In the chapters that follow this question is examined 
in detail. 
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Chapter 3 

The Model Formulation 


In the background of the large number of analytical studies and simple process 
models, it can be expected that convective internal dynamics determines a large 
part of the anomaly circulation in the tropics. This convective internal dynamics, 
of course, would depend on the embedding mean state. Since the mean state 
embodies the effects of various large-scale forcings as well as boundary conditions, 
proper inclusion of the mean condition and convective internal dynamics should 
enable us to simulate large part of tropica! circulation. This is the working 
hypothesis for developing our model of anomaly circulation in the tropics. 

The spatial similarity of apparently disconnected modes like ISO and ENSO 
of tropical atmosphere (Lau and Chan, 1988; Wieckmann, 1983; Wang and Mu- 
rakami, 1988) is rather intriguing. So is the phase-locked nature of the variabilities 
at the low frequency end of the spectrum. It has been prbposed (Lau and Chan, 
1988) that ISO might trigger the low frequency mode like ENSO. There might 
exist a common dynamical link which could selectively cascade energy from one 
mode to other. In view of this, it is essential to identify the connection between 
various components of these variabilities. This is possible only if all these vari- 
abilities are brought under a common dynamical framework. There have been a 
few theoretical and modeling (Lau and Slien, 1988; Hirst and Lau, 1990) attempts 
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made in this direction. An alternative view of interaction of convection and dy- 
namics in the tropics is proposed to develop a. unified dynamical framework for 
the oscillations at different time scales. 

A characteristic of the geophysical fluid is their large horizontal dimension 
compared to vertical extent. So they behave like shallow fluids. The shallowness 
can be quantified by aspect ratio, 5 = DjL. D and L are typical vertical and 
horizontal scales of the atmosphere respectively. For a synoptic situation in the 
tropics with horizontal scale of 1000 km. and vertical scale comparable to the 
height of tropopause, the aspect ratio is the order of 10" 2 . This has profound 
effect on the vertical dynamics of the geophysical fluids and has led to a class of 
models known as shallow water model. 

3.1 The Shallow Water Approximation 

Perhaps Laplace was the first to apply shallovwwater equations on a rotating 
sphere to study oscillations of the atmosphere. Since then this model had been 
used with appropriate modification to study the dynamical processes of relevance 
to both ocean and atmosphere. Gill (1980) was able to capture the wind and 
pressure distribution patterns similar to those observed in the tropics when the 
model was (breed by a. prescribed heating. 

3.2 The Basic Equations 

The horizontal components of momentum equation can be written as, 


du 

7k 


~ _ du du _ du 

+ Jx + u— + v— + iv— + 
dx dy dz 


1 dp 
p dx 


+ F x 


dv 

dt 


~ _ dv _ dv . dv 1 dp - 

+ f y + u— h V— h io— f- —— h F y 

dx dy dz p dy 


= 0 

= 0 . 


(3.1) 

(3-2) 


where u, f\ and w denote x, y and z components of velocity respectively, whereas 
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p, the pressure and p the density of air. The x and y components of momentum 
dissipation are given by and /'), respectively. The terms f r and /,, are the 
com|)onents of Coriolis force. The vertical momentum equation reduces to a 
simple hydrostatic balance under shallow water approximation. 

The nonlinearity in the advective term seems to play a secondary role (Gill 
and Phillips, 1986) in many of the dynamical scenarios in the tropics. Zebaik 
(1990) evaluated the terms in the vorticily budget for selected monthly surface 
wind anomaly fields over the tropical Pacific Ocean. He arrived at the conclusion 
that nonlinear advection terms are of second order importance. So the equations 
are linearized using conventional method wherein the total fields are partitioned 
into a mean term and a perturbation term. Exploiting the smallness of the 
perturbation terms, the equations 3.1 and 3.2 are simplified as 


du 

In 

dv 


. rr du du 1 dp 
Byv + U— + V— + --f + F x = 0 

ox ay p ox 


. r dv , , dv 1 dp _ 

at +li « u + i ai + '% + + F ’ 


0 


(3.3) 

(3.4) 


The mean zonal wind (U), meridional wind (V) and vertical velocity (IF) are 
characterized by U — V = V(x,y,l) and IV 7 = 0. In the case of p, 

the Houftsinscsq approximation is evoked wherein if is treated as a constant in 
horizontal momentum equations and as a function of z in the vertical momentum 
equation. As the terms involving meridional and zonal gradients of mean wind 
are in general small, they are neglected. The zonal and meridional components 
of anomaly wind are denoted by u and v. The F x and F y stand for x and y 
components of horizontal friction. The Coriolis parameters in an equatorial j3~ 
plane can be approximated as —flyv and flyu respectively. 

The energy equation is given by ! 


£ 

dt 


= Q 


(3.5) 
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Here Q is the source term of total energy with Q as the perturbation heating. Lin- 
oariza! ion of (‘((nation about a. basic stab* having mean potential temperature 
(#(-)) gives: 


do Tr d0 jr de dd ^ 

— + U — b V — f W— b Fg 

of ox oy dz 


Q, 


( 3 . 6 ) 


The horizontal gradients of mean potential temperature are neglected due to its 
small value. The cooling term is represented by Fg. At this stage a spatially and 

(IT) . 

temporally uniform value of — is assumed. 

dz 


3.2.1 Incompressibility and Hydrostatic Balance 


Although in general the atmosphere is compressible, it can be shown that for large 
scab' dynamics incompressibility is a reasonable approximation. The continuity 
equation is, then, linearized to obtain 


du dv dw 

JiT + dTj + 7h 


(3.7) 


The observation that scale of motion in vertical plane is very small compared to 
that in the horizontal plane helps to simplify the vertical component of equation 
to a hydrostatic relation which can be written as 


+ P9 — 0? (3*8) 

where g is the acceleration due to gravity. 

A simple linear relationship between perturbation potential temperature and 
pressure used in earlier studies is adopted here. The details of the derivation and 
assumptions which lead to this simplification are given be!ow r . 

The equation of state can be written as 
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p + p 

p 


ip + nW + T) 

PT 


(3.9) 


Afirr sonic cancellation and neglecting the product, pT, equation 3.9 gives 


P T 

~ + rp =0 , (3.10) 

p J 

when we assume, for simplicity, that the fractional changes in temperature are 
much greater than those in pressure. 

The Poisson’s equation may be linearized to obtain 


o _ T_n_p 

0 - T ~ Cpp 


(3.11) 


The last term in the right hand side of this equation is dropped, based on the 
above assumption. Now, hydrostatic relation (3.8) may be written as 


<Kp + P) 


{p + p) 9z 
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The equations 3.11 is used to modify equation 3.13. 


1 dp 0 

Jidz 9 9 
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We may write the equation 3.11 as, 


(3.12) 


(3.13) 
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dl 
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(3.15) 



1* s i rl her approximat ion is possible only if oscillat ions having vertical wavelengths 
less t han half of (.lie scale height, of density are considered. Under this condition 
the magnitude of the second term is negligible. Now we may rewrite equation 
•5.15 as 


3 

Tz 




= 0 


3.2.2 The Vertical Structure 


(3.16) 


In the present formalism the tropical atmosphere is. represented by the shallow 
water equations describing the horizontal structure of the first baroclinic mode on 
an equatorial /3-plane. Though a. number of modes may be required to tackle gen- 
eral problems, a remarkable degree of success had been achieved in understanding 
certain aspects of tropical atmosphere behavior by using even single moth' (Kr- 
ishnamurthi el al., 1976; Anderson, 1976 and liannon 1979). llantel and Baader 
(1978) indicated that, when suitably averaged in horizontal and time the heating 
rate increases from zero at the surface to a maximum at about 500 mb, and then 
fall off above that. The wind in lower troposphere is oppositely directed to that 
in the upper troposphere. The horizontal velocity is zero in mid-level, whereas 
the vertical velocity is maximum there. This characteristics of first baroclinic 
mode can be mathematically represented as the following vertical structures of 
the model fields: 


u = u n (x,y,t)cos(z/H) 
= v a (x,y,t)co$(z/I-I) 
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(3.17) 

(3.18) 
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( 3 . 19 ) 

( 3 . 20 ) 


= iv a (n,yJ)sm(z/H) 

(O.Q) = (0 a ,Q a )sin(z/l!)t 

"o 

p = p a co.s{z/H )^r (3.21) 

Po 

We have neglected higher-order modes, a.s the first baroclinic mode provides a. 
reasonable approximation for tropics. This vertical structure has been used in 
several studies (l)avey, 1985; Lau and Shen, 1988; Davey and Gill, 1987; Goswami 
and Goswami, 1991). The 0 a and J> 0 arc' diarac.teristir values of 0 and ]>. This 
structure is bounded by rigid horizontal plane at 2 = nil and 2 = 0. Substituting 
3.19 and 3.21 in 3.10 


Pa 

Po 


= -q iPj- 

0o 


The perturbation vertical velocity (to) expressed as 


(3.22) 


10 a = ~H 


dUg 8Va 
dr dy 


(3.23) 


is obtained by removing vertical structure from the continuity equation by using 


3.17, 3.18 and 3.19. 

3.2.3 Moisture Dynamics 


Although the importance of moist dynamics in the tropical variabilities has been 
emphasized in many works, most of them do not consider explicit moisture dy- 
namics. In these studies, unlike in the present case, convective heating is deter- 
mined by a balance between column convergence and surface evaporation. 

The equations 3.3, 3.1 and 3.6 describing the evolution of perturbations of zonal 
and meridional wind components at the lower layer and the mid level potential 
temperature perturbation 0 ri respectively are supplemented by a fourth equation 
describing the evolution of depth integrated moisture: 


29 



()(j f ()u a <)i 


()i 


+ 


Ox (hj 


= E - P 


(3.2-1) 


Hern is the gain due to evaporation from the surface and P is the loss due to 
precipitation in the column and 


file do* 

V = —H / f- f(z)dz. (3.25) 

J 0 f/~ 

In the above expression, 1I 0 represents a height relevant for moisture dynam- 
ics. r I'he unit of q* is kgm~ 3 . The f(z) is a non-dimensional structure function 
representing the vertical structure of the dynamical fields. Since our model is 
formulated for the first baroclinic structure, we consider H 0 = H. For simplicity, 

-p is assumed to be independent of height and its value is estimated by assuming 
dz 

q*(H) = 0, and a climatological value for q*(z = 0). In the tropics moisture con- 
centrations tends to fall off with height over a vertical distance small compared 
with that of the first baroclinic mode. It should be noted that the explicit time- 
dependence of the moisture variable in our study allows us to avoid assuming 
a permanently convective region assumed in many studies ( c.q ban and Shell, 
1 988; (loswami and (losvvami 1991). In our formalism a non-convective (convec- 
tive) region can evolve to a convective (non-convective) region by the dynamics 
of moisture variable'. 

3.2.4 The Convective Time Lag and Precipitation 

The models of tropical dynamics are sensitive to how convection is handled. Most 
frequently used theories for interpreting intraseasonal oscillations are based on 
Wave- (1 IS K (ban and Peng, 1987; Miyahara, 1987; Hendon, 1988; Chang and 
bim, 1988; Wang, 1988; Kirtman and Vernekar, 1993; and Blade and Hartmann, 
1993) and evaporat ion wind feedback mechanisms (Neelin, 1987; Emanuel, 1987). 
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In majority of those works precipitation was estimated from a steady stale mois- 
ture budget equation. The linear versions of those theories suffer from scab' 
selection problem. Moreover, Wavc-CJISK mechanism often produces eastward 
propagating wave's with higher phase speed compared to the observed one. The 
EWF implicitly depends upon the presence of mean easterlies. In fact, wind is ob- 
served to be weak easterly or even westerly in locations (if. Indian and western 
Pacific oceans) where intraseasonal signal is strong. The studies by Numagati and 
Ilayashi (1991) and Numagati (1993) found out that the convective adjustment 
scheme produces more realistic simulation of convective complexes and its prop- 
agation characteristics. Though the parameterization of convective heating used 
in the present study can be classified as convective adjustment scheme, conceptu- 
ally it is different. In the present formulation, precipitation brings the moisture 
variable q to its equilibrium value q over a time scale r; thus, it can be expressed 
as 



In other words, the moisture variable does not relax back to its equilibrium value q 
immediately but, does so over a moisture relaxation timescale r. Several studies 
have emphasized the role and importance of a convective time-lag in tropical 
circulation (Betts, 1986; Betts and Miller, 1986; Emanuel, 1993; Goswami and 
Rao, 1994; Goswami and Mathew, 1994; Neelin et al, 1994). The concept and 
the effect of a moisture relaxation time scale were also discussed in GIL The 
basic philosophy is that for organized large-scale precipitation in the tropics, the 
process cannot be considered instantaneous; it has a characteristic relaxation time 
scale which can be typically a few hours. Seager et al . , (1995) have also used a 
time-lag of a few hours, although in a different context. 

It should be noted that the GTE embodies our hypothesis and it is not possible 



to assign a precise value to it at present. However, we shall assume that the CTL 
time scale of a few horns. It is obvious from the above equation that as model 
moisture evolve in time, regions of negative precipitation may appear. In the 
numerical version of this model, this limitation has been overcome by setting the 
value of precipitation to zero wherever it is less than zero. 

3.2.5 The Evaporation Wind Feedback 

Using well known bulk aerodynamic formula, the perturbation evaporation 
may be expressed as 


E = p s CdW{q s - q 0 ), (3.27) 

where p s is air density in the boundary layer, Cj the moisture llux coellicient, 
W the surface wind speed, and q s the saturation specific humidity at sea surface. 
The q a represents specific humidity at the surface, which is determined by sea 
surface temperature SST (in units of K) via following empirical relation: 

q 0 = 10~ 3 [18 -f 0.67(55T - 300)]. (3.28) 

This empirical relation has been used by Wang and Li (1993) and is derived from 
monthly mean Comprehensive Ocean-Atmosphere Data Set (GOADS) for the 
tropics where SST is higher than 295 K. This expression for evaporation can be 
further simplified. For this purpose the Claitsius-Clape} ron expiation is integrated 
from a reference temperature T r to SST, which yielded 


q.i{ SST) = q s {T r )<\T 


R V T? 


T r ) 


where /?,, is gas constant for water vapor. Inserting 3.28 and 3.29 into 3.27 we 



E = Pa C d I< q W(SST-l\), (3.30) 

where coefficient A 7 = 6.95 x 10~ 4 A —1 and T, = 293.2 K. This equation implies 
1 hat the evaporation rate increases with SST. Because the empirical relation was 
used, strictly speaking, equation 3.30 can be applied only to the tropical regions 
where SS I exceeds 22'T '. The evaporation rate is assumed to be independent of 
SST and is given by \.8CnWK q , when SST is below 293.2 K. 

A closed sel of model equations constructed by using equations 3.3. 3.1 3.6, 
3.22, 3.23, 3.25, 3.24 and 3.26 is given below. 
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E-P (.1.31) 


4 hc P'u, E v and b'o are expressed as liu a , Rv a and R0 respectively 
The atmospheric heating proportional to precipitation is in the form 

Q a = yP; (3.35) 

i 'A y 

where, ?/ = - — . 

(~ p l I pa 

'The L v and C v denote latent heat of vaporization of water and specific heat of 
air respectively. While density of the air is represented by p a , the value of the 
tropopause height is given 7 r times H . 


3.3 The Basic State 

In this study, the incorporation of effects of advection by mean wind and the effect 
of mean SST through evaporation would enable us to investigate the model’s 
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response to temporal and spatial inhoniogeneity present in the mean stale of the 
tropical af mosphere. The climatological monthly fields for meridional and zonal 
rompoiKMiis ol the wind and SS I’ were derived from monthly mean OOALhS from 
the year !!)7!) to 1987. I lie anomalies for individual year* had boon calculated 
by deducting (9-year) climatological mean from the total field of the each year. 

I hose anomalies form the initial fields for the control runs. The moan SST fields 
also decide the temporal and spatial variation EYYF. 

3.4 Numerical Scheme 

Owing to the presence of non-linearities in our equations it is not easy to obtain 
analytical solutions in general. Besides, we would like to investigate the detailed 
spatial and temporal evolution of the fields in a general setting, lienee, the 
model equations have to be numerically integrated with a suitable scheme. This 
has been done by using explicit finite difference methods in both space and time. 
Symbolically this can be written as, 

Tp 1 - + 2 A (3-30) 

where represents any one of the dependent variables , At the time step and $ 
the generalized forcing term. The time level and spatial grid point are identified 
by the subscript n and superscript k respectively. The numerical integration is 
initiated by using the forward difference scheme. This is inevitable, as Leap- Frog 
scheme can bo utilized only if values of dependent variable are supplied at all time 
levels. The spurious mode inherent in this time-integration scheme is eliminated 
by using Asseline (I 2 1 in lime) filter (Shapiro, 1970) at 1 every 18th time step. 
The model leapfrogs with time step of 500 sec. To avoid two-grid interval noise, 
the variables are staggered on the Arakawa-C grid. The same resolution is used in 
longitude and latitude: ie., 2°. While boundary conditions prescribed in the zonal 



direction arc periodic, the fields are forced exponentially to zero in meridional 
directions I lirough spongy layers implemented over 5 points near the boundaries. 
The integration domain in the meridional direction is ?)(F on either side of* the 
equator. In our formalism, the moisture variable ran assume any positive or 
negative values allowed by the equation d.2 I. But there are physical hounds on 
the magnitude' of tin* moisture' anomaly. Since' a positive q would imply a rise 
above (near) saturated reference state, allowing very large value would make it 
unphysical supersaturation. On the other hand, permitting too large' a negative 
value would make' the total moisture content very small or even negative, which 
is unphysical. We, therefore, impose this bound through the condition 

q- if q < q- 
q+ if q > q+ 

Although the lower and the upper bounds can in principle be different, we have 
taken them to be symmetrical about zero. Besides, it was found that a symmet- 
rical bound with q+ = </_ = 5 kgm~ 2 produces the daily rates of precipitation 
which were in best agreement with observation. 
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Chapter 4 


The Monthly Anomaly Fields: 
Observed vs Simulated Structure 


The anomaly fields in the tropics show a high degree of spatio-temporal vari- 
abilities. They also exhibit characteristic organization at various spatio-temporal 
scales. While there has been a large number of attempts to simulate specific vari- 
abilities like 30-50 clay oscillations, it is a prerequisite of a model of atmospheric 
circulation to be able to simulate the structure of the anomaly fields. The dy- 
namics of the anomaly field (as against the structure of a specific oscillation) are 
subjected to a number of forcings determined by geophysical conditions and mean 
state. Moreover, tropics are never dynamically isolated from the extra- tropics, 
especially at large spatio-temporal scales. However, the major forcing for the 
anomaly circulation in the tropics is still likely to be the convective forcing. The 
convective forcing itself is not independent of the mean conditions; for example, 
the convective activity seems to closely linked to the underlying SST fields. It 
is reasonable to expect that a proper representation of convective forcing along 
with the effect of mean state should explain a fair portion of tropical variabili- 
ties. To examine this hypothesis and to validate our model, a comparison of the 
spatio-temporal structures of the simulated vs the observed fields is made. The 
convective forcing in our model is determined by precipitational heating which 


36 


Ac for x=0.5 hrs. with uniform EWF 



Longitude (°) 


Longitude (°) 


Longitude-latitude structure of the Aa„ field with uniform EWF for 


Figure 


depends on large-scale convergence, EWF and CTL. However, many of the pa- 
rameters (in particular CTL) cannot be determined precisely. In addition, we 
need to identify and assess relative roles of various moist processes and CTL in 
the quality of our model simulation. However, the anomaly fields are also influ- 
enced by the mean circulation and underlying SST distribution. Finally, there is 
the question ol the sensitivity of model to initial conditions, ihus there are at 
least three ma jor issues that need to be addressed: 


Role of CTL 


2. The effect of SST-induced distribution of EWF 



x* / ▼ ] 




tiki* 


3. Sensitivity to initial conditions 
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Figure 1.2: Same as the Fig 4.1, but for r 


The first of these issues, viz. the role of GTL, is also related to the model 
calibration. While in earlier analytical studies a CTL of about six to twelve 
hours was used, these studies were aimed at investigating mechanisms of specific 
low-frequency variabilities like 30-50 day oscillation and 3-5 day oscillation. As 
mentioned earlier, the basic idea, had been that .specific ISO arises from convective 
forcing organized at specific scales. However, while considering the dynamics of 
the anomaly fields the possibility remains that even smaller period organizations 
also contribute to the genesis and the dynamics of the variabilities. Indeed, since 
the 1 arger scale organizations are dependent on the organization at smaller scales, 
it is possible that the smallest organizational time scale determines the subsequent 
hierarchy of organization. In this scenario, the relevant convective timelag should 
be comparable to the cumulus time scale. Since we still do not have a way of 
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igure 4.3: Same as the Fig 4.1, but for 


3.0 hr 


precisely determining fix' value of CTL, a. mnnlx'r of values for CTL, from flu 


cumulus scale to the meso-scale, is considered. A comparison of quality of the 
simulations for the various values of CTL then provides a way of choosing an 
appropriate value of CTL. The second issue, ie. the role of mean conditions, 


is significant in view of the fact that many of our model processes like EWF 
depend (implicitly as well as explicitly) on the mean conditions. For example, the 
quantity Aq appearing in for EWF depends on the underlying SST through the 
dependence of surface humidity ou SST. This dependence can introduce complex 
spatio-temporal variations in the strength of EWF. As we shall see, this affects 
the structure of the variabilities to a large extend. In addition, the presence 
of a space-time dependent advective mean wind has considerable effect on the 
evolution and structure of the variabilities. 
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Figure 4.5: Same as the Fig 4.1, but for A 


180 and sy mmetric about the ecjuator. 
on a 2 ° x 2° grid, with a time step of 500 
into the following two sets: 


All the simulations were also carried out 
sec. The experiments can be categorized 


Setl: Uniform strength of EWF, mean wind prescribed from monthly-mean 
wind for the period 1979-1987 for each month. 


2. Setll: Non-uniform strength of EWF determined by the distribution of 
monthly-mean SST; mean wind prescribed from monthly-mean wind for 
the period 1979-1987 for each month. j 
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i* igtire 4.6: Same as the Fig 4.5. but for 


4.2 Method of Analysis 


LU dLe U1C momm ^ ailomal y fields, the longitude- latitude structures of the 
simulated Melds wee written every five days and then averaged over a month. 

The simulated structure of standard deviation is obtained by compositing model 
simulations lor eael, month from 16 one-year model integrations. Out of these 
10 model simulations, eight were with initial conditions taken from the observed 

anomaly fields for the month of January and the other eight wil l conditions 

derived Iron, the observed anomaly fields for the mould, of for the respective 

yrar ' *" Uk ' ™ M, “' ljoll ‘ U» observed and the composite simulated fields were 
noi mali/cd to then respective maximum, so that any error due to difference in 

is fill '- r " d -n.» allows us to disenss exclusively the rpurlity of 
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Figure 4.7: Same as the Fig 4.5, but for r = 3.0 hr 


l lie simulation in terms of spatio-temporal distribution of the variabilities. Since 
the magnitudes of the simulated fields depend on the strengths of various model 
parameters, many of which cannot be determined precisely, we do not address 
the issue ol magnitudes of the simulated fields. However, for reasonable values 

ol the model parameters the distributions of the simulated fields are dose to the 
observed. 


To assess the quality of simulation wc define and use the normalized magnitude 
ol the difference between the observed and simulated standard deviations given 





Latitude (°) Latitude (°) Latitude (°) 


Aa v for t=6.0 hrs. with uniform EWF 
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Figure 4.8: Same as the Fig 4.5, but for r = 6.0 hr 


where, y, in) represents the standard deviation of the observed monthly 

anomaly field £ (for the period 1979-1987) at the point (x,y) for the month m. 
The corresponding simulated quantity is represented by a^ f (x. y, in). It should be 
noted that since the difference between the observed and the simulated standard 
deviations is divided by the local value of the observed standard deviation, the 
range of A a lies between zero and infinity. While Aa^(x,y, z) provides a. quanti- 
tative and objective measure of the ability of the model to simulate the tropical 
variabilities, an important question concerns the ability of the model to respond 
to different initial condition. This aspect is analyzed by comparing model simu- 
lations for different initial conditions. 







Latitude (°) Latitude (°) Latitude (°) 


for x=0.5 hrs. with SST-dependent EWF 
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Figure 4.9: Longitude-latitude structure of the Acr ?J field with SST-dependent EWF 
for r = 0.. p ) hr. 

4.3 Results 

For t'arli of t lu' sols of experiments I and IF a composite-simulation is constructed 
from I (> simulations corn's ponding fo winter and summer initial conditions for the 
right, year period extending from 1979-1987. In the lol lowing we shall present the 
longitude-latitude structure for the six months, three winter months of December, 
January, February and three summer months of June, July and August. For each 
of the fields u and v the results are presented for the four values of ( !TL: (.1.5, 1.5, 
5.0 and ().() hours. 

Simulation with uniform EWF: Figures from 1.1 fo 1.1 show the longitude- 
latitude structure of the difference between standard deviations (A a u ) of the 
normalized observed and composite simulated u-fields for the period 1979-1987. 

a n 
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Aa u for 1=1 .5 hrs. with SST-dependent EWF 
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Figure 4.10: Same as the Fig 4.9, but for r = 1.5 hr. 


A constant value of EWF corresponding to a. spatially and temporally averaged 
mean condition is used. As indicated in the figures, the results are shown for 
four values of CTL, from 0.5 hours to 6 hours. The corresponding results for the 
meridional wind are shown in figures 4. 5-4.8. In all these Figures, the darkest shade 
represents difference between the observed and the simulated variabilities equal 
to or greater than 100% of the local observed variability. A careful examination 
of the figures from 1.1 to 4.8 reveals that over most of the tropics, and especially’' 
in the summer months, the error is greater than 100%. It is intriguing that the 
simulations in the winter months are better than those in the summer mouths. A 
similar observation holds also for the v. This perhaps indicates that our param- 
eterization of convective forcing is not adequate for strongly convert ive summer 
months. In particular, it appears necessary to take into account the structure of 
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A o u for i=3.0 hrs. with SST-dependent EWF 
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Figure 4.11: Same as the Fig 4.9, but for r = 3.0 hr. 

the underlying SST field and its effect on the strength and distribution of moist 
processes. 

Simulation with SST-dependent EWF: As discussed in context of our model 
formulation. SST can implicitly affect moist processes in our system at least in 
two ways: through its ('(feet on the amount and distribution of the mean column 
integrated moisture and through its effect on the magnitude and distribution of 
A q, which controls the strength of EWF. Although it is difficult to evaluate the 
relative importance of these two effects a. prior, we shall focus on the second effect 
because a large number of studies have shown the significant role of EWF in the 
genesis and dynamics of tropical oscillations. 

The longitude-latitude structure of the A a u field for winter and summer 
months for four values of (JTL are shown in figures 1.9 to 1.12. As can be 
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Ao v for t=0.5 hrs. with SSI-Dependent EWF 
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Figure 4.13: Same as the Fig 4.9, but for A a v 
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between the observed and the simulated dynamical lields should arise from an 
actual year to year variability in the simulated fields. This ability to respond 
to different initial conditions is essential for the model of tropical climate. This 
issue is Jience, addressed by comparing simulations for two different years with 
corresponding initial conditions. For this purpose, two simulations corresponding 
to 1979 (normal year) and 1982 (an El Nino year) had been chosen. The initial 
condition for each one-year simulation is given a.s monthly anomaly fields for the 
month of January of the respective years. So far we have focussed only on the 
dynamical fields. The very important issue of the simulation of the precipitation 
field has not been addressed. In this section we address these two issues together. 
However, in view of the fact that it is only an SST-dependent EWF and a value 
of r = 0.5 hr. gives rise to the best simulation of the dynamical fields, only these 
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igure 4.14: Same as the Fig 4.13, but for 


simulations are analysed. Figures 4. IT and 4.18 compare the simulated structure 
of monthly precipitation fields for 1979 with corresponding observed structures 
for the same year. The results for the year 1982 are presented in figures 1.19 and 


1.20. In each of these cases only positive part of the observed anomaly is con- 
sidered, as t he simulated precipitation is only positive. In each of the figures the 
left panels are for simulated precipitation while corresponding observed monthly 
anomaly precipitation (positive part) is shown in the right panels. Once again, 
only the distribution of the fields are emphasized. Hence the figures are pre- 
pared according to their respective maximum and minimum and no comparison 


of magnitudes should be made 


A distinct aspect of the simulated precipitation fields, as can be seen from the 
figures 4.17 to 4.20, is their clear shift from the southern to northern latitudes 
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Figure I .LG: Same as the Fig 4.13, but for r = G.O hr 


This hints at the need for a more detailed and sophisticated implementation of 
the spatio-temporal distribution of convective forcing. 


4.4 Conclusions 


r 1 k ' main conclusion that emerges from the above experiments is that convection 


-induced interna.! (anomaly) dynamics can explain a. large pari of the tropical vari- 
abilities. This is particularly true when SST-inducetl spatio-temporal variations 


of the strength of KWh' is included. As we have em 


asr/e< 


between the observed and the simulated (composite) variabilities in this cas e is 


less than 50% of the local observed variability over most of the tropical belt, and 


especially over the equatorial bolt of 1 5 (, S to I5 (, N. The larger error in the higher 
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Figure 1.17: Comparison of simulated precipitation (left panels) with anomaly precip- 
itation (right panels) for the summer months of 1979. Only the spatial distributions 
for each of the months are considered, as the maximum’ and the minimum values and 
contour interval are different for different cases. 

latitudes is not altogether surprising since our model with its /7-plane approxi- 
mation is not expected to be valid at higher latitudes. There are regions even in 
the equatorial belt where the model fails (with error more than 100%) in a sys- 
tematic way. In particular, these regions are essentially continental. Obviously, 
our model with its all-ocean lower boundary, fails to capture the variabilities con- 
trolled by land processes. The fact that this happens systematically raises hope 
lor improving performance by including ground hydrology. 

Another interesting fact arises from these experiments is that for both the 
representations of strength of FAVF (Rxpt Sol, I and Set II), the simulations show 
less errors for the summer hemispheres than for the winter hemispheres. This 
points to the fact that convection plays a greater role in the structure of tropical 
variabilities in the summer than in the winter. Since our model emphasizes 
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Figure 4.18: Same as figure 4.17, but for winter months 

convection-induced dynamics of the tropical variabilities, it fails to account for 
the winter-time forcing which would account for (a large part of the) remaining 
variability. One way of improving our model would be, therefore, to identify and 
parameterize this component of the forcing. From a comparison of the results from 
vSet I and Set II, it is clear that the spatio-temporal variations in the strength of 
EWF induced by mean conditions should be taken into account . In the present 
case, we have represented this dependence of EWF on mean conditions through 
its implicit dependence on SST according to equation 3.29. Considering the 
profound effect of the distribution of EWF on the quality of simulation, it is very 
important to properly formulate the representation of EWF in terms of the mean 
conditions. Also, it will be interesting to carry out the above experiments with 
mean conditions derived from a GCM in place of monthly mean data as used in 
the present case. 
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Simulated and Observed anomaly precipitations for summer months (1982) 
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Figure 4.19: Same as figure 4.17, but for 1982 

One purpose of the present study has been to investigate the effect of the value 
of the CTL in the quality of the model simulations. As has been emphasized ear- 
lier, the CTL embodies the underlying hypothesis of our dynamical model, and 
a precise value of CTL cannot be assigned. Although earlier analytical studies 
of Gil and GM had used a CTL of a few hours, these studies were aimed at 
the mechanisms of particular intraseasonal oscillations. Following this argument 
outlined above, we have considered four values of CTL, from half an hour, corre- 
sponding to cumulus organization scale to six hours, corresponding to meso-scale 
convective systems. Indeed, as can be seen from the results of Set 1 and Set 11, a 
much smaller value of CTL less than 6 hours is needed for capturing the detailed 
structure of monthly anomaly fields. It would appear that a more appropriate 
value of CTL lies around an hour, or so. 

An important conclusion that arises from the present study is that SSI 1 can 
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dgure 1.20: Same as figure 4.17, but for winter months 1982 


play a significant role in the structure of tropical variabilities even at monthly and 
shorter time scales. This is evident from the fact, that a correct distribution of 
the convective heating over the equator is simulated only when the SST-induced 
spatio-temporal distribution of KWT is taken into account.. Although SS I may 
not ailed t he tropica! dynamics at monthly or shorter time scales as a direct heat 
source, it can ailed the atmospheric convective heating. However, for this it is 
necessary t hat the parameterization of convect ion is sensitive to SS I distribution. 

Another success of the present model is its ability to capture the year-to-year 

0 

variabilities in the observed anomaly fields, in particular precipitation. Thus, for 
example the observed anomaly precipitation field for 1979 August, shows a much 
stronger held in the northern hemisphere than the corresponding field in 1982. 
Those differences in the observed fields are present in the simulated fields to a 
largo extent. However, the simulated fields do not possess as much variability 



as in the observed lie ' 1 < I s both temporally as well as spatially. In view of the' 
rrucia.l roll' a yod by 1 1 u ' spat io-temporal structure ol moss! processes further 
ini provemeni in the formulation of moisi processes is likely to Ik* necessary to 
bridg(* the gap bel ween observation and simulaf ion. 

The present study shows the convection induced internal dynamics, governed 
by a. CTL, is cpiite successful in accounting for a large part of the tropical variabil- 
ities in terms of their monthly structure. However, the tropical atmosphere also 
(exhibits a number of characteristic variabilities at intraseasonai and longer time 
scales. Since these oscillations are an important feature of tropical circulation, a 
good model of tropica! variabilities should be able to support these oscillations. 
We investigate this aspect in the next chapter. 



Anomaly Circulation 


! lie anomaly eirenla.l ion in the tropics exhibits a pronounced annual eyrie wit h 
distinct zonal characteristics. The annual cycle's of convection in tin' tropical 
Indian Ocean and Pacific Ocean follow their own characteristic seasonal march. 
The moan annual cycle's present in the wind and SST can significantly affect 
the anomaly circulation. In our model with moist feedbacks, and its implicit 
dependence on mean SS I , the annual cycle in SST can also affect the anomaly 
fields. 

A characteristic feature of the mean annual cycle in the tropics is the association 
of the convective maximum with the Intertropicai Convergence Zone (ITCZ). 
This convective maximum is often seen to move from north to south and hack 
following the annual cycle of solar heating. However, this annual march also 
exhibits pronounced region ality, especially over the tropical Indian Ocean and the 
Pacific, due to the effects of land-ocean distribution and local ocean-atmosphere 
interaction. Thus, over the tropical Indian Ocean- Australia area , the convective 
maximum moves not only meridionally but also zonallv (Lati and Chain 1983). 
Since the dynamics of the -anomaly circulation in our model arc significant I y 
affected by the mean wind and SST, we expert the anomaly circulation also to 



Time-Latitude structure at 90°E: t=0.5 hrs. 
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Figure 5.1; Time-latitude structure of the unfiltered model fields. The negative 
contours are dashed 

exhibit a characteristic annual cycle. Since there is no external forcing in our 
model related to annual cycle, the model annual cycle should closely match the 
mean annual cycle. Therefore, the model annual cycle is first compared with the 
mean annual cycle. The question of an annual cycle intrinsic to the convection 
induced anomaly circulation will be taken up subsequently. 

Analysis of eight-year (1974-1982) monthly mean Outgoing Longwave Radia- 
tion (OLR), obtained from monthly means of daily averages on a 2.5 X 2.5 grid, 
by IVIcehl (1987) reveals the several interesting features about the spatial struc- 
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Time-Latitude structure at 160°W: x=0.5 hrs. 
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Figure 5.2: Same as the figure 5.1, but for 1G0°W 


tures of the annual distributions of variables like OLR, clouds and precipitation. 
Although in our model there is no explicit computation of OLR, the changes in 
OLR are often used to study implied relative fluctuations in the deep convection 
and associated rainfall in the tropics (Liebmann and Hartmann, 1982; Lau and 
Chan, 1983; Weickmann et al, 1985). Hence, the observed structures of both 
OLR and precipitation are compart'd with our simulated precipitation fields. 

To examine the model simulations for the annual structure of the variabilities 
the composite one-year model simulations generated with initial conditions from 
observed anomaly fields for 1979 to 1987 are used. Once again, the simulated 


60 











Time-Latitude structure at 80°W: x=0.5 hrs. 
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Figure 5.3: Same as the figure 5.1, but for 80"W 


annual cycles for different mean and dynamical conditions arc analyzed to identify 
and assess the relative roles of various processes in the genesis and structure of 
the anomaly annual cycle. The issues of the structure of the annual cycle specific 
to certain well-known oscillations will be taken up in the later chapters. 

As we have seen earlier, only a value of CTL between 0.5-1. 5 hours can generate 
the observed structure of the variabilities with reasonable success. To determine 
the role of the value of CTL on the structure of the simulated annual cycle, we 
have examined the four composite one-year simulations, one for each value of 
the CTL. It is found that only r =0.5 or 1.5 hours gives rise to the observed 
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Time-Latitude structure at 90°E: t=1.5 hrs. 
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Figure 5.4: Same as the figure 5.1, but for r = 1.5 hr. 


characteristic. Therefore, the results presented here are only for these two values 
of CTL. In addition, as we have found, a uniform strength of EWF, implying 
a uniform SST, does not simulate the structure of the variabilities. In partic- 
ular, SST dependent EWF is necessary to simulate the correct distribution of 
the convective heating. Hence, only cases of non-uniform strengths of EWF are 
considered to examine its role in the anomaly annual cycle. 
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Time-Latitude structure at t60°W: t=1 .5 hrs. 
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Figure 5.5: Same as the figure 5.4, but for 160°W 

5.1 Time-Latitude Structure 


The annual cycles in the form of time-latitude diagrams at three zonal locations 
are shown in figures 5.1 to 5.3 for r = 0.5 hours. The corresponding results for 
r =1.5 hours are shown in figures 5.4 to 5.6. A comparison of the results for 
the two values of OTL shows that although there are certain differences between 
the two simulations, the essential characteristics are very similar. However, these 
differences are more marked over the off-equatorial areas and especially over the 
northern tropics. 
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Time-Latitude structure at 80 °W: t=1 .5 hrs. 






Figure 5.6: Same as the figure 5.4, but for 80°W 


A remarkable feature of the time-latitude structure for both the values of CTL 
is a distinct zonal character of the annual cycle. In particular, the precipitation 
field over the monsoon region shows a clear shift from the Southern Hemisphere 
to the Northern Hemisphere and back with the march of the season. Such a 
cross-equatorial movement is much less pronounced over the Central and the 
Eastern Pacific. From an analysis of eight-year monthly mean OLR, Meehl (1987) 
showed that the dominant convective activity in the northern latitudes in the 
tropics is found during the monsoon. The monthly anomaly (zonal means covering 
45 to 155°E and 155°E to 90° W; annual average zonal mean subtracted from 
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Time-Longitude structure at 10°N: x=0.5 hrs. 
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Figure 5.7: Time-longitude structure of the unfiltered model fields. The negative 
contours are dashed 

the zonal mean for each month) fields of OLR and cloud also show a relative 
convective maximum north of the equator in July-August and south of the equator 
in January-April. Thus the monthly anomaly fields also tend to follow a cycle 
similar to the annual course of the sun. Such a trend is clear in the simulated 
annual cycles of the zonally averaged fields over the monsoon sector centered at 
90°F (Figures 5.1 and 5.4) and also over the Pacific Ocean centered at 160°W 
(Figures 5.2 and 5.5). 
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Time-Longitude structure at 0°: t=0.5 hrs. 





Figure 5.8: Same as the figure 5.7, but over the equator 


5.2 Time-Longitude Structure 

The similarities in the simulations with the two values of CTL can be seen also 
in the time-longitude diagrams for the unfiltered anomaly fields for one year 
period. The time-longitude diagrams for the model fields for r = 0.5 hours for 
three latitudinal locations, equator, 10°N and 10°S, are shown in figures 5.7 to 5.9, 
respectively. The corresponding results for r = 1.5 hours are shown in figures 5.10 
to 5.12 respectively. 

For both the values of CTL, there is a strong annual cycle over 10°N between 
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Time-Longitude structure at iO°S: t=0.5 hrs. 
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Figure 5.9: Same the as figure 5.7, but for 10°S 

40-120° and 240-300° (in all x-t plots the longitudes continuously increase from 
left to right), marked by strong westerlies and southerlies in the summer months. 
However, over other longitudinal positions the annual cycle is rather weak. In 
particular, there is no reversal of wind. We also note that the annual cycle of 
zonal and meridional winds do not quite coincide, although they overlap, in terms 
of their longitudinal characteristics. The annual cycle of both u and v at 10°S 
for r =0.5 hours (figure 5.9 are confine to the zonal location between 60-240°, 
a characteristic also shared by precipitation. However, evaporation is present 
throughout the year, although with temporal variation in .strength. The precipi- 


67 



Time-Longitude structure at 1 0°N: x=1 .5 hrs. 
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Figure 5.10: Same the as figure 5.7, but for r = 1.5 hr. 


tation seems to occur only when the strength of evaporation reaches a threshold 
value. The former is always accompanied by strong southerlies (northerlies) in 
northern (southern) hemespheric summer (figures 5.7, 5.9, 5.10 and 5.12). This 
clearly brings out the role of cross-equatorial wind in the anomaly annual cycles. 
In the 10°N the precipitation, is seen only with westerlies, but in the 10°S it is 
present in easterlies also. This behavior is common to both values of r used. 


f 
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Time-Longitude structure at 0°: t=1 .5 hrs. 
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Figure 5.11: Same the as the figure 5.10, but over the equator 

5.3 The Role of the Mean Annual Cycle 

An important question at this point is how much of the annual cycle of the 
anomaly circulation is induced by the annual cycle of the mean fields. For 
example, the monsoonal region is characterized by the well-known seasonal rever- 
sal of wind, which may induce an annual cycle also in the anomaly circulation. 
Similarly, the annual march of the SST, through its influence on EWF, can signif- 
icantly affect the anomaly annual cycle. On the other hand, it is possible that the 
convection induced internal dynamics generate an annual cycle that, is intrinsic to 
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Time-Longitude structure at 10°S: t=1.5 hrs. 
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Figure 5.12: Same the as figure 5.10, but for 10°S 

the tropical dynamics. To examine these questions we analyze the annual cycles 
in the one-year model simulations with the annual cycles of mean wind and SST 
suppressed. We shall consider the time-latitude structure of the model fields at 
three longitudes (ie. 90°E, 80°W and 160°W ) to examine the effects of the mean 
annual cycle on the anomaly circulation. 

5.3.1 Role of the Annual Cycle of the Mean Wind 

To study the role of the annual cycle of the mean wind on the anomaly annual 
cycle, we have analyzed the annual cycles from the two one-year simulations 
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Time-Latitude structure at 90°E 



Figure 5.13: Time-latitude structures for summer mean wind (first two panels from 
the top) and winter mean wind (third and fourth panels).- Negative contours are 
dashed. 

with winter and summer mean winds. Thus in the first experiment (of summer 
mean winds) the mean winds were prescribed by averaging the 0 years of (1979- 
87) monthly mean winds in the summer months from June to August, while in 
the second experiment (winter mean winds) the mean winds were obtained by 
averaging the 9 years of climatological monthly mean winds in the winter months 
from December t.o February. In both of these experiments the annual cycle of 
SST was kept intact. In particular, the northward march of high SST with season 
could still affect the model dynamics through its effect on the strength of FAY F. 
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Time-Latitude structure at 160°W 



Figure 5.14: Same as the figure 5.13, but for 160°W 


Figures from 5.13 to 5.15 show the time-latitude structure of the model zonal 
wind and precipitation for the cases of only-summer-mean-wind and the only- 
winter-mean-wind for three zonal locations respectively. It can be seen that the 
pronounced annual cycle over the monsoon sector disappears in the absence of the 
annual cycle of the mean wind. This is especially true in case of the winter mean 
winds, which shows a.n equator centered structure throughout the year. However, 
in case of the only-mean-summer-wind, there is a weak annual cycle, although 
the pronounced north-south shift seen over the monsoon sector with the mean 
annual cycle is no longer present. We note that in presence of the summer mean 
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Time-Latitude structure at 80°W 
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Figure 5.15: Same as the figure 5.13, but for 80°W 

winds the model dynamics respond to the location of higher SST in the northern 
latitudes in the summer months, while no such seasonal change appears in case 
of the only-winter-mean-wind. Another distinction between the model fields with 
summer mean winds and with the winter mean winds is that the fields with 
summer mean winds tend to exhibit a higher degree of spatio-temporal variability, 
with shorter periods and larger north-south asymmetry. Indeed, in the case of 
the only-winter-mean-wind, the fields are quite periodic, with a periodicity of 
about two months. This issue of the difference in the periodicities will be taken 
up when I, ho characteristics of the infraseasonal oscillations in these cases are 
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Time-Latitude structure at 80°W 
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Figure 5.16: Time-latitude structures for summer SST (first two panels from the 
top) and winter SST (third and fourth panels). Negative contours are dashed. 


discussed. Another difference between the u-field from the two simulations over 
90" E is the appearance of the westerly anomalies at 16°N during the summer 
months persisting into the winter. In fact, there are westerlies centered at the 
equator all throughout the year in the case of the winter mean winds. However, 
in the case of summer mean winds these equator-centered westerly anomalies 
appear only during the first three months. A comparison of the u-field with 
only-summer-mean-wind with the corresponding structure with the annual cycle 
of both mean wind and SST show that the main difference between two cases is 
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Time-Latitude structure at 160°W 



Figure 5.17: Same as the figure 5.16, but for 160°W 


the absence of the equatorward shift in the post monsoon months in the absence 
of the annual cycle of the mean wind. Thus the north-south migration of the 
westerlies with season may appear to be primarily due to the annual cycle of the 
mean wind. However, this conclusion is not true, as we shall see when we examine 
the structure in the absence of the annual cycle of the mean SST. Both the zonal 
wind and precipitation are localized too far north in the summer months in the 
absence of the mean wind. 

In contrast to the monsoon region, the annual cycle of the mean wind does 
not seem to influence flic latitudinal distribution of either the zonal wind or 
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Time-Latitude structure at 80 °W 



Figure 5.18: Same as the figure 5.16, but for 80°W 


precipitation over the Eastern or Central Pacific. However, it seems to have 
much greater influence on the temporal distribution. 

5.3.2 Role of the Annual Cycle of SST 

The relevance and significance of SST-dependent EWF in the context of simu- 
lating the structure of anomaly fields have been noted earlier. Now we examine 
the role of the mean SST annual cycle in the annual cycle of the anomaly fields. 
Tire time-latitude diagram of the model zonal wind and precipitation are shown 
in figures 5.16 (90"E), 5.17 (J60"W) and 5.18 (80°W) for summer and winter SST 
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Time-Latitude structure at 90°E 
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Figure 5.19: Time-latitude structures for perpetual summer (first two panels from the 
top) and perpetual winter simulations (third and fourth panels). Negative contours 
are dashed. 


Precipitation 


simulations as indicated in the panels. A comparison of the annual cycles with 
summer and winter mean winds (5.14) at 90° E with the corresponding results 
for summer and winter SST simulations (5.16) shows that the northward shift 
seen over the monsoonal region in the summer months is essentially due to the 
annual cycle of SST. However, the fields do not return to southern latitudes when 
the annual cycle of mean wind is absent. The winter simulations are, however, 
generally much more periodic and equator centered, figure 5.18 shows that zonal 
wind and precipitation are predominantly centered in the Northern Hemisphere, 
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Time-Latitude structure at 160 °W 
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Figure 5.20: Same as the figure 5.18, but for 160°W 

although the SST maximum would be in the Southern Hemisphere. In contrast, 
over 90 ° Fj and 160°W the fields for winter SST simulations are more or less equa- 
tor centered. These features point to the conclusion that the annual cycle of the 
total anomaly fields in the tropics is a complex combination of effects due to the 
annual cycles of the mean wind and SST. 
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Time-Latitude structure at 80°W 



Figure 5.21: Same as the figure 5.18, but for S0°W 


5.3.3 Perpetual Summer and Perpetual Winter Simula- 
tions 

The time-latitude structures of the model simulation when annual cycle of both 
mean wind and mean SST are suppressed are given in figures 5.19 to 5.21, for the 
three longitudes. Each of the figures shows zonal wind and precipitation for the 
perpetual summer (summer mean wind and summer SST) and perpetual winter 
(winter mean wind and winter mean SST). 

As expected, in the absence of both the mean annual cycles, the anomaly cir- 
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dilation shows no significant change with season. However, the simulated time- 
latitude structures of the model fields are quite different over three longitudes, 
emphasizing the role of the spatial structure of the mean fields in the anomaly 
circulation. In particular, both the monsoon and the Eastern Pacific regions are 
characterized by equatorial westerlies for perpetual summer simulation, while the 
hastern Pacific region is characterized by anomaly easterlies for perpetual sum- 
mer. Another interesting feature is a reversal of the zonal wind from westerly 
for perpetual summer to easterly for perpetual winter at 160°W, with a similar 
but reversed effect at 80°W. In contrast, there is no appreciable change at 9 (HE 
where the zonal wind is predominantly westerly for both perpetual summer and 
perpetual winter. A scrutiny of the zonal wind for perpetual winter over 160°W 
(5.20) and for perpetual summer over 80° W (5.21) shows the signature of a sea- 
sonal scab' variability in the form of a reversal of the (dominant easterly) wind 
to westerly winds. The significant fact is that these westerly bursts coincide with • 
episodes of precipitation over these locations. Since for perpetual summer and 
perpetual winter there is no external contribution to the temporal variability of 
the anomaly fields, this seasonal cycle is essentially intrinsic to the convectively 
forced tropical circulation. 

5.4 Conclusions 

The annual cycle of the anomaly circulation in the tropics appears to be essen- 
tially determined by the annual cycles of the mean wind and SST. Indeed, as we 
ha, vc seen with our perpetual summer and perpetual winter simulations, there 
is no significant change in the circulation pattern with the march of the sea- 
son in the absence of the mean annual cycles. This is not surprising, since our 
model does not contain any other seasonally varying forcing. This implies that 
variabilities associated with the 30-50 day oscillations or QBW do not have any 
overall contribution to the seasonal variations. I his point is examined further in 
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connection with the structure of intraseasonal oscillation. 

An important point brought out by our simulation is the important role played 
by the annual cycle ol the mean SST. The GCM simulations indicated that rela- 
tively small SS I anomalies over the equator can give rise to large diabatic heating 
through the release of latent heating (Rowntree, 1979; Julian and Chervin 1978; 
Kesavamurthy, 1982). 1 he atmospheric response to this diabatic. heating can also 
induce an annual cycle in the anomaly circulation. While the mean wind and SST 
may influence each other to produce the observed mean anomaly cycle, the an- 
nual cycle of the anomaly circulation is determined by a combined effect of the 
two mean annual cycles. Neither of the mean annual cycle alone can give rise to 
a. characteristic annual cycle. However, the mean SST and the mean wind in the 
summer and the winter affect the anomaly circulations differently. An interesting 
point is the higher degree of variabilities in presence of the summer mean fields. 
This has an important implications for the predictability in the tropics. 

It can be seen from figures 5.19-5.21 that there is , however, a weak annual 
cycle intrinsic to the convection induced internal dynamics. This is particularly 
prominent over the Pacific regions, with a rather prominent seasonal-scale (200 
days) variabilities. 
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Chapter 6 


The Spectrum of Oscillations I: 
The Broad Spectrum 


The existence of a broad spectrum of oscillations spanning the intraseasonal to 
interannual time scale makes modeling of the tropical atmosphere a challenging 
task. These oscillations are an integral part of the "tropical dynamics, and they 
have important implications for the climate and predictability in the tropics. The 
simulation of these variabilities, therefore, has been a. major testing ground for 
both (J(!M and for models of intermediate complexity. Nevertheless, most of the 
efforts in simulating the tropical intraseasonal oscillations have focussed primarily 
on 30-50 day oscillations. As we have seen, there are at least two other significant 
oscillations at intraseasonal time scale in the tropical atmosphere. In addition, 
the tropical atmosphere exhibits several longer period oscillations, prominent ones 
being the quasi-biennial oscillation (QBO) and the El Nino-Southern Oscillation 
(ENSO) type interannual oscillations. 

The detailed evolution and structure of these oscillations with widely separated 
time scales are, no doubt, influenced by external and internal forcings of relevant 
time scab's. For example, it is quite likely and is suggested by several studies that 
the evolution and the structure of the interannual variabilities are significantly 
affected by the ocean-atmosphere coupling. Indeed in the currently held scenario, 
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Period (days) Period (days) Period (days) 

Figure (). I : Power spectra of unfiltered time series of the model u field from composite 
annual simulation for uniform EWF plotted with vP(v) as the ordinate and period 
as abscissa. The solid line for r =0.5 hr. and dashed line represent r =1.5 hr. The 
corresponding red noise spectra are marked with symbols. 

the ENSO variability is an ocean-atmosphere coupled instability. However, the 
question of (/enrsift of these variabilities can be addressed separately. In partic- 
ular, if is possible that the different oscillations are excited by the same basic 
mechanism, while their detailed evolution and structure are modulated by the 
forcings operative at spatio-temporal scale of particular oscillations. 

We investigate these questions with our model of convective tropical atmo- 
sphere. Previous analytical studies have already shown that convection-induced 
internal dynamics in presence of a CTL of a few hours can selectively excite 
infra, seasonal oscillations at 30-50 day (Goswami and Rao 1993, heiealtei OR), 
10-20 day and 3-5 day scales (Goswami and Mathew 1994, hereafter GM). How- 
ever, for our model simulation of the monthly anomaly fields, it was found that 
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Figure 6.2: Same as figure 6.1, but for meridional wind 


the most appropriate value of CTL is between 0.5 and 1.5 hours. This is not 
surprising, since the monthly anomaly fields contain contributions not only from 
the low frequency waves but also from the convection organized at smaller scales. 
Indeed, it was found that while the most appropriate value of CTL for selective 
excitation of Kelvin wave at 30-50 day time scale is about twelve hours (GR), 


the corresponding value for the higher frequency ISO is about six hours (GM). 
Besides, the spectrum and the nature of oscillations in presence of the non-linear 
heating processes and the observed space-time dependent mean states may be 
quite different. An important question is whether convection induced internal 
dynamics could support any longer period oscillation without ocean-atmosphere 
coupling (but may be with SST-induced variations in distribution of convective 


feedbacks). In particular, if our model is to provide an effective simulation of the 
tropical anomaly climate, a single value of CTL (chosen from a range) should be 
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Figure 6.3: Same as figure 6.1, but for precipitation 

capable of simulating the entire range of observed oscillations under appropriate 
dynamical conditions. Therefore the question of the existence of oscillations at 
various time scales and the relative roles of various processes in their genesis and 
structure is investigated. 

6.1 Design of the Experiments 

Given that the observed circulation in the tropics exhibits a wide spectrum of 
oscillations, an important question to be answered concerns the minimal and rel- 
evant physics that can support these oscillations. In particular, it is important to 
identify and ensure mechanisms that can support certain well-known phenomena, 
like the 30-50 day eastward propagating oscillation. While the detailed structure 
of each major oscillations in the model will be addressed and described subso- 
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Figure 6.4: Power spectra of unfiltered time series of the model u field from composite 
annual simulation for uniform EWF plotted with vP{v) as the ordinate and period 
as abscissa. The solid line for r =3.0 hr. and dashed line represent r =6.0 hr. The 
corresponding red noise spectra are marked with symbols. 


cpiently, our primary goal hero is l,o examine the existence of these oscillations 
and to identify the primary mechanism of their genesis. Most of the earlier stud- 
ies have simulated, as we have mentioned earlier, only a particular oscillation 
or at the most a part of the spectrum at a given time. But for a model to be 
successful in simulating the tropical climate it must be capable of simulating all 
the relevant observed periodicities. This requirement is getting stronger in view 
of the interconnection among the oscillation with widely separated time scales. 
Hence, a large number of experiments to identify the effects of various dynam- 
ical and mean conditions on the genesis and structure ol tiopical oscillations is 
carried out. The analyses were done using model simulations fiom both one-yca.i 
integrations and much longer 38-year simulations. A brief description of these 
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Figure 6.5: Same as figure 6.4, but for meridional wind 


two sets of experiments is given below. 


6.1.1 Composite Annual Simulation 

To examine the existence and the characteristics of the intraseasonal oscillations, 
we have carried out a number of (sixteen) one-year simulations with different 
initial conditions. The initial conditions are obtained by subtracting the clima- 
tological mean monthly fields from the monthly mean fields for a particular year. 
For the one year simulations we have used the 9-year (1979-1987) data set from 
Comprehensive Ocean Atmosphere Data Sets (COADS). The monthly mean SST 
fields were also taken from the same data set. The sixteen initial conditions for 
the sixteen one-year simulations were obtained from eight winter (January) and 
eight summer (corresponding June) monthly anomaly fields obtained fiom the 
observed data. The results of the sixteen simulations were then composited with 
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Figure 6.6: Same as figure 6.4, but for precipitation 

equal weights to arrive at a composite one- year simulation. This procedure is 
repea, ted for each of the four values of CTL considered in our study. 

6.1.2 Long-term Simulations 

To generate a model simulation suitable for investigating the broad spectrum of 
tropical oscillations, we have integrated the model lor thirty eight years with a 
given observed initial condition. The observed initial fields are obtained by sub- 
tracting the :l8-year mean for January from the monthly mean January field for a 
particular year (1950). The mean fields of zonal and meridional wind components 
are given as the 38-year mean monthly winds obtained by' taking the mean or ob- 
served monthly mean surface fields for zonal and meridional wind components 
from the 38-year (1950-1987) COADS. The mean SST field is also taken from the 
same data, set. 
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Figured. 7: Power spectra of unfiltered time series of the model u field from composite 
annual simulation for SST dependent EWF plotted with vP(v) as the ordinate and 
period as abscissa. The solid line for r =0.5 hr. and dashed line represent r =1.5 hr. 
The corresponding red noise spectra are marked with symbols. 


We have used two different data sets for obtaining the mean states and initial 
conditions for the annual and long-term simulations with the purpose of obtaining 
a more objective evaluation of the model simulations. 

Both the annual (composite) and the long-term simulations were repeated to 
examine the existence and structure of the oscillations under different mean and 
dynamical situations. The detailed design of flic simulations with emphasis on 
the different aspects is discussed below, 
a) Role of Convective Timelag 

As we have argued earlier, the values of CTL used in earlier analytical studies 
may have to be changed in presence of non-linear processes and complex mean 
fields. Indeed, we have already seen in chapter 4 that it was with a smaller value 
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Figure 6.8: Same as figure 6.7, but for meridional wind 

of CTL; from 0.5 hour to 1.5 hour, that the detailed structure of the anomaly 
fields are better simulated. However, the adopted value of CTL must also be able 
to support the observed spectrum of tropical variabilities. Therefore, four values 
of CTL, from 0.5 hour to 6.0 hours, are considered. For the annual composite 
simulation, this gives 16x4=64 one- year simulations, as the model is integrated 
for sixteen different initial conditions for each value of CTL. For the long-term 
simulation, this results in four 38-year model integrations, 
b) Uniform Strength of Evaporation-Wind Feedback 

A spatially and temporally uniform strength of EWF represents a mean condi- 
tion in which the underlying SST field is also spatially and temporally uniform. 
In the case of one-year integration, this would mean the absence of the annual 
cycle and the spatial variation of SST field. In the case of long-term simulation 
this would imply spatially uniform SST field devoid of any seasonal as well as in- 
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Figure 6.9: Same as figure 6.7, but for precipitation 


terannual variabilities. As we have already seen, a uniform strength of EWF can 
not simulate the spatial structure of the monthly anomaly fields well, signifying 
the importance of the structure of the SST field in the simulation of the tropical 
variabilities. We shall examine this issue in the context of intraseasonal and in- 
terannual oscillations. Thus sixteen-integration-composite one-year simulations 
were generated for each of the four values of the CTL with a uniform strength of 

EWF. 

c) Role of SST-induced Spatio-Temporal Variations of EWF 

One of the moist processes that has been shown to play a significant role in the 
dynamics of tropical variabilities is EWF. This is also true in the present study, in 
the context of the structure of the anomaly fields. We have seen that the inclusion 
of dependence of EWF on underlying SST, through dependence of A q on SST, 
lias profound effect on the structure of simulated variabilities. The corresponding 
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Figure' 6.10: Power spectra of unfiltered time series of the model u field from long 
term simulation for total SST dependent EWF plotted with vP(v ) as the ordinate 
and period as abscissa. The top row is for zonal wind, middle for meridional wind 
and the last row for precipitation. The solid lines for r =0.5 hr. and dashed lines 
represents r =1.5 hr. The corresponding red noise spectra are marked with symbols. 


effects on the genesis and structure of ISO are investigated here. Both tire annual 
and the interannual integrations were therefore repeated to investigate the role 
of SST-induced spatio-temporal structure of strength of EWF on the genesis of 
the oscillations. 

d) Role of the Annual Cycle of the Mean Wind 

The mean circulation in the tropics is characterized by a pronounced seasonal 
reversal of winds over the monsoon region. It is known that this annual cycle of 
the mean wind significantly affects the anomaly circulation in the. tropics. To ex- 
amine the effects of the seasonal variation of the mean wind on the intraseasonal 
oscillations, we have carried out two one-year simulations with winter (l)ecem- 
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Figure 0.1 1: Power spectra of unfiltered time series of the model u field from long 
term simulation for SST dependent EWF plotted with vP(v ) as the ordinate and 
period as abscissa. The solid line for r =0.5 hr. and dashed line represent r =1.5 hr. 
The corresponding red noise spectra are marked with symbols. 


her, January and February mean) and summer (June, July and August mean) 
background winds with a given initial condition, 
e) Role of Interannual Variability of SST 

In the current scenario, many of the low frequency (longer than a year) variabil- 


ities of the tropical atmosphere are viewed as resulting from ocean-atmosphere 
interaction. While there can be little doubt that the ocean-atmosphere inter- 
action plays a crucial role in the structure and evolution of these low-frequency 
variabilities, it is worthwhile to investigate whether the question of genesis of 
some of these oscillations can be addressed in the context of (convection induced) 
atmospheric dynamics alone. To address this question, we have generated long 
term simulations without the effect of the interannual variabilities of SST. 
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Figure 6.12: Same as figure 6.11, but for’ meridional wind 

6.2 Method of Analysis 

'L’o verify the presence of both intraseasonal and longer term oscillations, the 
time series from the various numerical experiments are subjected to a power 
spectrum analysis. The significance of peaks are tested against the correspond- 
ing red noise spectra. Both the 38 -year and one-year simulations are subjected 
to this analysis. While the long term simulations are primarily used to investi- 
gate longer period oscillations, we also look for certain intraseasonal oscillations. 
These results from the long-term simulation on intraseasonal oscillations are then 
compared with the results for the same oscillation from the composite one- year 
simulations. This provides a cross-validation of the two simulations with different 
mean fields, initial conditions and periods of integration. For all the simulations, 
the time series were constructed by creating zonal area aveiages over three le- 
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Figure 6.13: Same as figure 6.11, but for precipitation 

gions: monsoon sector (60-120° E), Central Pacific (170°E-130°W) and Eastern 
Pacific ( I 1 0-50" W). These zonal averaged time series were then recorded at three 
latitudinal locations: centered around 16°N, the equator and 16°S. 

6.3 Results 

6.3.1 Intraseasonal Oscillations 

In view of the significant role played by the ISO in tropical dynamics, the ability 
to simulate some of the characteristics of these oscillations can be considered as 
a. prerequisite lor any model for tropical dynamics. We, therefore, first analyze 
our model simulations in terms of existence of the intraseasonal oscillations. A 
more detailed analysis of their zonal and meridional structures will be presented 
subsequently. Figures from 6.1 to 6.3 show the power spectra from the composite 
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Figure 6.14: Power spectra of unfiltered time series of the model u field from long 
term simulation for climatological SST dependent EWF plotted with uP(v) as the 
ordinate and period as abscissa. The solid line for r =0.5 hr. and dashed line represent 
t =1.5 hr. The corresponding red noise spectra are marked with symbols. 

one-year model simulation with a spatially and temporally uniform strength of 
EWF at different locations as indicated in the panels. The solid and the dashed 
curves in each panel represent the power spectra for the two values of CTL, 
0.5 and 1.5 hours respectively. The red noise spectrum corresponding to each 
value of CTL is represented by the corresponding curve with solid circles. The 
corresponding figures for r =3.0 and 6.0 hours are shown in figures 6.4 to 6.6. As 
can be seen from these figures, neither value of CTL can give rise to significant 
power at relevant inf seasonal periods. Especially, the power above the red noise 
level is not seen at 30-50 day scale- the most prominent intraseasonal periodicity 
in the tropics. As we have seen in chapter 4, a uniform distribution or strength of 
EWF fails to simulate the observed structure for any values of CTL. Along with 
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Figure 6.15: Same as figure 6.14, but for meridional wind 


our present finding of the power spectrum, it now appears that a proper simulation 
of the tropical variabilities must involve a spatial and temporal distribution of 
the strength of EWF. 

The power spectra with a (SST induced) space-time dependent strength of 
EWF are shown in figures from 6.7 to 6.9. Once again, the solid line and the 
solid line with solid circles represent, respectively, the power spectrum for r = 0.5 
hour and the corresponding red noise spectrum. The corresponding spectra for 
T - 1.5 hour are represented by the dashed line and the dashed line with filled 
circles. As can be seen from these figures, the power spectrum shows significant 
power in a number of frequencies, corresponding to 3-7 days, 10-20 days and 30- 
GO days. However, the power is not always above red noise level in all locations, 
except for the 30-50 day wave, which exhibits a strong signal everywhere. The 
quasi -biweekly (or 10-20 day) oscillation also shows considerable power in several 
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Figure 6.16: Same as figure 6.14, but for precipitation 

locations. However, there appears to be only marginal power, although it is above 
the red noise level in many locations, at 3-7 day range. However, considering that 
our time series are obtained by averaging over rather large zonal extents, a process 
which might have weakened the signal, we shall consider the detailed structure 
of the oscillations at all these observed periods. 

6.3.2 Interannual Oscillations 

The tropical atmosphere is also characterized by a number of long period oscil- 
lations. Prominent among these are the quasi-biennial oscillations with a period 
between 18-36 months and the ENSO-type interannual oscillations with a period 
between 3-6 years. There is also evidence for decadal and interdecadal oscillations. 
We have, therefore, examined our 38-year model simulations for existence ol os- 
cillations at various observed frequencies. To validate our long-term simulation 
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Figure 6.17: Power spectra of unfiltered time series of the model u field from one 
year simulation with winter and summer wind plotted with vP(v) as the ordinate and 
period as abscissa. The solid line for summer wind+annual cycle in SST and dashed 
line for winter wind-f annual cycle in SST. The corresponding red noise spectra are 
marked with symbols. 


and to chock consistency with our one-year simulations, we have first examined 
the time series (written every two days) from our long-term simulation for the 
existence of intraseasonal oscillations. The results of this analysis are presented 
in figure 6.10, for two values of r = 0.5 hour (solid line) and 1.5 hour (dashed 
line). These results represent the power spectra of unfiltered time series model 
zonal wind (top panels), meridional wind (middle panels) and precipitation (bot- 
tom panels) at the equator. As indicated in the' figures, the t ('suits ate shown foi 
the three zonal locations. Only the part of the power spectrum with a period less 
than annual time scale has been emphasized to examine the existence of piorni- 
nent peaks at intraseasonal period above red noise level (lines with solid symbols). 
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Figure G. 18: Same as figure 6.17, but for' meridional wind 

A |> s'ominenf feature of the power spectra is the existence of a broad-band peak at 
infr seasonal frequency for the zonal wind and precipitation. However, the power 
in I. lie meridional wind is generally below the red noise level at intraseasonal time 
seal « 3 s. This is consistent with the results from the composite one-year simulation 
whi <■' li show weak power for the meridional wind over most locations. The broad 
ban els for the zonal wind and precipitation are, of course, indicative of substantial 
[„!, ( w , -annual variations of the characteristics of 30-50 day oscillations, a point we 
iio(,< ..5 for future investigation. However, the important conclusion lor the piesent 
is IT i at the two simulations with their different initial condition and mean fields 
gen *rate simulations with a high degree of consistency. However, the long-teun 
sin, m „ ilation also exhibits a strong annual cycle with very large power. To delin- 
(ak _ the longer term oscillations with relatively weaker power, corresponding to 
(|ll;i:2 ^i-biennial, interannual (3-7 years) and longer term (7-8 years) we have taken 
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Figure 6.19: Same as figure 6.17, but for precipitation 


a 10-month running average of our model time series. The results are shown in 
figures from 6.11 to 6.13. Most of these peaks are prominent only for one value 
of CTL, viz 0.5 hour. This is consistent with our finding that the value of r = 0.5 
hour also generates the best simulation of the monthly structures of the tropical 
variabilities. 

6.3.3 Role of the Mean Annual Cycle 

There are two ways in which the mean annual cycle influences the dynamics in our 
model. While the annual cycle of the climatological mean winds is incorporated 
through the advection, the annual cycle of SRT enters through its effect, on the 
strength of EWF. In the following section we shall describe the relative roles of 
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Figure (>.‘20: Power spectra of unfiltered time series of the model u field from one 
year simulation with winter and summer SST plotted with uP(u) as the ordinate and 
period as abscissa. The solid line for summer SST-|-annual cycle in the wind and 
dashed line for winter SST+ annual cycle in the wind. The corresponding red noise 
spectra are marked with symbols. 


the annual cycles of the mean wind and SST. 


6.3.4 Role of Annual Cycle of Mean Wind 


The results of the experiment to determine the role of the annual cycle of the 
mean wind are shown in figure 6.17 to 6.19. The nine panels in figure 6.17 show 
the power spectra for the zonal wind for one-year simulations with only winter 
(December, January, February) and summer (June, July, August) averaged mean 
winds at different locations as indicated in the panels. The solid line represents 
the power spectra, lor the summer mean wind run and the solid line' with symbol 
represents the corresponding red noise spectrum. The powei sped. in. and the 
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Figure 6.21: Same as figure 6.20, but for’ meridional wind 


corresponding red noise for the winter mean wind simulation are represented by 
the clashed line and the dashed line with symbol, respectively. 

As can be seen from figure 6.17, there are prominent peaks at intraseasonal 
frequencies for both the summer and the winter simulations. However, in almost 
all cases the peaks in the power spectra are located around distinct frequencies for 
the winter and the summer simulations. Figure 6.17 shows a remarkable and per- 
sistent winter-summer asymmetry in the genesis of the intraseasonal oscillations. 
In particular, while the peak at the power spectra occurs almost always around 
50 days for the winter simulation, it is located at around 30 days for the summer 
simulation. The differential response to the winter and the summer mean winds, 
therefore, appears to be a primary mechanism for the broad band nature of the 
30-50 day oscillations. A more detailed analysis of the effect of the annual cycle 
of the mean wind on the other characteristics of the 30-50 day oscillations will be 
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Figure 6.22: Same as figure 6.20, but for precipitation 

discussed in connection with the structure of 30-50 day oscillations. The effect of 
the annual cycle of the mean wind, however, is much more prominent on the os- 
cillations of the meridional wind component, shown in figure 6.18. In particular, 
there is very little power above red noise at 30-50 days at most locations. I his is 
especially true for the winter mean wind (dashed line), and over the eastern Pa- 
cific region. The effect of the annual cycle on the tropical oscillations can be seen 
from figure 6.19, which shows the power spectra of precipitation with summer 
mean wind (solid line) and winter mean wind (dashed line). The corresponding 
red noise spectra are represented by the lines with symbols. As can be seen from 
this figure, the 30-50 day signal in precipitation is excited by only the summer 
(winter) mean wind in the monsoon sector (eastern Pacific). This feature can be 
verified by analyzing data over the two sectors specifically for the winter and the 
summer months. Indeed, it is already known that the 30-50 day signal over the 
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monsoon region is prominent only in the summer season. 

6.3.5 Role of the Annual Cycle of SST 

Although it is often argued that the oceanic processes do not have significant 
influence over the oscillations at intraseasonal time scales, such a scenario is 
not supported by observations. For example, Krishnamurti et al. (1988) found 
significant modulation of the air-sea fluxes at 30-50 day time scales. Significant 
oscillations at 30-50 day time scales of the oceanic variables, in particular SST, 
have been reported by a number of authors. Since in our formalism SST plays a 
significant role through its modulation of the strength of EWF, it is necessary to 
examine the effect of the spatio-temporal distribution of SST on the structure of 
the intraseasonal oscillations. The results of the experiment to determine the role 
of the annual cycle of the SST are shown in figures 6.20 to 6.22. The nine, panels 
in figure 6.20 show the power spectra for the zonal wind for one-year simulations 
with only winter (December, January and February) and summer (June, July and 
August) averaged SST at different locations as indicated in the panels. The solid 
line represents the power spectra for the summer SST-run and the solid line with 
symbol represents the corresponding red noise spectrum. The power spectra and 
the corresponding red noise for the winter SST simulations are represented by the 
dashed line and the dashed line with symbol respectively. The winter-summer 
asymmetry in terms of the positions of the peak, apparent in winter-wind and 
summer-wind runs, is not present in this case. There is clear dominance of power 
for the winter-SST run over summer-SST in most of the locations. However, both 
u and v (figures 6.20 and 6.21) spectra continue to exhibit much weaker power 
over the off equatorial positions. 

Role of Interannual Variability of SST 

At the interannual time scale, the most important effect in our simulation is 
likely to be the interannual variability of SST. The analyses of our model re- 


105 



suits from long-period simulation show that the ability of the model to generate 
iiiterannual oscillations at observed interannual periodicities. However, these in- 
terannual oscillations are likely to be the response of the atmospheric dynamics 
to the interannual variations in the moist feedback induced by the interannual 
variability in the SST, which in turn may be due to ocean-atmosphere interaction. 
Perpetual Winter and Perpetual Summer Simulations 

As we have seen above, the annual cycle of mean wind and SST effect in 
the genesis of ISO in somewhat similar but still has its own characteristics. To 
investigate the effect of absence of both the annual cycles (mean wind and SST), 
we have generated two one-year simulations with only winter mean wind and 
winter SST (perpetual winter) and with only summer mean wind and summer 
SST (perpetual summer). The power spectra at different locations exhibited 
(not shown here) peaks at intraseasonal time scale at several locations. However, 
the periodicity of the perpetual summer was found to be generally around 20 
days, while that for perpetual summer was about sixty days in many cases. This 
is consistent with our findings with the summer and winter mean winds that 
summer (winter) mean conditions favor short(long) period oscillations. 

6.4 Conclusions 

A noteworthy feature of the simulations discussed here is the persistent appear- 
ance of an oscillation at 30-50 day time scale for diverse dynamical conditions. 
It appears that the genesis of 30-50 day oscillation is not crucially dependent 
on the details of the mean field. However, the observed 30-50 day oscillation is 
characterized by several other properties, in particular, its well known eastward 
phase propagation. While the excitation of a 30-50 day and other intraseasonal 
oscillations may not crucially depend on the structure of the mean fields, the 
detailed properties of the oscillations may. We will discuss this in more detail in 
subsequent chapters. 
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Anothei effect of the annual cycle of the mean wind is that while there are 
significant peaks at higher frequencies for summer mean wind, there is very little 
power above 30-50 day time scale for the winter mean wind. Thus, the existence 
of ISO at 3-7 day and 10-20 day time scales may crucially depend on the annual 
cycle of the mean wind. 

The model is capable of supporting interannual oscillations at observed interan- 
nua.1 variabilities also. This is essentially in response to the interannual variation 
in the SST-induced convective forcing; the interannual variabilities disappear 
when the interannual variability in SST is removed (figures 6.14 to 6.16). While 
the question of excitation of interannual variability in SST has to be addressed, 
which is likely to be determined by ocean-atmosphere coupling, the success of 
the model in simulating these interannual oscillations without any explicit ocean- 
atmosphere coupling implies that our atmospheric model can respond to the SST 
induced forcing so as to give rise to the observed variabilities. 
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Chapter 7 


The Spectrum of Oscillations: 
The 30-50 Day Waves 


One of the most well investigated intraseasonal oscillations (JSOl) in terms of ob- 
servation is the 30-50 day oscillations. The 30-50 day oscillation was first detected 
by Madden and Julian (1971, 1972) in the spectral and cross-spectral analyses 
of time series of tropical station data and is now also referred to as the Madden 
and Julian Oscillations (MJO). Subsequently, a large number of observational 
studies reported the presence of a spectrum of oscillations with periods ranging 
from 10 to 90 days. Until recently most of the studies on ISOl were based on 
dal, a. sets which were sparce and insufficient for a systematic analysis of phenom- 
ena, on intraseasonal time scales. However, these studies had helped to bring 
out the broad characteristics of the ISO 1. As revealed by these observations, 
the 30-50 day oscillations take the form of an eastward moving circulation cell 
in the equatorial zonal plane and are most prominent in the zonal velocity and 
surface pressure, with the amplitudes maximum near the equator. The frequency 
responses are rather broad, with convective signals having comparable power in 
the time periods ranging from 30-95 days and the zonal wave numbers lying be- 
tween 1 and 3. Comparatively, the dynamical fields possess much narrower peaks 
concent rated in zonal wave number 1 and time periods of 50-00 days. We shall 
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Power spectrum with period centred at 40 days: u-field with SST-dependent EWF 



Period (days) Period (days) Period (days) 

Figure 7.1: Power spectra of zonal wind filtered at 30-50 days plotted for four values 
of CTL with vP{u) as the ordinate and period as abscissa, where P is the power and 
v is the frequency. For each value of CTL the fields have been normalized to their 
respective maxumum. 

refer to these eastward propagating oscillations with time period ol several weeks 
as ISO 1 which is, observationally, the most prominent intraseasonal oscillation in 
the tropics. Its other feature is that the wind perturbation appears to propagate 
around the globe but strong convection appears to be confined within the Indian 
Ocean/western Pacific sector. Detailed analysis of several recent data sets has 
helped to bring out the rich structure and dynamics of ISOl. Using the cross 
covariance between outgoing long wave radiation (OLR): wind and temperature, 
Hendon and Salby (1994) constructed a composite life cycle of ISOl. This dis- 
plays characteristics of both forced and radiating responses, lhe lorced response 
is basically a coupled Ilossby- Kelvin wave (Gill, 1980) that migrates eastward 
with the convection anomaly at about 5 vis across the eastern hemisphere. 
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Power spectrum with period centred at 40 days: v-field with SST-dependent EWF 



Period (days) Period (days) Period (days) 

Figure 7.2: Same as Fig. 7.1, but for meridional wind. 

The radiating response appears as a Kelvin wave, which propagates eastward at 
about 10 ms~ l into the western hemisphere away from the convective anomaly. 
Energetics deduced from the covariances between convection, temperature, and 
vertical motion indicate a close relationship between the relative phase of these 
quantities and amplification of ISOl. During the growth-stage of anomalous con- 
vection, the heating, temperature and vertical velocity are positively correlated. 
This is consistent with the relationship needed for growth via production- and 
conversion of eddy available potential energy (EAPE). Nevertheless, the heat- 
ing perturbation is nearly in quadrature with the temperature anomaly, which 
implies little production of EAPE after its maturity. 

As a modulator of precipitation, the ISO have a special appeal to the meteo- 
rologists belonging to monsoon regions where agriculture is the main source of 
income. Also, the fact that midlatitude circulations are telocoimectod to ISOl is 
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Power spectrum with period centred at 40 days: precipitation with SST-dependent EWF 
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Figure 7.3: Same as Fig. 7.1, but for precipitation. 

corroborated by many observations. There are some evidences that the formation 
of incipient vortex related to tropical cyclone may be related to ISOl (Nakazawa, 
1986). Considering its potentiality in influencing global weather systems, it is 
not surprising that special attention has been devoted to modeling and simulat- 
ing these tropical variabilities. However, even alter many years ol ellorts, the 
regional and global manifestations of ISOl are still evasive to modeling studies. 
Considerable insight has been gained on the mechanism ol selective excitation at 
preferred time scales of the oscillation and on the interaction of convection and 
the large-scale dynamics associated with it. However, a clear picture of organized 
convection and its relationship to the circulation in general is still lacking. These 
modeling studies have shed considerable light on the mechanism of genesis and 
evolution of ISOl. 

The most distinctive aspect of the oscillation, namely its slow eastward move- 
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Power spectrum with period centred at 40 days: summer u-field with SST-dependent EWF 
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Figure 7.4: Same as Fig. 7.1, but for summer months. 


merit, and equatorial emplacement of its maxima, prompted earlier investigators to 
search for an explanation and analytical formulation of the phenomenon in terms 
of Kelvin wave dynamics. Parker (1973) interpreted the oscillation as a Kelvin 
wave in an equatorial /7-plane mode studied by Matsuno (1906), although the ver- 
tical scale of low-frequency Kelvin wave is very small in the absence of heating. 
Murakami et al. (1984) studied the three dimensional structure and phase prop- 
agation of the 40-50 day oscillation and concluded that its major energy source 
is the conversion of potential energy to kinetic energy, llayashi and Sumi (198G) 
have studied 30-40 day oscillation in a GCM for ocean covered earth. They have 
used a flat lower boundary with a specified zonally symmetric temperature. The 


precipitation and zonal wind anomaly have a large eastward propagating wave 
number 1 component. This study suggests that the zonal inhomogeneities in the 
climatological tropical convection and orography are not essential to the oxis- 
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Power spectrum with period centred at 40 days: summer v-field with SST-dependent EWF 
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Figure 7.5: Same as Fig. 7.1, but for meridional wind in summer months. 

fence of the ISC) 1 . Using a fully stratified tropical model which includes cumulus 
friction and Hadley cell basic state circulations, they showed that some eastward 
propagation of tropical circulation features can occur even when the anomaly 
heating remain fixed in one position. 

A fundamental problem in invoking the Kelvin wave dynamics to explain ISOi, 
the apparent disparity in the phase speed, was addressed by Swinbank et al 
(1986). They have come up with some evidence that this may be understood in 
terms of moist Kelvin wave dynamics, if effect of latent heat release on stability is 
properly taken into account. A simple model-calculation leads to the conclusion 
that precipitation associated with the Kelvin wave reduces static stability and so 
decreases the phase speed. These conclusions were further confirmed by two GCM 
experiments in which latent heat-release was increased, first by increasing SST 
and second by doubling the value ol latent heat constant. In both ( a.ses the wave 
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Power spectrum with period centred at 40 days: summer precipitation with SST-dependent EWF 
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Figure 7.6: Same as Fig. 7.1, but for precipitation in summer months. 

period was increased. The role of the moist feedbacks, and in particular EWF 
on the stability and dynamics of the Kelvin wave was investigated by Neelin et 
al. (1987) and Emanuel (1987). Though this mechanism destabilizes the Kelvin 
wave, the selection at observed period could not be attained. 

One of the most successful mechanisms of genesis of ISOl was provided by 
Goswami and Rao (1994) involving Kelvin wave dynamics in presence of an ex- 
plicit moisture dynamics governed by a convective time lag. In particular, the 
model could successfully predict excitation of the Kelvin wave at 30-50 day time 
scale with observed spectral characteristics. Thus, our present model with a 
convective time lag is well suited for investigation of ISOl. However, as men- 
tioned above, the recent analyses reveal a. structure o! ISOl that has strong 
spatio-temporal dependence. These observed characteristics provide powerful 
constraints for model simulations as well as for any mechanism pioposcd for 
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Power spectrum with period centred at 40 days: winter u-fieSd with SST-dependent EWF 
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Figure 7.7: Same as Fig. 7.1, but for winter months. 


the 30-50 day oscillations. We, therefore, first analyze our model simulations in 
terms of existence and structure of the 30-50 day oscillation. As we have seen 


earlier, only an SST-dependent EWF can generate a satisfactory simulation of 
the tropical monthly anomalies. The significant power in 30-50 day oscillation is 
present only for SST dependent EWF, although the high frequency part of the 
spectrum peaks above red noise level even for a uniform strength of EWF . 1 hus it 
seems the SST-dependent EWF is more relevant for the dynamics of ISOl in gen- 
eral. Besides being considered the SST-dependent EWF, a few special cases will 
a, Iso be discussed to isolate and identify the relative role of the mean wind. We 
first consider the composite one-year model simulations for different dynamical 


conditions. A certain specific issue will then be addiessed using both inlet annual 


and one-year simulations from specific experiments. 
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Power spectrum with period centred at 40 days: winter v-field with SST-dependent EWF 




Figure 7.8: Same as Fig. 7.1, but for meridional wind in winter months. 

7.1 The Characteristics of the Power Spectrum 


As we have seen, our model of the tropical atmosphere suppoits a piomincnt 
signal at 30-50 days at different locations and under different dynamical condi- 
tions. Although this global scale phenomenon is present throughout the year, its 
strength and behavior vary with location as well as with the season. We shall 
address this issue first in terms of the power spectra of the composite model sim- 
ulation filtered at 30-50 days at different locations. The unfiltered time series 
' have been generated by averaging the fields over a zonal extent of 60° centered 
at 90°E (Monsoon sector), 80°W (Central Pacific) and 160°W (Eastern Pacific). 
We shall consider three latitudinal positions, viz. the equator, ICPN and 16°S. 

(a) Analysis of the 12-month data 

Figures 7.1 to 7.3 show the power spectra of the model u, v and precipitation 
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Power spectrum with period centred at 40 days: winter precipitation with SST-deoendent EWF 
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Figure 7.9:. Same as Fig. 7.1, but for precipitation in winter months. 

fields from the composite one-year simulation filtered at 30-50 days respectively. 
The results represent an analysis of an entire year data of at different locations. 
The four curves in each panel represent the power spectra for four different values 
of CTL. As can be seen from these figures, the model in general supports ISO 
with 30-50 day time scale. However, it is evident that only r = 0.5 hr. gives rise 
to the most significant signal at 30-50 days time scale. A much quoted feature 
of this ISO is its presence in the dynamical field throughout the globe. Indeed, 
as can Ire seen from figure 7.1 there are almost equally strong signals for all the 
three zonal sectors, especially over the equator. Another significant characteristic 
of the power spectra is its larger amplitude over the equator compared to that 
at higher latitude, implying an equatorially trapped nature. I he corresponding 
results for the model v are shown in the figure 7.2. I he conclusions siinilai to 
those for the model u also hold for the v. Hut the powet of v is i datively less 
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Time series of area-averaged model fields filtered around 40 days: SST-dependent EWF 
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Figure 7.10: The time series of zonal and meridional wind components and precipita- 
tion filtered around 40 days. The solid, dashed and solid with symbol lines represent 
10"N, 0° and 10°S respectively. 


than the corresponding value for the u. In contrast to the u and t>, the power 
spectra for precipitation show pronounced zonal as well as latitudinal variations. 
In particular, the signal over the monsoon region is much stionger than over 
the other two regions. Similarly, in the meridional direction the strongest signal 


appears over the equator. 


(b) Analysis of Summer and Winter Fields 

The activity of ISO shows pronounced seasonality. This seasonality is under- 
standable, if its major energy source is the large-scale convection associated with 
the seasonally migrating ITCZ. It was found that the related cloud activity favors 


ie summer hemisphere and 1, lie general location ol ITCZ (Zangvin, 1975; Wang 


and Rui 1990; Weickmann et al 1985). 


To examine this we have separately 
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Figure 7.1 1: The time-longitude structure of model fields from a composite one-year 

simulation at 10°N for r =0.5 hr.; filtered at 30-50 day. Negative contours are dashed 

1 

studied the power spectra of time series of the filtered model fields for the winter 
(January to May) and the summer (June to October) months. Figures 7.-1 to 
7.6 show, respectively, the power spectra of the filtered model u, v and precipita- 
tion fields for the summer months at different zonal and latitudinal regions. The 
corresponding results for the winter months are shown in figures 7.7 to 7.9. 

A comparison of the figures 7.4 and 7.1 shows that the equatorial trapping seen 

l 

for the u for the complete year is not supported by the summer fields, which show 
comparable and sometimes even more power in the higher latitude than over the 
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Figure 7.12: Same as Fig 7.11, but over the equator. 


equator. This difference is much more pronounced over the monsoon region than 
over the Pacific sectors. This observation has been made by Nishi (1989) and 
Handon and Sal by (1994) who noticed that the amplitude of the signals do not 
drop away from the equator as rapidly as the simulated tropical circulation does. 
These anomalous signals show high degree of symmetry when the convection 
neighbors the equator. In addition, more the anomalous convection away from 
the equator, greater its asymmetry at the equator. 1 his la.rt could be pointed out 
as one of the reasons why the power at higher latitudes in the monsoon region is 
more than over the equator in summer. 
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<'igure 7.13: Same as Fig 7.11, but over 1Q°S 


Murakami el al (1986) nUltaod OLR date over a tropical belt between 45°N 
and 45*S during the eight years of 1978-77 and 197943. They found that over 

the Northern and Southern I lemispheres in the monsoon region, modes m the 
20-90 thvy period band were more pronounce, 1 during the summer than winter. 

Tire characteristics of 3040 day oscillations also show consirlerahle zonal varia- 
tions. The analysis of OUl data (from 1974 to 1986) carried out by Lau and Chan 
(1988) revealed the following facts. The peaks in spectra of Obit in 30-50 day 
range are mostly eonlinod to the tropical area west of dateline (80“ -120*li). Along 
the euuator the signal is most pronounced in the Indian Ocean region and drop 
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i ime-Longituae structure at 0 
x=1 .5hrs; Filtered at 30-50 days 




Figure 7.15: Same as Fig 7.14, but over the equator 

observed period was found to change from 50 days in boreal winter to about 35 
days in boreal summer over the Indian ocean regions (Hartmann el al, 1992). A 
comparison of figures 7.4 and 7.7 dearly shows that, over the equator and 16°S 
in the monsoon sector the spectral peaks of zonal wind shift from a higher period 
in winter to a lower one in mm.ll.er. In particular, over the equator this shift is 
from about 50 to 115 days, in accordance with Hartmann et al. 's finding. The 
most remarkable difference in the power spectra of the summer precipitation field 
from those for the full year is the absence of any significant power in the southern 
Hemisphere. I, ike in mean convection, the anomaly convection is dominant III the 
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Figure 7.16: Same as Fig 7.14, but over 10°S 

s 1 1 m rn e r h ern ispliere. 

(c) Time Series of Model Fields 

To investigate the phase relationships between various model fields, we have 
examined the time series of the filtered model fields at different locations as shown 
in figure' 7.10. lu eaeli panel, we consider variables at three latitudes: equator 
(dash line), 16°N (solid line) and 16°S (solid line with symbol). As before, the 
model fields have been averaged over three zonal regions as indicated in the 
panels. The top, middle and the bottom panels represent, respectively, the zonal 
wind, meridional wind and the precipitation field, as indicated in the figures. It 
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Figure 7.17: The time-latitude structure of model fields from a composite one- 
year simulation at the 90°E for r =0.5 hr.; filtered at 30-50 day. Negative contours 
are dashed 


is evident from the figure 7.10 that zonal wind at all the latitude in the Central 
Pacific and Eastern Pacific are almost in phase. But, no such clear relationship is 
seen among the zonal wind at three latitudes in the monsoon sector. In meridional 
components the striking feature is the out of phase nature of win’d at 16°N and 
16°S. This statement is more valid in the monsoon sector and Central Pacific. 
This implies a converging/diverging pattern with a periodicity of 30-50 days. 
The amplitudes of all three variables undergo modulation with period more than 
100 days. This amplitude modulation is very clear in the Central Pacific. 
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Time-Latitude structure at 160°W 
T=0.5 hrs; Filtered at 30-50 days 



Time (months) 

Figure 7.18: Same as the Fig 7.17, but over 160°W 

7.2 Time-Longitude Structure 

After presenting the power spectra and filtered time series, it would be more 
revealing to look at the time-longitude structure of the various dependent variable. 
An analysis using the Movrn filler diagram would make the case stiougei, if the 
simulated tropical disturbances with maximum power in the period mentioned in 
the last section agree; also with observed- direction ol pha.se' propagation. 1* igtues 
7.11 to 7.13 show the time-longitude structure's of the filtered model fields for 
r = 0.5 hr. at three latitudes: equator, 16°N and 1G°S. The figures 7.14 to 7.16 






Time-Latitude structure at 80°W 
1=0.5 hrs; Filtered at 30-50 days 



Time (months) 

Figure 7.19: Same as the Fig 7.17, but over 80°W 

are for r = 1.5 hr. All the variables are bandpass filtered at 30 to 50 days periods. 
From a scrutiny of the figures from 7.11 to 7.16 the following observation can be 
made. The fields at the equator show characteristic oscillation throughout the 
year, with nearly equal strength. In contrast, the fields at the northern and the 
southern latitudes show a distinct seasonal variation. In particular, the fields in 
the southern hemisphere in the (northern hemispheric) summer (winter) months 
are rather weak. This is consistent with the results shown in the power spectrum 
analysis. The modeled variabilities propagate both westward and eastward. This 
agrees well with what Wang and Rui (1990) had reported. They studied OLR 
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Time-Latitude structure at 90°E 

x=1.5 hrs; Filtered at 30-50 days 




Figure 7.20: Same as the Fig 7.17, but for r = 1.5 


data and summarized various behaviors of large-scale tropical cloud complexes. 
About 122 tropical intra.sea.son a, 1 convection anomalies were seen during a 10-yeai 
period (1975-85). Of those, the majority moved cast (77), 27 moved northward 


with no connection to an eastward moving waves, and 18 moved westward. 

Several studies (ban and Oh an 1983, 198(5; Weickman cl al, l9So ) had con 
tributed toward better understanding of 30-50 day wave. Nakazawa (1988) stud- 
ied the structure of the eastward-moving cloud masses with 3 hourly geostation- 
ary out data. Ilis main findings were the following. The eastward moving cloud 
systems were composed of several super cloud clusters (SCO) which moved m the 
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Tsl .5 hrs; Filtered at 30-50 days 





Figure 7.21: Same as the Fig 7.20, but over 160°W 


same direction. Each SCC, in turn was, composed of a number of smaller cloud 
cluster (CC) that propagate westward. The area of each SCC was of the order 
of 1.0 :i km. and that of CC was an order of magnitude less. The life time of CC 
was only 1-2 days. The genesis of CC took place to the eastern side of matured 
one. He noted that this behavior was confined to within 15° of the equator. Be- 
yond that the cloud movement on all spatial scales tends to be westward. The 
OLR anomalies follow upper level divergence and are strongest over the Indian 
and Pacific Oceans. They are very weak over the cooler water of the Eastern 
Pacific and Atlantic Oceans and weak over south America and Africa. The OLR 
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Time-Latitude structure at 80°W 
t= 1.5 hrs; Filtered at 30-50 days 





Time (months) 


Figure 7.22: Same as the Fig 7.20, but over 80°W 


anomaly is found to be strongest in the summer hemisphere. Nakazawa (1995) 
confirmed the hierarchic structure of the tropical convection by examining its be- 
havior during TCKIA-COARE lOP from November 1992 to February 1993. The 
eastward propagation speed of super cluster was about () deg/da,y. Apart from 
eastward moving super clusters, also seen over central Pacific along the equator 
was westward moving super clusters. The phase speed ol these cloud-complexes 
were about 3 dog roe/ day. Another interesting feature revealed by using twice- 
daily .JM A/C ANAL global analysis data set was that active convective regions 
correspond well with lower-level westerlies and eonvei genre (uppci level caster 








Winter mean wind 




Longitudef*) LongihideO 


Figure 7.23: Time-longitude plot for winter and summer mean wind simulations but 
with annual cycle of SST. Only zonal wind and precipitation over the equator are 
shown. Negative contours are dashed 


lies and divergence). The ISO signals were strong in the region between Indian 
Ocean and the dateline. However, if planetary scale ISO was dominant, the signal 
was found all over the globe. Using space-time spectral analyses, Guber (1974) 
and Zangvil (1975) showed the evidence of eastward-propagating cloudy areas in 
the equatorial region. The latter also demonstrated that this was more promi- 
nent over the summer hemisphere. The x-t structures of precipitation shown in 
the figures 7.11 to 7.13 and the observational structures mentioned above stud- 
ies match remarkably well. A paper by Weickmann and Khalsa (1990) presents 
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Time-Longitude structure at 0°: t= 0.5 hrs. 


Summer SST 



Longitudef) Longitudef) 

Figure 7.24: Time-longitude plot for perpetual winter and summer SST simulations, 
but with annual cycle in mean wind. Only zonal wind and precipitation over the 
equator are shown. Negative contours are dashed. 

an analysis of the behavior of the 30-60 day oscillation between September 1981 
and April 1982 that uses raw OLR data and dynamical quantities subjected to a 
5- day running mean. The time-longitude diagrams of OLR anomaly clearly show 
three, regions of convection; a. large area covering Indian Ocean- Western Pacific 
region and two small regions over Africa and South America. This geographical 
confinement of convection is also very clear in the model simulated precipitation. 
What is very clear from Weiekmann and K balsa’s (1990) observational analysis 
and model simulation is that 30-60 day periodicity in the Indian Ocean- Wester n 
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Figure 7.25: Time-longitude plot for perpetual summer and winter simulations. Only 
zonal wind and precipitation over the equator are shown. Negative contours are 
dashed 


Pacific region is not determined by the time it takes for a propagating convective 
disturbance to make it around the world. Convection dies out over the Eastern 
Pacific. The times of intense convection over the Indian Ocean- Western Pacific 
region are separated by relatively quiescent periods that are approximately 30-60 
days apart. We also found that the eastward movement of the pattern is not 
a smooth propagation with almost constant speed. Instead, shifting of patterns 
may occur during the evolution. In the Northern hemispheric summer monsoon 
region between 60° and 160°E a series of time-clustered, partially space over- 
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lapping disturbances in association with 40-50 day disturbances have also been 
repotted (Murakami cl al . , 1984). 

besides eastward propagating signal, strong standing oscillations are also ev- 
ident in Western Pacific, Central Africa and South America where active con- 
vection is persistent (Hsu cl. al, 1990; Knutson and Weiokmann 1987). I-Yoni 
figures 7.11 to 7.16 we see not only this regionality ol the stationary mode but 
also the propagating signals that radiate from them. Other noteworthy features 
are a) the concomitancy of precipitation and westerly wind, b) at 10°N and 10°S 
evaporation is out of phase with zonal wind whereas at the equator they are in 
phase. 

7.3 Time-Latitude Structure 

In addition to zonal propagation, the poleward and equator ward movement of 
the signals were also revealed in many observational analysis. Wand and Rui 
(1990) further classified the eastward moving complexes into three, based on 
their propagation characteristics. First, those confined within ±15° of the equator 
during their life time, which were named as ME mode. The NE and SE modes 
move eastward along the equator until 100° E, whereupon they traverse in north- 
east and north-south direction respectively. Third type (EN) were seen over 
the Indian and /or the Western Pacific, which exhibited complex eastward and 
northward movements. The. EE mode occurred during the (i months period from 
December to May. Members of the second group that moved to the southeast near 
100" IC (Shi) were' reported only from November through April, indicating relation 
to the ITCZ and Australian monsoon. Those that moved to north-east (NE) 
occurred primarily, but not exclusively, in northern summer, again associated 
with the 1TOZ and the. East Asian monsoon. The major formation region is 
west, -centra] equatorial Indian Ocean. There is also, a secondary source just west 
of equatorial Africa,. The maximum intensification rate of the 1 1C A were in the 
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central Indian Ocean, with secondary areas of intensification near 160°E north of 
the equator and the one extending from Australia to the date line south of the 
equator. 

To find how well the model can capture this aspect of propagation we also give 
in figures from 7. L7 to 7.22 the time- latitude structures of the filtered model fields 
at three longitudinal locations: 90°E, 160°W and 80°W. Out of this six figures, 
first three are for r =0.5 hr. and the rest for r =1.5 hr. All the variables are 
bandpass filtered at 30 to 50 day peidods. Figures 7.17 to 7.22 reveal that all 
the fields are equatorially trapped. The simulations show distinct seasonal cycle 
at all the three longitudes. The meridional migration is largest in the Monsoon 
sector and decreases from east to west. In the Eastern Pacific this always remains 
in the northern hemisphere, showing that the anomaly zonal wind, to a certain 
extent, follows the mean seasonal cycle. 

The maximum of the zonal wind is closer to the equator in winter than in 
summer, in particular in the Monsoon sector and Eastern Pacific. During January, 
in the Central Pacific the zonal wind shows northward propagation emanating 
from 5°S. In summer the zonal wind moves equatorward in the monsoon sector 
whereas in the Central Pacific and the Eastern Pacific it also shows northward 
propagation. A very distinct intrusions of waves from the subtropics are visible 
in the Eastern Pacific. These characteristics are, to a large extent, applicable to 
all other variables shown in the figures. 

7.4 The Role of Mean Fields 

An important question of the source of energy for the sustained cumulus heating 
(at the subseasoned time scale) is usually not addressed. Convection for such long 
period requires a steady supply of moisture from the ocean. It is found that on this 
time scale the latent heat flux from the ocean to atmosphere is of considerable 
amount. The study by Krishnamurthi et al. (1988) indicates that oscillation 
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of both the Wind and the SST contribute to a strong coupling of atmosphere 
and the ocean at this time scale. The amplitudes of SST variation at this time 
scale are of the order of 0.5°(! to !.()''(! and most pronounced SST oscillations 
occur over the equatorial Western Pacific Ocean and the Bay of Bengal. It is 
conceivable that a one degree enhancement of SST over these climatological warm 
oceans can affect the convection and their propagation. Thus the spatio-temporal 
distribution of the SST fields is likely to play a crucial role in the structure 
of 30-50 day oscillation. In our model the SST-distribution affects tiie model 
dynamics through its implicit effect on the strength of EWF. In addition, the 
spatio-temporal structure of the mean winds is expected to have considerable 
influence on the regional characteristics of the ISO I - 

To investigate this aspect in detail, we have carried out two simulations, first 
with only the annual cycle of the mean wind (ie., with no annual cycle of SST) and 
then with only annual cycle of SST (ie., with no annual cycle of mean wind). In 
the former case we carry out two simulations, one with winter (January to May) 
mean SST and the other with summer (June to October) mean SST. Similarly 
in the case with only the annual cycle of SST we analyze two simulations: one 
with winter mean wind and the other with summer mean wind. In addition, two 
more simulations are carried out in the setting of perpetual winter (winter mean 
wind and winter SST) and a perpetual summer (summer mean wind and summer 
SST) mean conditions. In all these cases only one value of CTL is used (r =0.5 
hr). As we have seen earlier, there is significant power above the red noise at 
30-50 day time scale, for all these special cases. Hence, we shall consider here the 
simulated fields filtered at 30-50 day scale lor further analysis. 

'The time-longitude diagram for zonal wind and precipitation at the cquatoi lot 
summer and winter mean winds are shown in 7.23. I he corresponding figures lor 
summer and winter SS r l’ and perpetual summer and perpetual winter are shown in 
7.24 and 7.25, respectively. It was found that the power spectrum had significant 



peaks around 50 days for winter mean wind and winter SST and around 30 days 
for summer mean wind and summer SST. The model fields for two fields have 
been accordingly filtered. The most striking feature in the figure 7.23 is the 
absence of any clear eastward propagation. Indeed, both for the summer and the 
winter mean winds, the fields now exhibit westward propagation, irrespective of 
the zonal locations. T lit absence ol zonal variations is particularly striking in the 
winter mean wind simulation. Thus the spatial structure of the mean wind alone 
is not sufficient to bring about the observed variations in the zonal characteristics 
of the 30-50 day waves. 

An interesting conclusion that emerges from an examination of figure 7.24 for 
summer and winter SST is a more or less stationary behavior of the fields. Thus, 
from a synthesis of the results from 7.23 and 7.24, we may state that the phase 
propagation itself is induced by the annual cycle of the SST (and in particular 
associated temporal structure of EWF) while the direction of phase propagation 
is determined by the annual cycle of the mean wind. 

7.5 Conclusions 

A large part of intraseasonal variability in the tropics occurs in the 30-50 day wave. 
It is also intimately connected with the variabilities of monsoon rainfall. As de- 
tailed analyses of observed data sets bring out richer spatio-temporal structure of 
the 30-50 day oscillation, the emphasis of the model simulations must include not 
only genesis but also these detailed properties of 30-50 day oscillation. It is note- 
worthy that our model using convection-induced dynamics with convective time 
lag succeeds in simulating many of the observed features of the 30-50 day oscilla- 
tion. This is consistent with and reflective of the fact that our model can simulate 
a large part of the variability of the monthly anomaly fields. However, there are 
also certain important aspects of 30-50 day oscillation which our simulations do 
not capture well. One such feature is the well known northward propagation 


137 



of 30-50 day oscillation over the monsoon region. As can be seen from figure 
7.17, the model fields filteied at 30-50 days show weak northward propagation 
over the monsoonal region in the summer months but equatorward propagation 
in the winter months. This feature is most pronounced in the meridional wind 
field, while it is nearly absent in the evaporation field. The fields at 90°E are 
also much more equatorially confined with no significant amplitude beyond 15°N. 
In conti ast, the fields at 80 W show much more of! equatorial distribution with 
pronounced equatorward propagation in the northern hemisphere. Verification 
of some of these features may need further analysis of observed data with higher 
spatial resolution. Another very important aspects of 30-50 day oscillation and 
tropical oscillation, in general, is the interconnection between oscillations at dif- 
ferent time scales. Our model with its ability to simulate the broad spectrum of 
tropical oscillations is suitable for investigating this aspect. Some of these issue's 
will be addressed in the subsequent chapters. 
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Chapter 8 


The Spectrum of Oscillations: 
The Quasi-Biweekly Waves 


The existence of certain prominent periodicities in the summer monsoon rainfall 
has been known for several decades. In particular, the Indian monsoon region 
exhibits a 10-20 day (or quasi-biweekly, QBW) oscillation. This oscillation is char- 
acteristic of the summer monsoon region and significantly affects the distribution 
of monsoonal rainfall (active and break conditions), and is not well documented 
over other parts of the tropics. Although the QBW is not reported in the winter 
hemisphere, it is possible that more detailed analyses of now available observed 
data sets for the global tropics may reveal new features. With our model it is 
possible to investigate the existence -and spatio-temporal structure of this oscil- 
lation in detail. The agreement of the simulated features with observations then 
provides a basis for further investigation of the predicted features in the observed 
data. 

The questions regarding the characteristics of the QBW range from one of the 
genesis of these waves to the variation of the characteristics of these waves in space 
and time. While there has been a large number of theoretical investigations of the 
30-50 day eastward propagating wave, very few studies have addressed the ques- 
tion of dynamical mechanism of 10-20 day wave. A consistent dynamical scenario 
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Power spectrum with period centred at 1 5 days: u-field with SST-dependent EWF 



Period (dava) Period (dava) Period (days) 


Figure 8.1: Power spectra of zonal wind filtered at 10-20 days plotted for four values 
of CTL with vP(u) as the ordinate and period as abscissa, where P is the power 
and v the frequency. For each value of CTL, the fields have been normalized to their 
respective maximum. 

for the genesis of 10-20 day wave as well a.s for another major intraseasonal wave 
( viz. the 3-6 day wave) was advanced by Goswami and Mathew (1994). In par- 
ticular, it was shown in Goswami and Mathew (1994) that the same mechanism 
that selectively destabilizes the MRG wave at 3-6 day time scale in presence of 
mean easterlies excites the 10-20 day wave in presence of mean westerlies. Since 
the mean westerlies are characteristic of the Indian summer monsoon region, it 
provides an explanation for the presence of the 10-20 day wave over the north- 
ern summer tropics. A significant finding in Goswami and Mathew (1994) was 
that the 10-20 day wave appeared as a mode distinct Irom either MRG wave or 
Kelvin wave. Indeed it was shown that this mode can be expected oidy in a moist 
tropical atmosphere in presence of mean westerlies. These analytical results thus 
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Power spectrum with period centred at 1 5 days: v-fieid with SST-dependent EWF 



Period (days) Period ( days) Period (days) 

Figure 8.2: Same as the fig 8.1, but for meridional wind 

provide a conceptual basis for examining the existence and the characteristics of 
the 10-20 day wave in our model with CTL and non-linear heating. In particular, 
we would like to address certain questions regarding the detailed structure of the 
10-20 day oscillations which cannot be addressed with simple analytical models. 
The availability of observational analyses from recent works make such an effort 
worthwhile. To analyze the general properties of the 10-20 day oscillation in our 
model simulation, we shall once again consider the composite one- year model sim- 
ulations for different dynamical conditions. In order to investigate the existence 
of any significant oscillation in the 10-20 day period, first the time series from 
the composite one- year model simulation has been analyzed. It was shown that 
with r =0.5 hour, there exists a significant peak in the 10-20 day range above 
the corresponding red noise level. We shall, therefore, investigate the charac- 
teristics of the 10-20 day oscillation in our model simulation with the available 
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Power spectrum with period centred at 1 5 days: precipitation with SST-dependent EWF 



Period (days) Period (days) Period (days) 

Figure 8.3: Same as the Fig. 8.1, but for precipitation 

observations acting as the constraints in the model. We have seen that, in the 
context of monthly anomaly fields that an SST-dependent evaporation-wind feed- 
back (EWF) gives rise to the best simulation. This point has been also verified in 
the case of simulation of the annual cycle. Hence, we shall consider only the case 
of SST-dependent EWF to investigate the structure of QBW in our simulation. 
Hut four values of CTL are used to examine the effects of (!Th and compare this 
to its effect on the 30-50 day wave. 

8.1 The Characteristics of the Power Spectrum 

While most earlier analyses focussed on the 10-20 day oscillation of the summer 
monsoon rainfall, recent analyses with more comprehensive data sets indicate that 
the QBW may be more generic for the tropical circulation. Preliminary analyses 
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Power spectrum with period centred at 1 5 days: winter u-field with SST-dependent EWF 



Figure 8.4: Same as the fig 8.1, but for winter months 

of our model simulation also indicate significant power at the QBW range above 
red noise at several locations. Hence the behavior of the 10-20 day wave in our 
model simulation is examined as a function of location and season. This issue is 
addressed in terms of the power spectra of the composite model simulation filtered 
at 10-20 days at different locations. In addition, we have also analyzed the model 
simulations separately for the winter (January to May) and the summer (June 
to October) seasons filtered at 10-20 days. As before, the unfiltered time series 
were generated by averaging the fields over a zonal extent of 60° centered at 90° E 
(Monsoon sector), 80°W (Central Pacific) and 160°W (Eastern Pacific). We shall 
consider three latitudinal positions, viz. equator, 16°N (NH) and 16°S (SH). 

Figures 8.1 to 8.3 show the power spectra of the model u, v and the precipitation 
fields from the composite one-year simulation filtered at 10-20 days. As seen 
from these figures, the model in general supports ISO with 10-20 day time scale, 
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Power spectrum with period centred at 1 5 days: winter v-field with SST-dependent EWF 



Period (dava) Period (davs) Period fdavs) 


Figure- 8.5: Same as the Fig 8.4, but for meridional wind 

although the power is relatively weaker than that of the 30-50 day waves in the 
model simulation. However, since our model fields are normalized to its maximum 
value for ease of comparing results for different values of CTL, the smallness of 
power is not necessarily an indicator of weakness of the spectrum in the absolute 
sense. From figures 8. 1-8.6 we see, once again, that it is r = 0.5 hr. that gives 
rise to the strongest signal at 10-20 days time scale. 

A striking feature in figure 8.1 and, to some extent also in figures 8.2 and 8.3, is 
the comparable power at 10-20 day period at all the longitudinal positions. This 
may a.t first glance appear to be. eotil.ra.ry to observational analyses which show 
the QBW more prominently over the mousooual region. However, it is possible 
that the variability in the QBW scales at particular location also depends on 
the seasons. This question can be examined by analyzing the model fields for 
winter and summer separately. The figures 8.4 to 8.6 present the results of power 











Power spectrum with period centred at 1 5 days: winter precipitation with SST-dependent EWF 



Period (davs) Period (davs) Period (davs) 

Figure 8.6: Same as the Fig 8.4, but for precipitation 

spectrum of variables u, v and precipitation for a 10-20 day for the winter months 
respectively. The corresponding results for the summer months are shown in 
figures 8.7 to 8.9. As indicated in the figures, the results are presented for three 
latitudinal regions, around 16°N (top panels), around the equator (middle panels) 
and around 16°S (bottom panels), and also for three zonal locations: monsoonal 
area (left row), Central Pacific (middle row) and the Eastern Pacific (right row). 
In both the figures, there is ample power around 10-20 days in the monsoonal 
region, but very small for the other zonal regions. Indeed, the power in general 
monotonically decreases as we move from the monsoon region towards the eastern 
Pacific region.' In addition, comparison of figures 8.4 to 8.6 with the corresponding 
figures for the summer months reveals that while in the summer months the power 
is higher in the Northern Hemisphere, it is higher in the Southern hemisphere in 
the winter months, as expected. This may imply that the 10-20 day oscillation 
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Power spectrum with period centred at 15 days; summer u-fieid with SST-dependent EWF 



Figure 8.7: Same as the Fig 8.1, but for summer months 


appears as a response of the anomaly circulation to the changing mean conditions 
over the summer monsoon region. In all the figures from 8.1 to 8.6, the results 
are shown for four values of CTL adopted in this study. It is noteworthy that 
only one value of CTL, viz. 0.5 hour, gives consistently good results. 

The figure 8.7 shows the power spectrum of the model u filtered at 10-20 days 
for summer months. As we see from this figure, the summer u shows strongest 
signal only over the monsoon region in general. A similar conclusion holds for 
the summer v, shown in figure 8.8, although the amplitudes over the Southern 
Hemisphere are somewhat weaker. In contrast to the u and the v fields, the power 
spectra for the summer precipitation, shown in figure 8.9, show marked variations 
with latitude. In particular, only two locations, viz. the Nil monsoon region and 
the equatorial Central Pacific region, support a significant signal in the summer 
precipitation field. We have not yet come across any detailed observational analy- 
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Power spectrum with period centred at 1 5 days: summer v-field with SST-dependent EWF 



Period (days) Period (days) Period (days) 

Figure 8.8: Same as the Fig 8.7, but for meridional wind 

sis for precipitation with which we could compare our finding. However, with the 
availability of global high resolution data sets these predictions can be verified. 

8.2 Time Series of Model Fields: 

To investigate the phase relationships between various model fields, we have ex- 
amined the time series of the filtered model fields at different locations as shown 
in figures 8.10 and 8.11. In each panel, we consider variables at three latitudes: 
16°N (solid line), equator (dashed line) and 16°S (solid line with symbol). As 
before, the model fields have been averaged over three zonal regions as indicated 
in the panels. The top, middle and the bottom panels represent, respectively, 
the zonal wind, meridional wind and the precipitation field. The filtered time 
series is shown for both winter and summer separately (figures 8.10 and 8.11). 
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Power spectrum with period centred at 1 5 days: summer precipitation with SST-dependent EWF 



Period (days) Period (days) Period (days) 

Figure 8.9: Same as the Fig 8.7, but for precipitation 


These figures reveal several interesting features. The zonal components of wind 
at all latitudes considered are in phase in all the three longitudes shown, albeit 
the amplitudes decrease from east to west. In winter months amplitudes are of 
comparable magnitude in all three longitudes. This is not true in the case of the 
summer months (figure 8.11). While in Eastern Pacific zonal wind at 16°N has a 
large value in the case of Monsoon sector it is the equatorial zonal wind that dom- 
inates. There seems be to less change with season in the Central Pacific. In terms 
of phase relationships among the components of zonal wind at different latitudes, 
the winter season shows a more consistent picture; all arc more or loss in phase'. 


The meridional components at 16°N and 16°S are out of phase irrespective of the 
season. Nevertheless, a, dramatic change takes place in precipitation Irom sum- 
mer to winter: In summer precipitation is weak in Southern Hemisphere wheicas 
it is so in Northern Hemisphere in winter. This shows the seasonal migration of 
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Time series of area-averaged model fields filtered around 1 5 days: SST-dependent EWF 



Time (months) Time (months) Time (months) 


Figure 8.10: The time series of zonal and meridional wind components and precip- 
itation filtered around 15 days for the winter months. The solid, dashed and solid 
with symbol lines represent fields at 10°N, 0° and 10°S respectively. 

convective anomaly with season, similar to the mean precipitation. Precipitation 
over equator in Monsoon sector and Central Pacific do not change much from 
winter to summer. 

8.3 Time-Longitude Structure 

The westward propagation of the 10-20 day monsoon mode first was reported by 
Krishnamurti and Ardunay (1980). There is now considerable information about 
the zonal phase propagation characteristics of the QBW. Chen and Chen (1993), 
for example, report clear westward propagation of 850 mb meridional wind filtered 
at 10-20 days at 5°N. However, at 20° N, and especially over the Monsoon sector 
in the summer months, eastward propagation is also seen between 60-120°E. In 
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Time series of area-averaged model fields filtered around 15 days: SST-dependent EWF 



Time (months) Time (months) Time (months) 


Figure 8.11: Same as the Fig 8.10, but for summer months 

fact, this eastward propagation in the analysis of Chen and Chen (1993) is much 
more pronounced between 80 and 120°E and during May- July. Although the 
analysis of Chen and Chen (1993) is only for a single summer (of 1979), and may 
not be typical, it provides a good testing ground for our model simulation of the 
characteristics of the 10-20 day mode. 

The time-longitude structure of the filtered model fields for the entire twelve 
months of the model simulation are shown in figures 8.12 to 8.M, (or three lat- 
itudinal locations. These figures indicate a mixture of eastward and westward 
phase propagation for the 10-20 day oscillation. In particular, the fields show an 
essentially westward propagation over most longitudes over the equator. 1 his Ma- 
ture is, thus, in close agreement with the finding of Chen and Chen (1993) whose 
analysis shows clear westward propagation over 5°N, i.e close to the equatoi. In 
contrast, the fields at 16°N show, especially during the summer months, and ovei 
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ime-Longrtude structure at 1 0°N 
i=0.5hrs; Filtered at 10-20 days 
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Figure 8.12: Time-longitude structure of the model fields from a composite one-year 
simulation at 10° N for r =0.5 hr. Negative contours are dashed. 

80- 160° E, shows a clear eastward propagation. This eastward propagation begins 
to reverse to westward propagation around the end of August. Our model simu- 
lations are, thus, in excellent agreement with the observed features reported by 
Chen and Chen (1993). Indeed, the conspicuous stationary belt centered around 
40°E at the equator and more prominently at 16 D N, can also be seen in the anal- 
ysis of Chen and Chen (1993). The time-longitude structure at 16°S (figure 8.14) 
shows that the fields are rather weak over the Southern Hemisphere, especially 
in the summer months, thus making the 10-20 day mode more characteristic of 
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Time-Longitude structure at 0' 

t=0.5hrs; Filtered at 10-20 days 
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Figure 8.13: Same as the Fig 8.12, but over the equator 


the Northern Hemispheric summer. As can be seen from 8.13, the model u a.t, the 
equator exhibits strong zonal variation in its strength, while the corresponding 
meridional wind is more or loss uniform in the zonal direction. The time-longitude 
structure of the evaporation field at the equator is very similar to that of the u, 
although it shows a, more temporally uniform distribution titan the zonal wind 
field. This implies that the evaporation field at the synoptic scale is largely de- 
termined dynamically. In contrast, the precipitation shows much higher degree 
of spatial and temporal uniformity in its distribution. A significant feature of the 
fields at the equator is their relatively weak strength, especially in the zonal wind 
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Time-Longitude structure at 1 0°S 

T=0.5hrs; Filtered at 10-20 days 
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Figure 8.14: Same as the Fig 8.12, but over KPS 

component, during the summer, months at most longitudes. The fields at 10°N. 
shown in figure 8.12, however, exhibit characteristics quite distinct from those 
at the equator. For one thing, the fields show a much stronger, and somewhat 
reversed, seasonality than those at the equator. While at the equator the u has 
stronger amplitudes during the winter months at most longitudes, at 10° N, both 
the zonal and the meridional winds are strongest during the summer months. 
This contrast is strongest in the distribution of the meridional wind, which shows 
a rather uniform distribution in time and zonal direction at the equator. Perhaps 
the strongest distinction between the fields at the equator and 10°N is in terms 
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Time-Longitude structure at 1 0°N 
x=1 .5hrs; Filtered at 1 0-20 days 
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Figure 8.15: Same as the Fig 8.12, but for r =1.5 hr. 


of their phase propagation characteristics. While the fields at the equator quite 
generally exhibit a westward phase propagation, the fields at 10°N exhibit an 
eastward phase propagation in a characteristic manner. In particular, the fields 
tend to show an eastward propagation whenever amplitudes are large. Since this 
feature is present at several longitudes (and, in particular, in non-monsoon region 
where a seasonal reversal of the mean winds takes place), it is unlikely to be due 
to the mean winds. Thus, this feature is likely to be due to the dynamics of the 
system. 

An interesting feature of the structure of the model fields is revealed by figures 
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Time-Longitude structure at 0° 

tsI ,5hrs; Filtered at 10-20 days 
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Figure S.16: Same as the Fig 8.15, but over the equator 


8.12 and 8.13. It can be seen that the u is weaker at -the equator than at the 
Northern Hemisphere over the eastern pacific region, which would be consistent 
with a MRG wave structure for the u over this region. However, this feature is 
not present at other locations, especially over the monsoon region. The other 
feature is the presence of certain intense episodes, most prominent in the zonal 
wind and evaporation fields, both over the equator and the Northern Hemisphere. 

The corresponding results for a value of r=1.5 hr. are shown in figures 8.15 to 
8.17. Comparison of the results for r=0.5 hr. with the corresponding results for 
r=1.5 hr. reveals no significant change in the phase propagation characteristics 
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Time-Longitude structure at 10°S 
x=1.5hrs; Filtered at 10-20 days 



Figure 8.17: Same as the Fig 8.15, but over 10°S 


for the two values of r. 


8.4 Time- Latitude Structure 

We have seen that the model fields exhibit a pronounced annual cycle at certain 
longitudes, in agreement with observations. More we examine the characteristics 
of the annual cycle specific to the 1 0-20 day wave. This will allow us to assess the 
contribution of the 10-20 day wave to the annual cycle of the anomaly circulation 
in the tropics. In addition, this will also provide ns an insight into the detailed 
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t= 0.5 hrs; Filtered at 10-20 days 





Figure 8.18: Time-latitude structure of the model fields from a composite one-year 
simulation at 90° E for r =0.5 hr. Negative contours are dashed. 

structure of the- 10-20 day wave. The time-latitude plots of the model fields ' av- 
eraged around 90°E and filtered around 10-20 days are shown in figure 8.18. The 
corresponding results at 160° W and 80° W are shown in figure 8.19 and £>.20, re- 
spectively. A few striking zonal characteristics of the 10-20 day wave are discussed 
below. For example, at 90°E, the u is largely equator centered, although it also 
has prominent amplitudes in the Northern Hemisphere around 20 N, especially 
during the summer months. Thus this feature does not agree with the observed 
MRG wave structure of the 10-20 day wave, which should have u amplitude zero 
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Time-Latitude structure at 1 60°W 


t=0.5 hrs; Filtered at 10-20 days 
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Figure 8.19: Same as the Fig 8.18, but over 160°W 


over the equator. On the other hand, this structure agrees rather well with the 
classical easterly wave, which shows significant amplitude around 20°N. Further, 
the amplitudes at 20 ° N are strongest during the monsoon season. We have' noted, 
in connection with our discussion with the time-longitude structures, that this is 
also the time and location when the 10-20 day wave exhibit a. prominent eastward 
phase propagation. The picture emerges, th(' re lore, is that our model simulation 
contains both classical easterly waves and the westward propagating 10-20 day 
wave. Over the Central and the Eastern Pacific, on the other hand, the model u 
is prominently off-equator, while the v is equator centered. Ihus ovei this region, 
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Time-Latitude structure at 80°W 
t=0.5 hrs; Filtered at 1 0-20 days 
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Figure 8.20: Same as the Fig 8.18, but over 80°W 

the 10-20 day wave is essentially of MRG wave structure. Similar conclusions 
also hold for the Eastern Pacific sector, shown in figure 8.19. 

The time-latitude diagrams show a prominent double-cell structure in all the 
model fields at 90°E (8.18 and 8.21). The two cells, one at the equator and the 
other at about 20° N, are related through an almost 90° phase difference in all the 
model variables. Such double-cell structure, however, is very weak or almost non- 
existent over the Pacific sectors. Over 80° W, this double-cell structure still exists 
to some extent, but the out-of- phase relation between the two cells prominent over 
90°E has been now replaced by an in-phase relationship (8.20 and 8.23). Indeed, 
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Time-Latitude structure at90°E 
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Figure 8.21: Same as the Fig 8.18, but for r = 1.5 hr. 


over the 160°W, it is then, reasonable to conclude that the two cells merge at the 
mid-point of their separation to give rise to a single off-equatorial cell (8.19 and 
8 . 22 ). 

The time-latitude structures of the model fields also reveal that the annual 
cycle of the 10-20 day oscillations is quite different from the annual cycle of the 
unfiltered anomaly fields in general. In particular, over 90°F, we can discern two 
distinct components. One mode, with its a, centered at the equator, shows very 
little or no change of its latitudinal distribution with the march of the season. I he 
other, off-equatorial component, has a more distinct annual cycle. 1 his (summer) 
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Time-Latitude structure at 1 60°W 
t=1 .5 hrs; Filtered at 1 0-20 days 
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Figure 8.22: Same as the Fig 8.21, but over 160°W 


component has significant activity at around 2CFN, and has appreciable ampli- 
tude only during the summer months. Indeed, this Northern Hemisphere summer 
component of 10-20 day oscillation can be said to merge with the equatorial com- 
ponent with the march of the season as can be seen from figure 8.18. It is this 
second component of the 10-20 day oscillation that gives rise to the double-cell 
structure in the latitudinal direction. Such an equatorial expansion (or double- 
cell structure) in the meridional wind between /5-95°E was reported by Chen 
and Chen (1993) from an analysis of data generated by Global Data Assimilation 
System (GDAS) of the European Center for Medium- Range Weather Forecasts 
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Time-Latitude structure at 80 °W 

t= 1 .5 hrs; Filtered at 1 0-20 days 
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Figure 8.23: Same as the Fig 8.21, but over 80°W 

(ECMWF). Although the analysis by Chen and Chen (1993) was carried out for 
only one year, and may not reveal the interannual variabilities ol the characteris- 
tics of the 10-20 day inode, it is interesting that several ol the observed features 
are also present in our composite model simulations. It can also be seen that the 
genesis of the two components with distinct zonal characteristics makes a signif- 
icant contribution to the annual cycle ol the total anomaly fields. In paiticulai, 
the two components together determine the precipitation variabilities ovei the 
monsoon region at 10-20 day time scale. 
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8.5 Conclusions 


A remarkable feature of the model simulation is its ability to simulate distinct 
features for the winter and the summer months, as well as for different zonal 
and meridional locations. In present model the distinctions between the summer 
and the winter as well as between different zonal and meridional locations come 
essentially through the mean fields. Hence, the results imply that detailed char- 
acteristics of the 10-20 day wave are to a large extent a response to the mean 
conditions. The mean conditions affect our model dynamics both dynamically 
and thermodynamically. Thus both the dynamical and the thermodynamical ef- 
fects of the mean fields are essential for simulating the detailed structure of the 
intraseasonal variabilities. 

To determine the role of the mean fields and the annual cycles of the mean wind 
and SST, the model simulations with only winter and only summer mean states 
are examined. In contrast to the 30-50 day oscillation which shows a prominent 
signal even in the absence of a mean annual cycle, the QBW disappears in the 
absence of a mean annual cycle. Thus not only the spatio-temporal structure but 
even the genesis of QBW is strongly linked to the mean annual cycle. 

It appears that the contribution of QBO to the anomaly annual cycle is not 
negligible. Thus an inadequate simulation of QBW can adversely affect the sim- 
ulation of the anomaly fields which, in turn, can change the characteristics of the 
total circulation in the tropics. The QBW. like the ISOl is, therefore, likely to 
play a vital role in the tropical climate. 
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Chapter 9 

The Spectrum of Oscillations: 
Inter annual Variabilities 

The most prominent oscillation of the tropical climate at longer than annual scale 
is the ION SO related variabilities. The atmospheric part of the IONSO scale vari- 
abilities, the Southern Oscillation (SO) arises from an interannual variations in 
the exchange of mass between the Indian and the Pacific Ocean and has been 
known to the meteorologists for a long time (Walker, 1932). The SO is known 
to have strong correlation with the oceanic variabilities known as the El Nino. 
Owing to its large climatic and associated socio-economic effects, the ENSO phe- 
nomena have received wide and sustained attention from both the communities 
of meteorologists and oceanographers. However, even from a purely dynamical 
point of view the ENSO phenomena, with its rich structure, play a crucial role 
in the tropical climate system. Although, following the conceptualization by 
Bjerkness (1969) and its subsequent development, ENSO is often regarded as an 
ocean-atmosphere coupled variability; it is possible that the tjcncfns ol the SO is a 
response of the tropical atmosphere to (SST-induced) convective forcing. There- 
fore, next the investigation is on the existence and structure of the interannual 
variability in our model simulation. 

The SO variability is quasi-periodic, with a time scale between 3 and 7 years. It 


164 



is in this background of SO oscillations that the ENSO ( warm and cold events) 
that biings about such drastic regional and global scale changes in the atmo- 
spheric and oceanic circulations is embedded (van Loon and Shea. 1985: Lau and 
Sheu, 1988). The SO signal is characteristic of the global tropics, although the 
oceanic counterparts of ENSO variabilities. El Nino and the La Nina that give 
rise to extreme warm and the cold events are more pronounced over the Pacific. 
However, contrary to earlier belief. ENSO phenomena are also more global in 
character. Thus warm events, in phase with ENSO events over the Pacific but 
with weaker amplitudes also appear over the Indian Ocean (Yasunari. 1987). We 
shall leave the question of the role of ocean-atmosphere coupling in the structure 
of ENSO variabilities to a later chapter. 

While the mechanism and structure of ENSO related variabilities have received 
wide attention in recent years, the tropical atmosphere also supports other long- 
term variabilities. For example, the existence of a biennial variability in the 
atmosphere has been now known to the scientists for at least three decades. A 
summary of a number of earlier (prior to 1960) works on the existence of a bien- 
nial oscillation in the stratospheric and the tropospheric variables was prepared 
by Landsberg (1962). As it often happens in the atmospheric oscillations, the 
biennial variability was found to be quasi-periodic, with a time scale between 2 
and 3 years. A more comprehensive analysis of the existence of a QBO in the 
tropospheric variables was provided by Landsberg et al. (1963) from an analysis 
of time series of surface temperature at a number of widely spaced locations over 
the globe. The subsequent years have seen a surge of research works emphasizing 
various characteristics of the QBO and its relation to various other oscillations. 

A characteristic of the QBO is its rather global presence in a diverse set of 
variables. Thus QBO signal has been observed in East African rainfall (Rodhe 
and Virji, 1976) and United States surface temperature (Rasmusson et al. 1981). 
Biennial variability in the surface pressure has been documented by Trenberth in 
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Figure 9.1: The time-longitude structure of model fields from a long term simulation 
at KPN for r =0.5 hrs.; filtered around 18-30 months. Negative contours are dashed 

a series of studies (Trenberth 1975, 1976a, 1976b, 1980). Other studies (Trenberth 
and Paolino, 1981 and Trenberth and Shin, 1981) also provide ample evidence 
for the existence of a biennial signal in the tropospheric variables. 

The QBO thus differs from many other tropical oscillations in that it is not 
necessarily confined to the tropics, (liven this lad, it is not a priori obvious that, 
our model, with essentially convective dynamics, is appropriate for investigating 
the genesis and structure of the QBO. Nor it is clear that a single mechanism 
generates the QBO the world over. We, therefore, wish to explore the possibility 
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Time-Longitude structure over the equator 
T=0.5hrs; Filtered at 1 8-30 months 
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Figure 9.2: Same as figure 9.1, but over the equator 

of a tropical QBO sustained by the convection induced dynamics in response to 
the mean fields. 

The focus of the present chapter is on the structure of the interannual variabil- 
ity present in our model simulation and their comparison with observed structure. 
Also addressed is the question of essential mechanism necessary for the existence 
of interannual variability in our model simulation. To analyze the general prop- 
erties of interannual variabilities in the simulation, the long-term simulations for 
different dynamical conditions are considered. To investigate the existence of any 
significant oscillation in the interannual scale, first the time series from the inter- 
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Time-Longitude structure at 10°S 
t=0.5hrs; Filtered at 18-30 months 
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Figure 9.3: Same as figure 9.1, but over 10°S 

annual model simulation is analyzed. It was shown that both for r =0.5 hour and 
r = 1.5 hour there exist significant peaks in the. 18-30 months and 3-7 years range 
above the corresponding red noise level. However, as we have seen in our earlier 
studies, the best simulations of the characteristics of the tropical variabilities at 
monthly and intraseasona.1 timescales are provided by a. value ol CTL=0.5 hour. 
We shall, therefore, investigate the characteristics of the interaimual oscillations 
also for t = 0.5 hour. 
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Time-Longitude structure at 0°: Monsoon sector 
xsG.Shrs; Filtered at 48-72 months 



60 80 100 120 60 80 100 120 



Figure 9.4: The time-longitude structure of model fields from a long term simulation 
at 10° N for r =0.5 hrs.; filtered around 48-72 months. Negative contours are dashed 

9.1 Longitude-Time Structure 

Having seen the existence of significant signals at the two interannual periods in 
our model simulation, we have examined the phase propagation characteristics 
of these waves in detail to compare with observed structure. Hence, the time- 
longitude structures of the model fields filtered at the two interannual periods viz. 
18-36 months and 3-7 years are examined. The data in these figures represent 
model simulations from the 38-year integration with an SST-dependent EWF. 
Once again, the results have been analyzed over three zonal regions, viz. the 


169 




Time-Longitude structure at 0°: Central Pacific 
t=0.5hrs; Filtered at 48-72 months 
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Figure 9.5: Same as figure 9.1, but over the equator 

monsoon - sector (60-120° E), the central Pacific (170°-130°W) and Eastern Pacific 
( 110-50° W) and also over three latitudinal positions: equator, 10°N and 10°S. 
However, we shall present only the selected examples of simulations. 

9.1.1 The Quasi-Biennial Oscillation 

Figures 9.1 to 9.3 show the time-longitude structure of the model fields at 18-36 
months. For clarity of presentation we plot only about twenty years (100-350 
months) of total model simulation of 450 months. The most prominent feature 
in the time-longitude structures at the equator and the 10°N (figures 9.1 and 9.2 
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Time-Longitude structure at 0°: Eastern Pacific 
x=0.5 hrs; Filtered at 48-72 months 





Figure 9.6: Same as figure 9.1, but over 10°S 


respectively) is a clear eastward propagation, especially over the region 60°E and 
160°W. In the southern hemisphere at 10°S (figure 9.3), the phase propagation 
is often westward. Besides, even over the equator and 10°N, the eastward phase 
propagation is not prominent over the zonal location beyond 160 W. 

There are not many detailed analyses of the structure of QBO for comparison of 
our simulations. However, Ropelewski et al. (1992) analysed the time-longitude 
structure of QBO by filtering analysed COADS shipwinds at biennial frequency. 
The equatorial time-longitude structure of the biennially filtered data for the 
period 1956-1985 over the Indian and Pacific Oceans (40°E-80°W) clearly shows 
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Figure 9.7: Power spectra of unfiltered time series of the model u field from long 
term simulation for climatological SST dependent EWF plotted with vP(v) as the 
ordinate and period as abscissa. 

an eastward phase propagation over this region. Our model simulations agree 
well with these observed structures. More detailed observational analyses are 
required for verifying the zonal and latitudinal distributions in the structure of 
QBO implied by our model simulations. 

Another feature of the simulated time-longitude structure is the episodic nature 
of QBO activity with a time scale of about ten years. It has been recognized in 
recent years that there is considerable phase coherence among the annual cycle, 
the QBO and the low frequency variability of SO (Rasmusan, et al. 1996). The 
model simulations also support an oscillation with 8-10 years at several locations. 
The interrelationships among the annual cycle, QBO and ENSO are discussed 
later in this chapter. 
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Figure 9.8: The time series of the u-fieid, filtered at annual (solid line), QBO (dashed 
line), and ENSO (solid with symbol) scales. 


9.1.2 The ENSO Scale Oscillation 


The time-longitude structure of the model fields for the entire 450 months of 
model simulation and over the equator are presented in figures 9.4 to 9.6 for three 
latitudinal positions. To describe regional structure of the 3-7 year variability in 
our mode simulations, we have presented the results for three zonal segments. A 
prominent feature of the simulated 3-7 year oscillation over the eastern Pacific 
is its strong eastward phase propagation (figure 9.6). This eastward propagation 
is also seen, but only weakly, over the other two zonal locations viz monsoonal 
region and Central Pacific (figures 9.4 and 9.5). There are, however, several other 
characteristics that distinguish the fields at three different zonal locations from 
one another. One such characteristic is the phase relationship between zonal wind 
and evaporation. While over the monsoon region, they are almost in phase and 
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Figure 9.9: Same as figure 9.8, but for meridional wind 


are nearly out-of-phase over the Pacific regions. In particular, over the Pacific 
regions strong evaporative forcing is associated, with easterly wind anomalies. 
Figures show that the 3-7 year variability also exhibits episodic nature, more 
pronounced over the monsoon region with a low frequency variability. 


9.2 Role of Interannual Variability of SST 

Several studies suggest that the interannual oscillations are ocean-atmosphere 
coupled modes. Given the large time scale of these oscillations it is reasonable to 
expect that their evolution and structure is significantly affected by the ocean- 
atmosphere interaction. To determine the role of the interannual variability of 
SST on the genesis of QBO, we have integrated the model for 1-5 years with a 
climatological mean annual cycle of SST. The climatological mean annual cycle of 
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Figure 9.10: Same as figure 9.8, but for precipitation 


SST was constructed by averaging the monthly mean SST fields for 38 years from 
GOADS. We have examined the role of the absence of interannual variability of 
SST primarily in terms of power spectra of the model fields at different locations. 
Figure 9.7 represents the power spectra for the model u at three different zonal 
locations over the equator as indicated in the panels. The solid line and the solid 
line with symbol in figure 9.7 represent the power spectra for r = 0.5 hour and 
the corresponding red noise, respectively. As can be seen from this figure, there 
is no peak in the power spectra at the interannul periods above red noise level 
for any of the variables at any of the locations. The most prominent signal is 
centered around the annual time scale, followed by smaller but significant peaks 
at intraseasonal frequencies. Since the role of SST in our model is to modulate 
the strength of EWF, the above results imply that the interannual variation of 
convective forcing is necessary for the genesis of interannual oscillations. This in 
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turn emphasizes the importance of ocean-atmosphere interaction in interannual 
oscillations. Since the interannual variability in SST crucially dependents on the 
forcing by the surface winds. 

9.3 The Interrelationships Among the Annual 
Cycle, QBO and ENSO 

A question of great importance is the interrelationships among various oscilla- 
tions with different time scales. This is related to the phase locking of various 
oscillations. There is a strong evidence of phase locking of the annual cycle with 
the ENSO cycle (Rasmusson and Wallace 1983). Similarly, there is evidence of 
phase locking between the QBO and ENSO (Lau and Sheu 1988). Brier (1978) 
and Nicholls (1978), on the other hand, suggested a mechanism of the tropo- 
spheric QBO involving a nonlinear interaction between the annual component 
of the surface wind, surface pressure and SST. Given the fact that our model 
supports strong signals at annual, QBO and ENSO scales, we have analyzed the 
phase locking among the three oscillations by considering time series for the model 
fields filtered at the three frequencies. 

Figure 9.8 presents the time series of the model u field filtered at annual (solid 
line), QBO (dashed line) and ENSO (solid line with symbol) scales. The fields 
have been zonally averaged (before filtering) over three zonal locations as de- 
scribed before. Once again, we show the results for three latitudinal positions: 
16°N (top panels), equator (middle panels) and 16°S (bottom panels). Figures 
9.9 and 9.10 represent the corresponding results for the model v field and precip- 
itation, respectively. 

As can be seen from these figures, there exists a clear relationship between the 
QBO and the annual cycle. Not only are the two oscillations phase-locked, but 
their amplitudes vary proportionately in most of the time. In particular, they 
go through periods of intermittancy in phase. This is also largely true about the 
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ENSO signal (solid line with symbol), especially in the zonal wind (figure 9.8). A 
scrutiny, however, reveals a much richer structure in terms of zonal and meridional 
variations of the relationships among these oscillations, and highlights the need for 
more detailed analyses of observed data. For example, over the monsoon sector 
the ENSO is often of opposite phase with the QBO (and hence also with the 
annual c\ cle ) . The lack of phase locking decreases as we move to central Pacific 
where, over the ecjuatoi, there exists a much greater phase coherence. However, 
the phase locking is generally weak over the off-equatorial regions. It should 
be noted that the zonal wind, especially over the Pacific sectors, exhibit strong- 
equatorial trapping, while the signal in the meridional wind is small in general in 
all the three oscillations. While these relations between the three oscillations are 
indication of inter-period interactions, more elaborate analyses of observed as well 
as model simulations are required to build up a clear and consistent dynamical 
scenario. 


9.4 Conclusions 

The above results indicate that the convective tropical atmosphere, under- the in- 
fluence of SST, can support interannual variabilities with observed characteristics. 
Since, we have used the observed SST, this variability is already present in the 
convective forcing in our model through its dependence on SST. What the above 
results show, therefore, is that our model with convective time lag can respond 
to this forcing and support the interannual oscillations with observed periodici- 
ties. The question, of course, arises as to how the interannual oscillation arises in 
the SST, and quite likely the genesis mechanism depends on ocean-atmosphere 
coupling. An important point, however, is that once the interannual oscillation 
is generated, its structure determined by the convective dynamics is quite similar 
to the observed one. 

A number of studies have emphasised the interconnection between annual cy- 
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cle, biennial vaiiabilitv and. the SO (Rasmussm et al. 1990). In particular, the 
biennial variability appears to be fundamental to SO itself (Ropelewski et al. 
1992). A number of studies have emphasized the biennial component of the SO 
variability (Rasmussm et al. 1990; Barnett, 1991; Yasunari, 1987. 1989; Van Loo 
and Shea. 1985). On the other hand, there is also considerable evidence of a 
close relationship between the 30-50 day intraseasonal oscillation and SO both in 
terms of genesis and structure (Lau and Chau, 1986. 1988; Lau. 1985; Lau and 
Sheu. 1988). Our model with its ability to simulate this entire spectrum of SO 
is particularly suitable to investigate these aspects. Indeed, as we have seen, the 
model simulated annual cycle, QBO, and ENSO scale variabilities exhibits the 
observed phase coherence; the question of phase coherence between intraseasonal 
and interannual vraibilities needs to explained in future. 

We have seen in the context of our discussion of the total spectra of simulated 
variabilities, there also exists a significant (above red noise) decadal scale oscilla- 
tion in the model simulations at several locations. Although we have not explored 
this in the present discussion, the decadal variability is an important component 
of the tropical circulation. It is very encouraging that our model can capture this 
entire range of variabilities with a single basic mechanism. This paves the ground 
for a detailed investigation of each these oscillations under different dynamical 
and regional condtions, 
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Chapter 10 


Convective Coupling and Tropical 

Variabilities 


Recent analyses reveal that ocean-atmosphere coupling (OAC) plays a significant 
role in a wide spectrum of oceanic and atmospheric variabilities in the tropics. 
Following the conceptualization by Bjerknes (1969), it has long been held that 
the most prominent interannual oscillation (IAO), El Nino-Southern oscillation 
(ENSO), is an ocean-atmosphere coupled instability. Importance of OAC in the 
annual cycle has been shown for the eastern Pacific and the Atlantic Oceans by 
observational studies of Mitchell and Wallace (1992) and Wang (1994). Wang 
(1994) used monthly SST fields for the period 1950-1987 obtained from the Com- 
prehensive Ocean Atmosphere Data Sets (COADS) to investigate the annual cycle 
in the tropical Pacific (between 30°S-30°N and 120 o E-76°W). His analysis reveals 
two distinct components of variations in SST with respect to the equator. One 
antisymmetric, extratropical component that explains about 70% of the SST vari- 
ance on the average and a symmetric equatorial annual cycle (SEAC): The SEAC 
is found to dominate the central-eastern equatorial Pacific, with amplitude de- 
creasing with distance away from the equator. Wang (1994) concludes that while 
the antisymmetric component is essentially a delayed response to solar radiation 
with a maximum in the late summer (minimum in the late winter), the SEAC is 
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primarily governed by OAC, with maximum (minimum) in late fall (spring). In 
an earlier study Wang (1993) demonstrated the important role played by OAC 
in the annual development of the cold tongue in the eastern Pacific. In the in- 
traseasonal scale the role of OAC has been less clear, although both the tropical 
atmosphere and ocean exhibit strong signals in the intraseasonal time scale (1ST). 
The best known example of atmospheric intraseasonal oscillation (ISO) is the 40- 
60 day Madden and Julian oscillation (Madden and Julian, 1971: 1972). It is 
now known that a prominent intraseasonal signal is present also in the oceanic 
circulation. Analysis of time series data from surface moored buoys in the eastern 
equatorial Pacific (between 110 and 140°W) for the period 1975-S3 by McPhaden 
and Taft (1988 hereafter MT) revealed a very prominent ISO in zonal current, 
temperature and dynamic height. These signals were found to have character- 
istics of Kelvin waves, with eastward propagation, zero meridional velocity and 
with zonal velocity and temperature nearly in phase. Enfield (1987; E87) showed 
that the sea level fluctuations at 1ST during 1975-84 along the west coast of the 
Americas are remotely forced by atmospheric ISO in the western Pacific. These 
forced signals with Kelvin wave structure propagate eastward across the equator 
to reach the eastern boundary, where they propagate as coastal Kelvin waves. 
Although the analysis of MT revealed the period of the signals to be 60-90 days 
as against 40-60 days found by E87, it is likely that both the signals represent the 
same physical process (MT). A scenario that provides a qualitative understanding 
for such a situation is discussed later using the present model. 

There have been several attempts to develop a dynamical mechanism to explain 
the coupled variabilities in the tropics, especially the ENSO. The basic premise 
of most of these studies is that while the linear atmosphere and the ocean may be 
by themselves stable, OAC can make the (coupled) system unstable. The OAC 
in these studies involves a mechanical coupling which determines a forcing of the 
oceanic current by the (lower) atmospheric wind stress and a thermal coupling 
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which determines the effect of SS I on atmospheric heating. In most of the studies 
mentioned above, the latter is prescribed as a proportionality relation between 
SST anomaly and perturbation heating of the atmosphere. In such a formalism, 
SST acts like a direct heat source for the atmosphere, and the link between SST 
and atmospheric convection is an implicit one. Of course, the direct heat source 
provided by SST implicitely involves the convective processes, in the sense that 
the changed low-level circulation will result in modified intensity and distribution 
of convection. Although there have been a few attempts to incorporate the effects 
of convection directly into the formulation of OAC (ban and Slien, 1988; Goswami 
and Selvarajan, 1991), these formulations can be termed as direct coupling where 
SST can heat the (lower) atmosphere, regardless of presence or absence of con- 
vection. The role of convection in these formulations is indirect or implicit. One 
reason for the indirect role of convection in these studies is that the moisture vari- 
able, and hence the precipitative heating, is determined by a balance equation. 
It should be noted that although Lau and Shen (1988 hereafter LS) begin with a 
dynamic moisture equation, their subsequent assumption of a convective regime 
effectively removes the time-dependence of the moisture variable. However, in a 
dynamical scenario the convective and non-convcctive regimes themselves evolve, 


one changing to the other depending on the local and time dependent forcing of 
the moisture variable. This is particularly true for a coupled environment, where 
SST plays a crucial role in the moisture variable. Indeed, recent analyses seem 
to indicate that the atmospheric convective activity is inseparably coupled to the 


oceanic processes. 


studies point to strong association 


between SST and 


the atmospheric convective activity at different time scales. In an analysis of data 


from International Satellite Cloud Climatology Project (1S0CP), Woare (1994) 


found significant correlation of SST with cloud variables both at 1ST and intei- 
annual time scales (I AT). The analysis of Weare (1994) further shows that not 


only SST but time rate of change of SST is also important for these processes. 
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From an analysis ol intrascasonal variability of the tropical 


I > ar ilic-Imlian oroan 


and atmosphere using 7 year (1986-93) daily gridded data for variables like outgo- 
ing longwave radiation (OLR), surface stress and SST, Hendon and Chick (1995) 
report significant role of intraseasonal (30-50 day) cloud activity, and resultant 
short wave anomalies, on SST, especially over the Indian ocean eastward of about 
165°E. The importance of the evaporative process in ENSO related variabilities 


has been emphasized by Myers et al. (1986). 

The above analyses thus point to the possibility that the mechanism of OAC, 
in general, is different from the one usually assumed to explain ENSO. While 
a strong association between SST and (low-level) convection is quite expected, 
the observed association between SST and cloud activity, in particular, shows 
that the entire (lower) atmospheric column and the atmospheric convective ac- 
tivity are directly involved in OAC. It is possible that the primary mechanism 
of OAC is through the effect of SST on the atmospheric convection. This would 
imply that, unlike in the conventional formalism where the atmospheric heating 
is directly controlled by SST, the (OAC induced) atmospheric heating is possible 
only through deep heating by the column precipitation. In particular, in such a 
scenario, OAC behaves like a (second order) threshold process, somewhat akin 
to first order phase transition that results in precipitation. In the same way that 
precipitation cannot occur until the moisture content exceeds a critical value, 
the OAC induced atmospheric heating in our scenario is not effective until SST 
induced convection reaches a critical value. We shall call such a representation 


ol OAC as convective coupling, to contrast it from the conventional (direct) cou- 
pling. Thus, if the basic mechanism of direct coupling can be summarized as 
the circular interaction between (low level) winds affecting the ocean which in 
turn afFects the winds, then the basic mechanism of convective coupling can be 
represented as convective forcing affecting the winds which aftect the ocean which 
in turn affects the winds and convective forcing. 



The purpose of this study is to formulate such a convert! vo coupling an,! to 
explore the dynamics of a convectively coupled system. A noteworthy feature of 
many of the observed variabilities in the tropics is their predominantly Kelvin 
wave character. The central role played by the Kelvin waves in ENSO has been 
emphasized and utilized in many studies (Hirst, 1988; Hirst and ban, 1990; Phi- 
lander, 1990). Similarly, the prominent atmospheric ISO with 40-GO day scale is 
also known to exhibit a prominent Kelvin wave character (Madden and Julian, 
1971; Lau and Chan, 1985; Knutson and Weickman, 1987; Weickman ct al. 


1985). Importance of the Kelvin waves in the dynamics of oceanic ISO has been 
also well recognized (MT; E87; Johnson and McPhaden, 1993; Hendon and Click, 
1995). These studies, along with the results quoted above indicate that Kelvin 
wave dynamics plays a central role in the tropical variabilities at 1ST, l AT as 
well as at intermediate time scales. We shall therefore consider convectively cou- 
pled Kelvin waves to explore a unified dynamical mechanism for origin of coupled 
oscillations at different time scales. 

While we confine our dynamics to Kelvin waves, we address several questions 
pertaining to coupled variabilities in the tropics. In particular, we show that the 
convectively coupled system can support a wide spectrum of oscillations charac- 
terized by the selectively destabilized waves (the mode with the largest growth 
rate in the lower frequency regime; SDW) with diflerent periods and wavelengths. 
Thus the two prominent oscillations, viz. ENSO and ISO, with two widely dif- 
ferent times scales appear as manifestations of a single dynamical mechanism. 
However, such a unified dynamical framework lor oscillations, with periods rang- 
ing from 1ST to I AT, must take into account the differences in the strengths and 
the relative roles of various processes at the widely separated time scales. 1 hus 
we consider two scenarios: one representative of IS I and the other representative 
of I AT. These scenarios, as described below, arc characterized by the strengths of 
certain physical processes appropriate for the respective time scale. In addition 
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we also consider scenarios which are region specific, such as the Indian ocean 
summer monsoon region. It will be seen that the same mechanism that gives 
rise to I AO generates ISO as one changes from the I AT to 1ST scenario. For 
intermediate scenarios the mechanism can generate oscillations of seasonal and 
annual time scales. 

We also address the question of global nature of ENSO type of variability, as 
revealed by recent analysis (Yasunari, 1987; Nigam and Shen, 1993; Tourre and 
White, 1995). In particular, it is now recognized that ENSO type of variabilities, 
almost in phase with the Pacific events, also appear over the equatorial Indian 
ocean. It is possible that both regional mean conditions as well as basin-specific 
SST dynamics contribute to the genesis of IAO with relatively short zonal scale. 
We therefore consider two ocean models, characterized by different dynamics of 
SST. We then explore whether either of the SST models can support interannual 
variabilities with relatively shorter zonal scale. 

In section 2 we describe the two components (atmospheric and the oceanic) of 
our basic model, the concept and formulation of the convective coupling and the 
two scenarios. Section 3 discusses the dispersion relation and the eigenfunctions 
for the two ocean models. In section 4 we discuss the mechanisms of ISO and 
IAO in the two ocean models and the relative roles of various processes. Section 
5 presents our conclusions. 

10.1 The basic model 

10.1.1 The atmospheric component 

Since in our scenario, the entire lower atmosphere is affected by OA(J, we consider 
an atmospheric component of the coupled model which describes the dynamics 
of anomaly circulation of the lower troposphere. This component is represented 
by the shallow water equations describing the horizontal structure of the first 
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baroclinic mode on an equatorial /5-plane. Thus the basic set of equations Col- 
our atmospheric model is constituted by equations from 3.31 to 3.36. However, 
for discussing the coupled dynamics in an analytical setting, these equations a. re 
lineal lzed about the mean state of rest. Further, lor reasons discussed above, 


only dynamics of the Kelvin wave in a coupled environment is 

investigated. For 

the Kelvin wave, the nondimensional forms of these equations ; 

are given by, 

du„ 80 

dt ~^ + RUa ~ 0 

(10.1) 

80 

pyu a - 7T = o 

av 

(10.2) 

80 8u„ _ 

m ox + m = Q 

(10.0) 

Os _ du„ 

dt + 1 dx L 1 

(10-4) 

q = q(l+s) 

(ll).rm) 


The subscript ’a’ denotes atmospheric variables. Here, R represents both 
Rayleigh damping and Newtonian cooling coefficients which we consider to be 
equal as their difference was not found to change the results qualitatively. 

The length, time and temperature scales used for nondiinensionalizingour equa- 
tions are given by 


1 0 


c a 0 


n = L\ = On = and c? = 


r d0 

Tz 


yllOt 


L ° ~ 2(3 d c a 0o'~°~ 2/3 d 6 0 '''" " dz a 0 

The values of these quantities arc given in Table 1 A. The perturbation evaporation 

E is controlled by both atmospheric and oceanic processes. While a high SS I (7 ) 
enhances potential evaporation, ventilation by the surface wind is necessary to 
make the evaporation process effective. The contribution from the atmospheric 
process to evaporation comes essentially from evaporation-wind feedback (LW b ). 
If we assume these two processes to work in an additive manner (LS), then we 
can write, 
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E = aT + An, 


( 10 . 5 ) 

where the nondimensional forms for the evaporation-SST feedback parameter <v 
and EWF parameter A are given in the following. 

and A = 

<1 q 

The various symbols with subscript M’ used in this paper denote the corre- 
sponding dimensional value of the parameters which are represented by the sym- 
bols. The sign of A (as given in table 1) can be either positive or negative 
depending on whether the mean background winds are (implicitly) assumed to 
be westerlies or easterlies (Neelin et al. , 1987). Although the mean background 
winds in the tropics are easterlies in general, there are situations (eg. Indian 
ocean summer monsoon) when the mean winds are westerlies. 'Thus we shall 
consider either of the signs of A depending upon the situation. The values of 
A q used are adopted from earlier studies (Weare et al. , 1981 and Zhang and 
McPhaden, 1995). 

As in GR, the atmospheric heating is considered proportional to precipitation 
which in turn is proportional to s, as can be seen from eq.(10.6). Thus the 
nondimensional form of heating is given by 


Q = v*-, 


(io.fi) 


. L v 13q g, 

where ’ ” = 

The constants L v ,C p ,p a and H are described in table 1. The nondimensional 
lorm of the other parameters in the model equations are H = Il,il a, E = /u t ’and 


P = PdLoTo 
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10.1.2 The oceanic component 

The basic ocean model describes the dynamics of an oceanic mixed layer using 
reduced gravity equations for Kelvin waves. In addition we consider two equa- 
tions which represent dynamics of SST in different ocean basins, or in different 
situations. 

Ocean Model I 

In the eastern Pacific, for example, there exists an observed positive correlation 
between the thcrmoclino depth h and SST (Philande, 1990). Therefore, in our 
first SST model we take 


T = (h (10.7) 

where £ is a constant of proportionality of the order of lO~ 2 Km~ l . This is in 
accordance with the value of ( used by Hirst (1986). Then the governing equations 
for model-I can written as 

= t x ( 10 . 8 ) 

= 0 (10.9) 

= 0 ( 10 . 10 ) 


Ocean Model II 

Unlike in the Pacific basin which has a shallow thermocline, the SSI, in general, 
follows its own dynamics, although it is significantly affected by upwelling and 
downwelling. Thus in our second model SST is governed by an independent 
dynamical equation similar to one used in LS. 


where £* 


(h 

0 o 




du 0 dT 
— + £ — + «„ 

dT 

(jyu 0 + e— 

dT 'dn . 

_ + ,_ + 6k „ + ( , ; 


CK 


0q dx 
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du 0 d<j) 

dt +£ & + “"“ = 

(10.11) 

= o 

dy 

(10.12) 

~ + Gu 0 + bT = 6^ 
vt ()l 

(10.13) 

dcj> du 0 

+ ^ = 0 

(10.1 1) 


irj'i 2 

where e = -ry-, G = -p— ^ and </> = 

Ln 0n a.r 8 0 u 


Lq d'T (j)6d 


o w o 

In cc] nations 10.10 to 10.10, a, h and d represent, respectively, damping coelli- 
cients for momentum, SST and the geopotential in nondimensional form (Table 
113). For convenience, we use same 


Table IB 

Description of the coupling and scenario parameters 


Parameter 

Symbol 

Nondimensional Value 

Intraseasonal 

Interannual 

SST- 
M 1 

SST- 

M2 

SST- 

Ml 

SST- 

M2 

Evaporation- wind feedback 

A (10 _1 ) 

5.0 

-1.0 

2.5 

2.5 

Convergence feedback 

r 

1.35 

1.2 

1.2 

1.0 

Zonal gradient of mean SST 

G(10- 2 ) 

-10.5 

12.5 

-8.5 

-7.5 

Evaporation-SST feedback 

a 

4.6 

4.6 

4.6 

4,6 

Coefficient of mechanical coupling 

7 (i0- 3 ) 

9.5 

9.5 

9.5 

9.5 


SST-Ml stands for SST-model I and SST-M2 stands for SST-modcl II 


188 



symbols (e and G) to -represent the similar processes in two SST- models, though 
the values and forms of their nondimeusional counterpart differ. The variables 
?/,„and T are. perturbations in zonal current and SST respectively. The thermo- 
cline height (in geopotential units) of the upper ocean is denoted by <f>. The 
parameters appearing in equations 10.1 to 10.16 arc described in tables 1A and 
IB. 

10.1.3 Ocean- Atmosphere Interaction (Convective Cou- 
pling): 

The nature of the ocean-atmosphere coupling depends on the prescribed func- 
tional dependence of the wind stress r and the atmospheric heating Q. for the 
mechanical coupling through the wind stress we have adopted the conventional 
parameterization, t x = j d u a where •y d represents the coefficient of mechanical 
coupling whose nondimensional form is given by 7 = 7 d T 0 . 

In the convective coupling that we propose SST cannot heat the atmosphere 
directly; it influences the atmospheric heating by controlling the column moisture, 
and hence column precipitation. Further, the atmosphere is driven solely by 
precipitational heating. As such, the influence of the oceanic variables on the 
atmospheric flow would be absent in the absence of precipitative process. This 
process is included through equations 10.3, 10.-1 and 10.7. 

10.1.4 The Interannual and the Intraseasonal Scenarios: 

One of the prime objectives of our study is to explore a common dynamical 
mechanism for the coupled variabilities in the tropics. However, while the basic, 
dynamical mechanism may be the same, the oscillations with widely dillerent 
timescales, like IAO and ISO, occur over different mean states. Although we 
have no explicit reference to a mean (low, we need to incorporate this difletence 
in the mean conditions. One of the eflects of the mean conditions is that the 
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Table 1A 


Parameter 

Dimensional 
symbol value 

Nondimcnsional 
symbol Value 

Time scale 

T 0 

0.85 days 

- 

- 

Length scale 

L 0 

4400 km 

- 

- 

Temperature scale 

0o 

21 K 

- 

- 

Moisture scale 

Qo 

40-50 kg in ~ 2 

- 

- 

Mean potential temperature 

0 

310 1< 

- 

- 

Mean thermoline depth 

h 

70 in 

- 

- 

Oceanic gravity wave speed 

C 0 

1 in .s _! 

- 

- 

Tropopause Height 

Hit 

16 bn 

- 

- 

Atmospheric dry gravity 
wave speed 

Ca 

200?n.s _1 

- 

- 

Convective time lag 

T 

0.25 days 

13 

3.4 

Upwelling coefficient 

$d 

0.02° A'?n -1 

8 

9 xl0~ 1 

Oceanic Rayleigh friction 
coefficient 

(hi 

80 days" 1 

a 

l()" 2 

Oceanic Newtonian cooling 
coefficient (for 4> and T) 

bd 

200 days -1 

b 

4.3x1 O' 2 

Atmospheric Rayleigh fric- 
tion coefficient 

Rd 

1-5 day - 1 

R 

0.17-0.85 

Atmospheric Newtonian 

cooling coefficent 

Rd 

1-5 days -1 

R 

0.17-0.85 
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strengths and the relative roles of various processes, especially those which de- 
pend on the mean conditions, are likely to be quite different for widely separated 
time scales. For example, selective destabilization of waves at ISO (IAO) scale 
cannot possibly be expected to occur for a parameter strengths which are rep- 
resentative of the 1ST (I AT). It is therefore important to take this aspect into 
account. In view oi paucity of appropriate data, it is not easy to quantify these 
changes in the parameters as we change from one time scale to another. We 
shall, therefore, attempt to identify parameters that can provide a qualitative 
classification of scenarios with widely different time scales. In particular, we shall 
choose these parameters so as tp be consistent in their characteristics with the 
respective time scale. Thus, for example, a value of A that is relevant for 1ST 
should only accompany a value of G relevant for the same time scale. Following 
this reasoning, we consider two scenarios viz. interannual and intraseasonal, for 
the same dynamical system described by the above equations. These two scenar- 
ios, for each ocean model, are characterized by parameters that reflect the mean 
conditions (and hence time scales). In general, quantities averaged over IAT are 
likely to be smaller than those averaged over 1ST. Therefore, we shall assume a 
value for a scenario parameter for IAT smaller than its corresponding value for 
1ST. Specific reasons for this can vary from parameter to parameter and will be 
indicated below. Discussion on intermediate time scales can be accommodated 
by considering parameters of intermediate strengths. We choose the following 
scenario parameters: 

Strength of EWF { A): This parameter is essentially determined by the time 
mean value of the air (anemometer level)-sea humidity difference. Following the 
reasoning outlined above, we assume that the interannual scale is characterized 
by a value of A smaller than that for the 1ST. This can happen, for example, 
because the air-sea humidity difference Aq in the IAT contains an average over 
summer as well as winter conditions, and hence would be smaller than the value 
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averaged over only summer conditions. In addition, for 1ST, as in the Indian 
Ocean summer monsoon region, the mean winds can be westerly. Thus for the 
ISO, both these cases are considered. For the values of mean wind assumed in 
our model ( 5 nis ), the observed \alue ol Aq is about 3 qtti. kg~^ (Zhang and 
McPhaden. 1995) which corresponds to the value of A 0.5. We consider this value 
to characterize our intraseasonal scenario. 

Zonal gradient of mean SST (G): It is well-known that zonal gradient of mean 

SST plays an important role in coupled dynamics (LS). However, the value of G is 

likely to be different for IAT and 1ST. In addition the value of G also varies from 

basin to basin. For example, the summer monsoon Indian Ocean is characterized 

by a positive east-west SST gradient, while, in the interannual scale, the Pacific 

is characterized by a negative east-west gradient of SST (Philander. 1990). We 

(IT 

adopt a value of — — = —5.0 x 10“ ' Tm’ 1 , similar to the one used in manv 
ax 

earlier studies (Hirst, 1988; Lau and Shen, 1988), for the interannual time scale. 

Strength of convergence feedback (T): The value of convergence feedback, whose 
nondimensional strength is measured by P, varies from 1ST and IAT. If we assume 
the 1ST to represent the summer situation, the value of P for the IAT will be less 
than that for the 1ST. For a variation of q(z = 0) from 0.002 to 0.015 gmkg~ l 
eq.(10.5.b) implies a range of 0.01<P<1.6. Thus we adopt a value of P between 
1 and 1.5 with the smaller limits to characterize IAT scenario. 

It is possible to name other variables, like the mean thermocline depth (h), 
strength of the mean wind and the mean lapse rate of the atmosphere, that are 
different for different time scales and for different regions. However, we shall 
confine ourselves mainly to these three parameters to characterize the scenaiios, 
roles of certain other parameters are discussed in appropriate context. Certain 
important physical quantities as well as the parameters that do not change with 
time scale are described in table 1A while the standard values of the parameters 
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for the two scenarios are given in table IB. The question of sensitivity of the 
results to these parameters will be taken up subsequently. 

10.2 The dispersion relation and the eigenfunc- 
tions 

For both the SST models, we seek separable solutions for the coupled systems in 
the form 

sU-.M) = (10.15) 

Here £ represents any of the variables u a , v a J.s. u 0 ,T and o. The conditions for 
obtaining nontrivial solutions then result in a dispersion relation for each ocean 
model. 

(a) The dispersion relations Ocean model I: 

ia'yrjkG* — [cr s ( k 2 - a 2 R ) - i(Tk - ii\.)r)k][a a crb + keG *] = 0 (10.16) 

where G* = (iG - C*k)- 

Ocean model II: 

a , jr]k(8ujk — iGcrd ) ~ db{i<JB{k 2 — cr R ) + (Tk — iA)rjk][a a a d — k 2 e] = 0 (10.17) 

where ctr = u + iR, as = uj + iB. a a = ^ + ia. ^ ^ and a d = ^ + id. 

Equations 10.18 and 10.19 can be expressed as polynomials in u of degrees 
respectively 5 and 6. The expressions for these polynomials and the coefficients 
were arrived at using the MATLAB symbolic processor. The resultant polyno- 
mials are solved for various values of k (real). Hence, in general ui is a complex 
quantity. For further analysis of the dispersion relations we need to look at the 
structure of the eigenfunctions. 

b. The eigenfunctions 
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Since the form of the equations governing the atmospheric part do not depend 
on the particular ocean model, the functional form of the atmospheric variables 
axe also independent of the ocean model. So the eigenfunction of the zonal wind 
can be written as 

u a (. r, yj) = exp(ikx - iut - Ay 2 ), ( 10.18) 

where A = -f — • (21) 

2ctr 

\s stated earlier. is in general complex and hence the meridional structure 
of the waves are trapped only when JSe(A) > 0. The eigenfunctions for the other 
atmospheric variables are given below in terms of u a - 


n °R 

(10.22) 

5 — 7 a a 

rjk 

(10.23) 

The oceanic eigenfunctions: SST Model I 


u 0 — ; Tr^ a 

abCfa + e kG 

(10.24) 

i'yG* 

T = T7hr^ Ua 

CTb^a + e/c(jr 

(10.25) 

SST Model II: 


i'l'CTd 

U ° - 7^#?* 

(10.26) 

- iGcrd) 

t — — rrr Ua 

o h {oi<y* - k 6 ) 

(10.27) 
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( 10 . 28 ) 


• h'k 

Wa ~ kh Ua 

The physical fields are the real parts of these complex functions. 

(c) General Characteristics of the Dispersion Relation and the Eigen- 
functions: 

We have solved the dispersion relations (18) and (19) subject to the consistency 
condition: Re{ A) > 0. Before going into the detailed study of the full dispersion 
relations, it is instructive to analyze one of its limiting cases: B = 0. In our model 
the parameter B controls the atmospheric heating; thus, B = 0 corresponds to 
a case when precipitational heating is absent. As noted earlier there is no direct 
contribution of SST to atmospheric heating and SST controls the atmospheric 
heating only bv controlling the moisture availability of the atmospheric column. 
Thus in absence of precipitational heating ( B = 0), there is no contribution of the 
oceanic variables to the atmospheric heating. This is reflected in the dispersion 
relations by the fact that the terms involving ay disappear for B = 0. The 
same argument also holds for EWF, which can affect the atmospheric heating 
only through affecting the atmospheric column moisture. Thus for B = 0. all the 
terms involving EWF parameter (A) dropout of the dispersion relations. In other 
words coupling cannot give rise to ocean-atmosphere instability in the absence of 
precipitational heating. 

The other extreme case is represented by r — *• 0 which implies P * oc. 
Although it is possible to analyse the case mathematically, it does not appear 
physically relevant. In the following we shall, therefore, explore the dispersion 
relations for the existence of selectively destabilized waves which can correspond 
to observed variabilities. Although we shall demonstrate that our framework can 
support a wide spectrum of coupled variabilities, our emphasis is on interannual 
and intraseasonal oscillations. In particular, we shall look for SDW with appro- 
priate time period and zonal scale for values of the parameters that belong to 
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Figure 10.1: Dispersion curves for the standard case (refer table 1A and IB). The left 
panels show the frequency (u> r ) as a function of wave number (k). The corresponding 
growthrates (ay) as a function of k are shown in the right panels. The top and the 
bottom panels represent the intraseasonal and the interannual scenarios, respectively. 
The solid and the dashed lines represent, respectively, the SST-model I and SST-model 
II. The allowed parts of the solutions are those starting from the circle marked at the 
smaller values of k. The arrow marks on the left panels indicate the wavenumbers 
corresponding to the SDW. The coordinates are in non-dimensional units. 

either of the scenarios. 

(d) Interannual and intraseasonal oscillations: The standard case 

Out of the allowed solutions of the dispersion relation, only the ones that 
have positive growthrates are physically relevant. The growing solutions of the 
dispersion relation for the Ocean models I and II are accordingly presented in 
figure 10.1. The top panels represent the intraseasonal scenario while the bottom 
panels represent the interannual scenario. The left panels show the variation of 
the real part of uj (frequency) with wavenumber k in nondimensional units. The 
corresponding variations of the growthrate (imaginary part of cu) are shown in the 
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Intraseasonal (SST-model I) 



Figure 10.2: Longitude-latitude plots of the atmospheric and oceanic eigenfunctions 
and the evaporation field (E) for the maximally growing wave for intraseasonal scenario 
for the ocean model I. The negative contour are dashed. The coordinates are in non- 
dimensional units. 

the IAT and 1ST. The significance of this fact can be appreciated from a scrutiny 
of the energetics of the perturbations. 

The bottom right panel in each figure shows the structure of the corresponding 
evaporation field. Once again, the correlations of the evaporation field with SST. 
moisture anomaly and the zonal wind are different for different time scales and 
for the two SST models. The evaporation field is negatively correlated with 
SST and atmospheric temperature (T and 9 respectively) for the intraseasonal 
case for SST Model I and interannaul case for SST Model II (figures 10.2 and 
10.5 respectively). This feature is also present for SST-model I in the case of 
1ST scenario (figure 10.2). Such a feature in the observed field was reported by 
Cornejo- Garrido and Stone (1977) who found negative correlation between SST 
and evaporation and atmospheric heating in the regions of the Walker circulation 
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Figure 10.3: Longitude-latitude plots of the atmospheric and oceanic eigenfunctions 
and the evaporation field (E) for the maximally growing wave for interannual scenario 
for the ocean model I. The negative contour are dashed. The coordinates are in non- 
dimensional units. 


near 10°S. Low correlation between SST and rainfall (s in our model) has been 
emphasized by also Ramage (1977). Murakami and Wang (1993) reported an 
out of phase relationship between OLR and SST from an investigation of annual 
cycle of the atmospheric circulation and SST over the equatorial Pacific and 
Indian ocean. Due to lack of extensive observational analysis, it is not easy to 
verify these relations. There are, however, some evidence of such differences in 
the SST-evaporation correlations for different scenarios (Liu, 1988). 

(e) Energetics: To gain further insight into the model behavior and elucidate 
the role of various physical mechanisms, analysis of the perturbation energetics 
is carried out. The total energy of the atmospheric and oceanic part are given by 






(10.29) 


<E a >, - Dg~^-<$ s > + D;a d <sT> - A d D 2 s<SUa > 


where E a = («„ + D'g& 2 + D 2 s 2 )/2 and T] d = 

SST- Model I : 


Ly 

C p Hp a ' 


•> dT _ 

— <ry 0>. T = £7i (10.30) 

where Eq — [uq + Dj-T 2 )/ 2 

SST-model II: 


<i? 0 > ( = 7d<u a ti 0 > - Dj.—<Tu 0 > (10.31) 

where i?o = («5 + )/2 

where the angle bracket denotes average over a domain much larger than the 
disturbances. All ’D’s appeared in above equations convert the dimension of 
corresponding term into dimension of energy. Exact forms of these conversion 
factors are not of much importance as we are interested in a qualitative study. 
The energy integrals of the oceanic equations show that, as found in earlier studies 
(Yamagata, 1985), one of the necessary conditions for the instability is the net 
positive correlation between atmospheric wind and oceanic current. It can be seen 
from the figures 10.2-10.5 that this condition is satisfied for both the models and 
scenarios. However, in the present case, this condition is neither necessary nor 
sufficient as evident from equations 10.30 to 10.32. The correlation represented 
by < Tu 0 > in the case of 1ST scenario for modei-II is negative which along with 
a positive SST-gradient gives a positive contribution to the growth rate. This 
correlation is positive for model-I lor 1ST scenario which along with a ^ < 0 
gives a positive contribution. For the IAT, ^ for both the cases is negative, while 
corresponding < Tu 0 > are generally positive, once again making the contribution 
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Figure 10.4: Longitude-latitude plots of the atmospheric and oceanic eigenfunctions 
and evaporation field (E) for the maximally growing wave for intraseasonal scenario 
for the ocean model II. The negative contours are dashed. The coordinates are in 
non- dimensional units. 

positive. Another important source of energy, which is characteristic of convective 
coupling is < su a >. which is different for different scenarios. It also depends 
on the SST modeling to some extent. In the 1ST scenario, for model-I, s and u a 
are almost out of phase whereas for the model-II they are more in phase. For 
the model-II sign of A corresponding to 1ST scenario is negative. This implies a 
positive contribution to growthrate. This is also true for model-I as sign of A is 
positive. The term < sT > which measures the contribution to the growth rate 
from SST-evaporation feedback is generally positive for both the scenarios and 
models. 
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Interannual (SST-Model II) 



Figure 10.5: Longitude-latitude plots of the atmospheric and oceanic eigenfunctions 
and the evaporation field (E) for the maximally growing wave for interannual scenario 
for the ocean model II. The negative contours are dashed. The coordinates are in 
non- dimensional units. 


10.3 Relative Roles of Various Processes 


In view of the uncertainties involved in estimating many of the model parame- 
ters. it is important to address the question of sensitivity of our results to some of 
these parameters. As the number of parameters is rather large, we shall confine 
our study to those parameters which represent crucial physical processes in our 
model. It, however, may not be realistic to investigate the variation of certain 
parameters alone over wide ranges. Since these parameters also represent cer- 
tain mean conditions their variations are likely to be accompanied by changes 
in certain other parameter. Hence, the characteristics of the SDW for different 
plausible combinations of these parameters are examined. Further, these sensitiv - 
ity studies will be aimed at exploring certain specific issues of observed coupled 
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variabilities in the tropics. 

(a) The spectrum of convectively coupled instabilities: 

As mentioned above, a basic premise of this study is that it is essentially the 
same mechanism that accounts for many aspects of the OAC in the tropics. A 
large pait of the variations in the characteristics of the oscillations at different 
time scales, we expect, is due to the variation in the strengths and the, relative 
roles of various processes. Representative samples of the spectra of convectively 
coupled oscillations for the two ocean models are presented in table 2A and 2B. 
respectively for different 

Table 2A 

Spectrum of convectively coupled oscillations 
(SST-model I) 




Scenario 

Parameters 


T 

(Days) 

A 

(10 3 hm) 

Time Scales 

A 

(to- 1 

) (10- 2 ) 

r 


5.5 

-12.5 

1.35 

107 

18 

Intraseasonal 







5.0 

-10.5 

1.35 

93 

19 


5.0 

-8.5 

1.35 

175 

38 

Seasonal 







4.0 

-10.5 

1.35 

150 

34 


2.5 

-10.5 

1.10 

331 

10 • 

Annual 




337 

10 

2.0 

-8.5 

1.10 


2.5 

-8.5 

1.00 

1493 

17 

Interannual 



1.00 

1527 

1.6 

2.0 

-6.5 


scenarios. It can be seen that for both the ocean models oscillations of appropriate 
period appear in a consistent manner as the scenario changes from one time scale 
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to another. Table 2 also shows that while the characteristics of the SDW depend 
on the scenario parameters, the existence of coupled oscillation for a particular 
timescale is not critically dependent on the choice of these parameters. 

(b) Selectively destabilized oscillationsrthe governing processes: 

Several processes are known to play significant role in the excitation and the 


characteristics of the SDW. Table 3 summarizes the relative roles- of various kev 
processes in the existence of 


Table 2B 

Spectrum of convectively coupled oscillations 
(SST-model II) 



Scenario 

Parameters 



Time Scales 

A 

G 

— p- 

T 

A 


(10- 1 ) 

(10- 2 ) 


(Days) 

( 10 3 km) 


5.0 

6.5 

1.4 

89 

26 

Intraseasonal 

6.5 

8.5 

1.4 

98 

34 


3.0 

-21 

1.5 

206 

33 

Seasonal 

3.0 

-12.5 

1.3 

202 

15 


3.0 

-8.5 

1.3 

336 

20 

Annual 

2.5 

-10.5 

1.3 

386 

27 


2.5 

-5.0 

1.2 

2310 

23 

Interannual 

2.5 

-7.5 

1.2 

1172 

20 


the SDW. It can he seen that there is no SDW in absence of OAC. The existence 
of an unstable mode (without scale selection) for a = 0 is due to the fact that 
unlike in most lineai ' analyses, the linear atmospheric component of our coupled 
model can by itself support unstable waves (GR). The important role played by 
east- west gradient of mean SST is also highlighted in this table. 
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(c) Effect of mean westerlies: 

A well-known aspect of tropical circulation is the presence of strong surface 
westerlies over the monsoonal region due to the seasonal reversal of the winds. 
This reversal is also accompanied by a strong positive east-west gradient of SST 
over a large part 


Table 3 

Role of various processes in existence of selectively destabilized 

waves 

(For respective standard set of values the other parameters) 


Process 

SST- Model I 

SST- Model II 

Scenario 

— 

Scenerio 

1ST 

IAT 

1ST 

IAT 

T A 

Days 10'* km 

T A 

Days 10 3 km 

T A 

Days 10 3 1cm 

T A 

Days 10 3 km 

a = 0 

* 

k 

** 

kk 

A = 0 

35 11 

kkr 

'k 

irk 

r = o 

* 

kr 

182 41 

kk 

G = 0 

77 23 

172 1.2 

kk 

224 2.8 

5 = 0 

Not app 

dicable 

91 2.6 

4064 2.3 ■ 


* Selection at the smallest scale ** No growing modes 


of the Indian Ocean. Since the ISO signal is strongest over the Indian summer 



monsoon region, it is important to investigate the effect of these changes in the 
mean condition, on the coupled oscillations at seasonal and mtraseaslal time 
scales. Table 1 shows the characteristics of the SDW for the 1ST (and seasonal! 
scenario for mean westerlies (A < 0) and different strengths of the other scenario 
parameters G and F as well as for different strengths 


Table 4 

Effect of mean westerly on the selectively destabilized waves 


Paramet ers 

SST-model I 

SST-model II 

n 

A 

(io--) 

1’ 

G 

GO" 1 ) 

T 

(days) 

A 

(10 3 km) 

T 

(Davs) 

A 

(10 3 km) 

4.5 

-1 I 

1.1 

1.25 

212 


15 

94 

24 

4.5 

-11 

1.1 

-1.25 


irk 



** 

4.5 

- 1 1 

1.1 

12.5 

47 


22 

39 

19 

4.5 

-11 

1.1 

-12.5 


* 



irk 

4.5 

-1 1 

1.2 

1.25 

356 


10 

90 

17 

4.5 

-1 1 

1.2 

-1.25 


irk 



irk 

4.5 

-1.1. 

1.2 

12.5 

48 


16 

38 

15 

4.5 

-11 

1.2 

-12.5 


* 



irk 

4.5 

-22 

1.2 

1.25 


-k 


111 

18 

4.5 

-22 

1.2 

-1.25 


irk 



** 

4.5 

-22 

1.2 

12.5 

57 


16 

43 

15 

4.5 

_22 

1.2 

-12.5 


it 



-hk 

9.0 

-1 1 

1.2 

1.25 

226 


8.0 

68 

16 

9.0 

-11 

1.2 

-1.25 


irk 



irk 

9.0 

-1 1 

1.2 

12.5 

35 


15 

30 

15 

9.0 

- 1 1 

1.2 

-12.5 


* 



Hr 


* Selection at smallest scale ; ** growing mode absent 




of tt, while Uio oilier parameters are kept at their standard values. It can be 


seen from this table that coupled ISC) appears only when mean 
accompanied by a positive east- west gradient of mean SST. This 


westerlies are 
is consistent 


with the fact that the monsoonal Indian ocean is characterized by a positive 
gradient of mean SST. Another noteworthy feature is that SST- model I which 
is more appropriate for the eastern pacific region does not, in general, support 


modes of period less than seasonal time scale for mean westerlies, except for very 
high values of ( 5. For SST model- 1 1, on the other hand, the period of the SDW lies 
in the range of 30 110 days for a wide range of values of the scenario parameters. 
The corresponding wavelength also stays in the range of 1.5,000 to 24,000 km. 
(d) Interannual oscillations with short zonal scales: 

As mentioned above, recent analyses show iiiterammal variability also over the 
equatorial Indian Ocean, almost in 


Table 5 

I AO with short zonal scale 
(SST-Modol II) 


Parameters 

T 

(Days) 

A 

(io 3 M 

O' 

A 

(nr'L 

(i 

no- 2 ) 

9.5 

-6.5 

8.5 

1201 

8.7 

9.5 

5.0 

8.5 

1 336 

9.1 

12.5 

6.5 

6.5 

1827 

8.7 

4.5 

3.5 

5.0 

1141 

8.1 

9 

5 

8.0 

1985 

7.1 


phase with the I'll Nino over the Pacific, hut with lower amplitudes, 
interannual signals are also observed over the Atlantic Ocean (Zcbi , ”) 

hi lire conventional delay, nl oscillator mechanism, such variabilities need large 


207 



zonal extent (Hirst, 1988; Battisti and Hirst, 1989). Thus the delayed oscillator 
mechanism cannot provide a mechanism for the interannual oscillations in the 
global tropics. In order to explore whether a convectively coupled system can 
support interannual variabilities of different zonal scales, we have considered a 
number of scenarios for the interannual scale. The results are presented in table 
5. It is interesting to note that while the SST model II can give rise to IAO 
with relatively short zonal scales for a number of scenarios, such a trend is not 
available for SST model I. As emphasized earlier, the SST model II, which does 
not assume a shallow therrnocline as in SST model I is more appropriate for the 
Indian Ocean region. For SST-model 11, interannual modes with period ranging 
from about three to five years appear for a very wide range of parameters. A 
comparison of table 2B and 5 reveals that EWF and zonal gradient of moan SST 
are two primary mechanisms that give rise to these 1A0 with short zonal scales, 
(e) Role of convective time lag: 

While the above investigations provide an extensive study of sensitivity of the 
results to the model parameters, we present below a summary of a conventional 
sensitivity study of the existence and the characteristics of the SDW with respect 
to model parameters. Of particular importance is the sensitivity of the results 
to the value of CTL. All of the above studies were done for a fixed value (0.25 
day) of CTL. It is also one of the parameters that cannot be assigned a precise 
value since it embodies our hypothesis of CTL. As argued in OR, its value should 
be a few hours. It is, however, desirable that the existence and characteristics 
of the SDW for the two SST models and for different scenarios do not depend 
sensitively on the value of t. The values of period and wavelength of the SDW 
for the two scenarios and for the two SST models are presented in table (i for 
a range of CTL. We note that our results do not depend critically on the value 
of r. Indeed, for a four fold change in the value of r, from .1 to 0.4 day, the 
time periods for the SDW for the 1ST scenario remain within 115 and 66 days for 
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SST model I and between 95 and 83 days for SST model II. The corresponding 
wavelengths are between 15 and 28 thousand km. For the I AT scenario, the 


time periods of the SDW for the SST model II are between 1280 and 114-1 days, 
with corresponding wavelengths between 15000 and 22000 km. For SST model I, 
however, the changes in the time period are much larger, from 663 days to 2308 


days. Also, the trends of change in the time period as well as wavelength in the 
case of SST-rnodel II for IAT scenario are 


Table 6 

Effect of convective time lag oil characteristics of the maximally 

growing wave’s 


Convective 
time lag 

SST- Mo del I 

SST-Modoi II 

Scenario 

... _ 

Scenario 

1ST 

EAT 

1ST 

IAT 

r 

B = 

a > 

T 

T 

A 

T 

A 

T 

A 

T 

A 

0.1 

8.5 


115 

21 

663 

15 

95 

28 

1280 

22 

0.2 

4.2 


101 

20 

1093 

16 

91 

26 

1253 

21 

0.3 

2.8 


85 

18 

2271 

17 

88 

25 

1158 

18 

0.4 

2.0 


66 

15 

2302 

18 

83 

24 

| 

1 144 

15 


opposite to those for other cases. Although the reasons for this are not clear, it 
was found - that the solutions of the dispersion relations enter different regimes 
depending upon the relative strengths of the parameters. 

(f) Effect of mechanical coupling and broad spectrum ISO: 
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Since the strength of the mechanical coupling (7) is controlled by the strength 
of the mean wind, it is likely to have different strengths for different time scales or 
loi different situations. The sensitivity of the coupled oscillations to the strength 
of the mechanical coupling (or, equivalently, the product 07, since a and 7 always 
appear together in eqs. 18 and 19) is shown in table 7. The values of the mean 
winds corresponding to these values of 7 are also gi ven in the table. We assume 
that the strongest value of 7 represents, qualitatively, the 1ST scenario, with 
stiong mean wind. It can be seen that as the strength of the mean wind (coupling) 


is decreased from .10 m s 1 (7 = 19 ) to 5 


111 s 


(7 = 9.6), the period of the 


coupled ISO increases, from 44 days to 87 days for SST model I, and from 69 days 
to 90 days for SST model II. The corresponding changes in the wavelengths are, 
however, marginal, especially for the SST model 1. It is interesting to note that 


in an analysis of SST data for the period 1979-85, E87 had recorded eastward 
propagating oscillations with period 40-60 days, while MT, in their analysis for 
the period 1983-86 for the same region recorded Kelvin wave like disturbances of 
period 60-90 


Table 7 

Effect of mean wind strength and broad spectrum ISO 


u 

TO S -1 

7 

(10- 3 ) 

"SST-’M 0 d el T* 

~SSTdvTdddril^" 

T 

(Days) 

A 

(TO 3 km) 

T 

(Days) 

A 

( IQ 3 km) 

10 

19.0 

44 

11 

69 

16 

7 

13.5 

65 

13 

79 

16 

5 

9.60 

87 

15 

90 

17 


* A=0.65, r= 1 . 3 , G=-0.12, ** A=-0.13, F=1.3, G=0.12 


days. As observed by MT, this could be due to the changes in the variables due to 
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the El Nino of 1982-83. Indeed, the analysis of E87 shows that subsequent to the 
El Nino the period increased to 70 days. MT concludes that the two oscillations 

Table 8 

Sensitivity of characteristics of SDW to scenario parameters 

(1ST) 


Parameter 

SST-Model I 

SST-Model II 

Symbol 

Value 

T 

(days) 

A 

(10 3 km) 

T 

(Days) 

A 

( 1 0 3 km ) 


2.5 

29 

16 

* 

* 


3.5 

53 

16 

* 

* 

A(10- J ) 

4.5 

87 

IS 

82 

24 


5.5 

104 

19 

101 

29 


6.5 

115 

20 

116 

35 


1.0 

316 

12 

113 

23 


1.1 

106 

6.7 

103 

22 

r 

1.2 

54 

2.9 

96 

22 


1.3 

71 

15 

93 

24 


1.4 

109 

22 

90 

25 


1.5 

. 125 

27 

87 

27 


-2.0 

82 

22 

165 

30 


-4.0 

87 

21 

114 

27 

G(10~ 2 ) 

-6.0 

92 

20 

90 

26 


-8.5 

93 

19 

73 

24 


-10 

93 

18 

62 

22 


-12 

94 

18 ! 

52 

20 


* Selection at smallest scale 


observed by £87 and MTO are essentially the same phenomenon. Table 7 provides 
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a qualitative support lor such a broad spectrum ISO due to changes in the mean 
conditions, although, in practice, more than one mean process (such as zonal 
gradient of mean SST) is likely to be involved. 

(g)Effect of EWF, mean SST gradient and F 
As we have emphasized earlier, an arbitrary and independent variation of pa- 
rameters like A and G is not justified as such a variation would imply a change in 
certain mean conditions and hence corresponding changes in the other scenario 
parameters. However, for the sake of completeness, we present in Table 8 and 
Table 9 variations in the characteristics of the SDVV for both the ocean basins 
with respect to A, G and F, for the 1ST and 1AT, respectively. In each of these 
cases only one parameter was varied when all other parameters were kept con- 
stant (at their standard values). As we have seen in table 3, the existence of SDW 
critically depends on the processes of evaporation- SST feedback (a) and EWF. 
There is no SDW for a = 0, while there is SDW only at ISO for A = 0 (Table 3). 
Further, as discussed in the context of interannual and intraseasonal scenarios, A 
significantly affects the characteristics of the SDW. This trend is further reflected 
in Table 8 and Table 9. For both the scenarios, A has a profound effect on the 
existence and the nature of the SDW. The sensitivity of existence of SDW on 
A, however, is more sensitively dependent in the case I AT (Table 9), especially 
for SST- model II. The existence of IAO is also sensitive' to F, as can be seen 
from Table 9. The characteristics of IAO are, relatively, less sensitive to varia- 
tions in G. A scrutiny of tables 8 and 9 reveals that the existence of SDW with 
ISO characteristics is supported by a much wider range of scenario parameters 
in the intra-seasonal scale than in the interannual scenario. However, as we have 
mentioned earlier, conclusions from such independent variation of the scenario 
parameters may have only limited validity. 
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10.4 Conclusions 


Tho (-(Mitral emphasis of the |>r('S(Mit work has been exploration of a miili('<l dy- 
namical framework for a wide class of tropical variabilities. In our scenario, there 
is a basic dynamical mechanism which, under the interplay of various processes 


Table 9 

Sensitivity of characteristics of SDW to scenario parameters 

(IAT) 


Parameter 

SST-Modcl I 

S ST- Model 11 

Symbol 

Value 

T^ 

(days) 

A 

( 10 :i km) 

I 1 A 

(Days) ( 10' km) 


1.5 


** 

320 16 


2.5 

2221 

17 

3655 24 

A (10 -1 ) 

3.5 

366 

14 

4- 


4.5 

223 

12 



5.5 

170 

11 

* 


1.0 

2986 

18 

* 


1.1 

192 

8.0 

* 

r 

1.2 



1100 30 


1.3 


* 

893 10 


1.4 


* 

673 48 


-6.7 

537 

15 

1394 20 


-7.1 

634 

15 

1289 20 

G(10- 2 ) 

-7.5 

779 

16 

1132 19 


-7.9 

1028 

16 

1009 18 


-8.3 

1493 

17 

910 18 


* Selection at smallest scale ; ** Growing mode absent 




involved, gives rise to selective excitations at a number of observed spatial and 
temporal scales. The detailed structures of these observed oscillations, of course, 
will be enriched by various other processes, especially nonlinear processes, that 
operate beyond the genesis of the waves. 


This dynamical framework for a class of tropical variabilities discussed here is 
dynamics of tropical Kelvin wave, and its selective excitation at various scales 
through internal forcing generated by OAC and moist processes. The motiva- 
tion for exploring Kelvin waves has come from the predominantly Kelvin wave 
character of many observed tropical variabilities. The mechanism of selective ex- 
citation is through OAC and moist processes. The two distinguishing features of 
the present framework is an explicit dynamics of the moisture variable controlled 
by a convective time lag, and a convective OAC. The inclusion of an explicit 
moisture equation implies that we do not restrict our system to be necessar- 
ily in a ’convective regime’ and assumed in many earlier studies (Goswami and 
Selvarajan, 1992; LS). The inclusion of CTL embodies the hypothesis that the 
selective destabilization (or high frequency cut-off) takes place through forcing 
by convection organized at certain time scales (characterized by CTL). The basic 
premise of convective coupling is that the atmosphere docs not respond directly 
to SST heat source, but only through the deep column heating associated with 
precipitation. OAC controls this precipitational heating through the oiled of 
SST on column moisture. However, as precipitation is a threshold process, con- 
vective coupling makes OAC also a (second-order) threshold process. Although 
more complex coupled dynamics have been considered with conventional direct 
coupling (Hirst, 1986; Hirst and Lau, 1990), the present study stands apart from 
these earlier studies for its emphasis on convective coupling which also incorpo- 
rates the concept of convective time lag. Formulation of such an alternative OAC 
parameterization is motivated by the fact that the conventional formalism cannot 
account for many of the features of the observed variabilities as revealed bv recent 
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This dynamical framework for a class of tropical variabilities discussed here is 
dynamics of tropical Kelvin wave, and its selective excitation at various scales 
through internal forcing generated by OAC and moist processes. The motiva,- 
tion for exploring Kelvin waves has come from the predominantly Kelvin wave 
character of many observed tropical variabilities. The mechanism of selective ex- 
citation is through OAC and moist processes. The two distinguishing features of 
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to SST heat source, but only through the deep column heating associated with 
precipitation. OAC controls this precipitation al heating through the edect ol 
SST on column moisture. However, as precipitation is a threshold process, con- 
vective coupling makes OAC also a (second-order) threshold process. Although 
more complex coupled dynamics have been considered with conventional direct 
coupling (Hirst, 1986; Hirst and Lau, 1990), the present study stands apart from 
these earlier studies for its emphasis on convective coupling which also incorpo- 
rates the concept of convective time lag. Formulation of such an alternative OAC 
parameterization is motivated by the fact that the conventional formalism cannot 
account for many of the features of the observed variabilities as revealed by recent 
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analyses. Chief among these are the global nature of the ENSO type variabilities 
and the rather broad spectrum nature of the coupled oscillations. 

One of the chief merits of the convective coupling is that it naturally supports 
coupled modes at a number of time scales, which correspond to the spectrum of 
observed variabilities. The same mechanism that excites the intraseasonal mode 
for the 1ST scenario selectively excites the interannual mode when the scenario 
changes to IAT, as characterized by the relative strengths of the processes. This 
is in line with findings from several studies that there are marked similarities 
between ISO and IAO (Madden and Julian, 1971; ban and Chan, 1985,1986; 
ban and ban, 1986; and Barnett, 1984). The close relationship between the 
ENSO and the annual and the seasonal cycles has also been emphasized. In our 
analysis, these different variabilities emerge as SDW of time and length scales 
governed by the mean conditions. Further, these variabilities can coexist with 
their respective characteristics as in nature, as they respond to different mean 
conditions. Indeed, for the 1ST scenario, only rarely an IAO signal is excited and 
vice-versa. Another interesting fact is that while the periods of SDW changes with 
scenario parameters, these changes are in bands, such as intraseasonal, seasonal 
and interannua.1, as we change the parameter regimes. 

The present formalism also supports interannual oscillation with relatively 
small zonal scale when the SST dynamics corresponds to a, situation with a deep 
thermocline. Such a framework may have relevance in the genesis of interannual 
oscillations over the Indian Ocean almost in phase with ENSO, as shown by re- 
cent observational analysis. Since the IAO with shorter zonal scales can coexist 
with interannual oscillations of comparable periods but larger zonal scales in a.n 
ocean basin with different SST dynamics, like the ENSO with larger zonal scale 
over the Pacific, the present formalism offers a unified dynamical framework for 
genesis of coupled IAO in the global tropics. 
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Ti. is needless to emphasize that the simple' mechanism of genesis of eotiph'd 
oscillations proposed in this study cannot account for the entire structure of the 
observed spectrum. It is perhaps equally needless to emphasize the many limi- 
tations that a simple model like ours suffers from. One of the major limitations 
of the present study is that it exclusively discusses the dynamics of the coupled 
Kelvin waves. Such an approach, however, may not always work as in a. real sit- 
uation other modes may be excited which may significantly affect the dynamics. 
Such studies, with conventional formulation of OAC have been already carried 
out by several authors. However, as emphasized earlier, the focus of our study is 
the proposed parameterization of OAC, for which an investigation with the sim- 
plest but a physically relevant dynamical system appeared a logical first choice. 
The results, we feel, provide valuable insight into the relevance and effect of con- 
vective coupling in the dynamics of coupled oscillations in the tropics. This, we 
hope, will pave the way to investigate convective coupling in a realistic modeling 
exercise. 
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Chapter 11 


Concluding Remarks 


We would now like to synthesize our results from the previous chapters to arrive 
at the main conclusions regarding the role of convective forcing and convective 
time lag in the dynamics and simulation of tropical variabilities. We would also 
like to identify the major successes and failures, and the need and the scope for 
future improvements. 

The first and the main conclusion that emerges is that a significant part of the 
tropical variabilities can be understood as forced by convection induced dynam- 
ics. We have examined this aspect from a number of angles, such as comparison 
of average structure of observed and simulated anomaly fields, the anomaly an- 
nual cycle and the existence and structure of various observed oscillations. These 
multi-pronged analyses have helped us to arrive at a robust and consistent dynam- 
ical picture of the genesis and structure of the tropical variabilities. In particular, 
the same dynamical condition as well as parameter range provide the best and the 
consistent simulations of various aspects of tropical variabilities like the structure 
of the monthly anomaly fields, genesis and structure ol ISO and the structure ol 
the annual cycle. 

An important parameter in our formalism is convective time lag, which rep- 
resents a time scale for interaction between organized convection and large-scale 



dynamics. It was shown in earlier studies in analytical setting that the inclu- 
sion of convective time lag provides a high-frequency cutoff to the tropical eigen 
modes destabilized by the moist feedbacks (Closwami and Rao, 1993, (loswami 
and Mathew, 1994). Such a scenario holds also in the present case, although, 
the relevant value for the CTL in the detail and the non-linear case is about an 
hour as against a few hours implied by the analytical studies. This, we argue, is 
because the detailed structure of the variabilities requires contributions from a 
hierarchy of convective organizations, the smallest of which is characterized by a 
time scale of about 0.5 hour. As we have seen in the context of coupled dynamics 
with convective time lag, the relevant value of (I'l’b is once again a. few hours in 
an analytical (linear) setting. 

A question that naturally arises at this point is a dynamical description of the 
effect of CTL: mechanism through which it effects a high frequency cutoff. There 
is perhaps more than one way of looking at this aspect. One possible explanation 
is in terms of what the convective time lag basically implies: an adjustment 
time scale over which precipitation, and hence convective heating, takes place. 
Thus oscillations with periods much shorter than this convective adjustment time 
scale would not have time to be forced by the convective heating. Thus a high 
frequency cut-ofF naturally results. 

A significant finding from the present work has been that SST can significantly 
affect the tropical variabilities even at monthly and other intraseasonal scales. 
This has been seen in the contexts of structure of the monthly anomaly fields 
as well as the genesis and the structure of the intraseasonal oscillations. While 
SST plays a central role in most dynamical scenarios of interannual oscillations, 
the significant role of SST in the dynamics of monthly anomaly fields and iu- 
traseasonal oscillations in our formulation may appear surprising. The primary 
dynamical reason for this, we believe, is that in the present model the of[c( f, ol 
SST is essentially through its modulation of the moist feedback, and not as a cli- 
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rod, boat source. Thus SST directly controls the convective process which is the 
one of the main driving mechanisms. We have seen that this concept of direct as- 
sociation between SST and convection formulated and applied more explicitly in 
convective coupling excites the entire spectrum of oscillations under appropriate 
mean conditions. 

There are several directions in which the present model needs to be and can be 
improved. These improvements can be in terms of both dynamics and physics. 
Improvements in terms of dynamics would primarily involve a model with vertical 
structure and global (meridional) domain with more advective nonlinearity. More 
important, we feel, would be the improvements in the physics. 'I'll us a phased 
programme for further development of our model of anomaly circulation in the 
tropics would involve improving certain physics in the present first-baroclinic 
model itself. There are several such factors whose relative importance in the 
simulation of the anomaly fields needs to be investigated. 

Since in our scenario the large part of the tropical variabilities is due to the 
response of the convection driven dynamics to the prevailing mean conditions, 
a major area of improvement would be to involve the mean conditions as com- 
prehensively as possible. For example, the present simulations were generated 
for a spatially and temporally uniform value of the mean value of the column 
integrated moisture, q. Since q in our model determines the distribution and 
intensity of precipitation, and hence convective forcing, inclusion ol spatial and 
temporal (seasonal or monthly) variations of q is likely to have significant effect, 
on the model simulations. Another possibility is to use not only the mean dy- 
namical fields but also the mean precipitation field in the dynamics of anomalies 
through an appropriate parameterization. In another approach adopted by Saji 
and Goswami (1996) a considerable degree of success in simulating the anomaly 
climate in the tropis was achieved by directly relating precipitation to SSI. 

Another area where our model needs improvement was evident from the sitnu- 
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Power spectrum with period centred at 5 days: summer u-field with SST-dependent EWF 



Figure 11.1: Power spectra of zonal wind filtered at 3-7 days plotted for four values 
of CTL with vP(u) as the ordinate and period as abscissa, where P is the power and 
„ is the frequency. For each value of CTL the fields have been normalized to their 
respective maximum. 


lation of the structure of the monthly anomaly fields. It was found, in particular, 
that the simulations showed consistently high errors over the land regions. This is 
understandable, since our model has no land physics. The importance of the land 
surface processes in the dynamics of tropical circulation has been emphasized i 
many works. While our simulations clearly bring out the deficiencies brought 
about by the absence of the land surface processes, if needs to be 
and whether inclusion of a (parameterized) land suilace piocess 
would shift the model simulations closer to the observed structure. An important 
first step, of course, is to consider a realistic land-ocean distribution, as 
possible to the observed distribution. This would automatically imply certain 
change in the structure of the lower boundary foicing and, in par tic 
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Power spectrum with period centred at 5 days: summer v-field with SST-dependent EWF 
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Figure 11.2: Same as figure 11.1, but for meridional wind 

feedback. In the second phase, an appropriate formulation ol the moist, lend I rack 
over the land region has to be incorporated. 

It is worth noting that the first-baroclinic mode model developed by us can 
support a rich spectrum of physics. Indeed, the power spectrum shows, although 
not at all locations with equal strength, peaks almost at all t he observed period- 
icities relevant for the model simulations. In observation also various oscillations 
are not seen at all locations and at all seasons with equal strength. I he scenario 
that emerges, therefore, is one in which a very broad spectrum ol the observed 
tropical oscillations is excited by the convection induced internal dynamics in the 
tropics. These oscillations then respond to the seasonally and spatially varying 
mean conditions and also to ocean-atmosphere coupling with varying strengths 
and characteristics. This variation of response of the oscillations of different pe- 
riodicities to these factors then result in the observed structure of the tropical 
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Power spectrum with period centred at 5 days: summer precipitation with SST-dependent EWF 
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Figure 11.3: Same as figure 11.1, but for precipitation 

oscillations. We have seen that such a scenario is also borne out by our investi- 
gation of the convectively coupled tropical dynamics. 

Although our model can simulate a very wide spectrum of tropical oscillations 
it is still not complete. There are well-known tropical oscillations both at the high 
frequency end as well as with longer periods than what we have covered here. At 
the higher frequency end of the intraseasonal oscillations lies the s\ noptic fre- 
quencv’ easterly waves. There are, however, two distinct classes of easterly waves. 
One class is characterized by an off-equatorial ( 20°N) maximum in the amplitude 
with a zonal wavelength of about 2000 km. The other class is characterized by 
a westward phase propagating (relative to the basic current), longer wavelengths 
than the first type, equatorially trapped structure. Most of our knowledge about 
the 3-7 day westward propagating wave with a zonal scale of about 6000 km. 
comes from recent analyses of analyzed global data sets (Liebmann and Hendon. 
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Power spectrum with period centred at 5 days: winter u-field with SST-dependent EWF 
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Figure 11.4: Same as figure 11.1, but for winter months 


1990, hereafter LH90). 

The discovery of these synoptic frequency waves with relatively large zonal 
scales raises several interesting theoretical problems. These questions range from 
one of genesis of these waves, related to the general driving mechanism, to the 
variation of the characteristics of these w r aves in space and time. While theLe 
have been a large number of theoretical investigations ol the 30-o0 day eastward 
propagating wave, a very few studies have addressed the question of dynamical 
mechanism of 3-7 day wave. A successful dynamical scenario for t he genesis of 


3-7 day wave was advanced by Goswami and Goswami (1991), who showed that 
inclusion of combined effect of EWF and convergence feed track and dissipation 


selectively destabilizes the mixed Rossby gravity (MUG) wave at 3-7 day time 
scale. Since the observed spectral characteristics of the 3-7 day wave were by 


then available from the analysis of LH90, it was possible to verify the predict, ions 
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Power spectrum with period centred at 5 days: winter v-field with SST-dependent EWF 
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Figure 11.5: Same as figure 11.4, but for meridional wind 

from the theory with observations. Thus it was shown that not only the period 
but also the zonal wavelength and the e-folding length (characteristic equatorial 
trapping length) were also in close agreement with observations. The investigation 
of genesis and structure of 3-7 day wave in presence of a CTL was carried out 
in Goswami and Mathew (1994). In particular, it was shown in Goswami and 
Mathew that the same mechanism that selectively destabilizes the MRG wave at 
3-7 day time scale in presence of mean easterlies excites the 10-20 day wave in 
presence of mean westerlies. These analytical results thus provide an excellent- 
motivation for examining the existence and the characteristics of the 3-i day 
wave presence of CTL and non-linear heating. Indeed there are peaks at 3-7 day 
range in our composite one-year simulation. However, as the associated power is 
rather weak, we have not carried out a detailed analysis of this oscillation in our 
present study. It is, nonetheless, interesting to note that the power spectra of 
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Power spectrum with period centred at 5 days: winter precipitation with SST-dependent EWF 
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Figure 11.6: Same as figure 11.4, but for precipitation 


the model fields filtered at 3-7 day window exhibit many characteristics similar 
to the observed. 


Figures 11.1-11.3 show the power spectra of the model u, v and the precipitation 
fields from the composite one-year simulation, with an SST-dependent LYYK 
filtered at 3-7 days. As before, the unfiltred time series have been generated 


by averaging the fields over a zonal extent ol 60" centered at 90" E (monsoon 
sector), 80°W (Central Pacific) and 160° W (Eastern Pacific). We shall consider 


three latitudinal positions, viz. equator, 16°N (Nil) and 16°S (Sll). 1 he result! 


are shown from analysis of the fields for the summer months (.June-October) 


at different locations as indicated in the panels. The four curves in each panel 


represent the power spectra for four different values of CTL. The corresponding 
power spectra for the winter months (January- May) are shown in figures 1 1 A - 1 1 -6 
As can be seen from these figures, the model in general supports ISO with 3-7 
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day time scale, although the power is relatively weaker than that for the 30-50 
day and the 10-20 day waves in the model simulation. Since our model fields are 
normalized to its maximum value for ease of comparing results for different values 
of CTL, the smallness of power is not necessarily an indicator of weakness of the 
spectrum in the absolute sense. Figures 11.1-11.6 show that r = 0.5/m gives rise 
to the strongest signal at 3-7 days. 

The first striking feature of the power spectra of the winter fields is their dis- 
tinctly different characteristics from the corresponding summer fields. Thus, in 
contrast to the summer u, the winter u shows much weaker power over the equa- 
tor than over the off-equatorial positions. This is especially true over the eastern 
Pacific region, where the the equatorial u has very little power. This character- 
istic is in consistent with an MRG wave structure of the 3-7 day wave as found 
by Liebmann and Hendon (1990). The MRG wave character of the winter fields 
is further borne out by the power spectra of the winter v (figure 11.5), which 
shows the highest amplitude over the equator. The precipitation, however, shows 
rather distinct zonal and latitudinal characteristics from the dynamical fields (fig- 
ure 11.6). It was found that many other characteristics of the 3-7 day wave, like 
its westward propagation, are also consistently reproduced by the model simula- 
tion. But a more detailed analysis of the structure of the high-frequency part of 
our model simulation is necessary. 

At the low frequency (but less than centennial scale) end of the tropical oscil- 
lation are the decadal and inter-decadal (about 25 years) oscillations. Very little 
is known about the structure and dynamics of these variabilities. Although our 
long-period simulation is upto 38 years and hence can not be used to investigate 
interdecadal oscillations, it is interesting to note that there exists a significant 
peak above the red noise level in the 8-10 year range at several locations. Once 
again, this aspect needs to be investigated through elaborate analyses of simula- 
tions for much longer period. 
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Another aspect of the tropical oscillations that is discussed here is the question 
of their phase speeds. Although it is possible to compute a measure of the wave- 
length associated with each periodicity, and hence an associated phase speed, this 
discussion has been left out in the the present work. This and other spectral as 
well as dynamical characteristics of the simulated fields must be carried out in a 
detailed and objective manner for further evaluation of the model simulations. 

Considering the model’s present level of success and its understandable and 
conceptually consistent failures, it is a potential candidate for a tropical atmo- 
spheric model with more detailed structure like vertical variations. This would 
allow not only a discussion of the three dimensional structure of the (simulated) 
anomaly fields, but inclusion of a more detailed structure of the mean fields. The 
latter is quite important in view of the significant role of the mean circulation in 
the dynamics of the variabilities. The success ol the present formulation in the 
context of tropical variabilities indicates that the concept ol a convective time lag 
may be useful in. parameterizing convective forcing in more complex GCMs. 
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Appendix 


Description of the model parameters 


t= 7-"'= r=- = — =z =r= : = — : tzztz : 

Parameter 

symbol 

value 

Mean potential temperature 

T) 

310 K 

Tropopause Height 

II 7T 

1 5 km 

Atmospheric., dry gravity wave speed 

Ca 

200 ms- 1 

Atmospheric Rayleigh friction coefficient 

R d 

1-5 days -1 

Atmospheric Newtonian cooling coellicent 

Ri 

1-5 days -1 , 

Density of air 

T> 

1 .225 k(j m -3 
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